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SECTION  1.0 


INTRODUCTION  AND  SUMMARY 


This  report  presents  results  of  a  study  of  radiation  effects  on  electronic 
materials,  devices,  and  integrated  circuits.  Emphasis  is  placed  on  determining 
the  basic  mechanisms  of  the  interaction  of  radiation  with  these  materials  and 
components  with  a  view  toward  gaining  understanding  of  benefit  to  developers  of 
radiation-tolerant  devices.  Contract  DNA001-80-C-0146  was  a  two-year  program 
which  spanned  calendar  years  1980  and  1981.  Work  performed  during  1980  was 
documented  in  an  Interim  Report,  authored  by  0.  R.  Srour,  et  al.,  entitled 
"Investigation  of  Radiation  Effects  on  Infrared  Focal-Plane  Arrays,  VLSI  struc¬ 
tures,  Optical  Fibers,  and  Charge-Coupled  Devices."  Results  of  the  1981  work 
are  reported  in  the  present  document  (Volume  1)  and  in  a  companion  report  (Vol¬ 
ume  2).  (Volume  2  is  authored  by  A.  H.  Kalma  and  S.  Othmer.)  In  Section  2.0 
of  the  current  report,  results  of  an  investigation  of  the  mechanisms  of  neutron 
effects  on  VLSI  circuits  are  described.  Section  3.0  presents  results  of  ther¬ 
mally  stimulated  current  measurements  on  irradiated  M0S  capacitors.  In  Sec¬ 
tion  4.0,  an  investigation  of  ionizing  radiation  effects  on  HgCdTe  MIS  devices 
is  described.  The  remainder  of  the  present  section  summarizes  the  major  find¬ 
ings  reported  in  this  document. 

Results  of  analytical  and  experimental  studies  of  the  mechanisms  of  neu¬ 
tron  effects  on  VLSI  structures  are  described.  The  issue  of  whether  a  single 
neutron-produced  defect  cluster  can  produce  a  significant  amount  of  permanent 
damage  ("hard  error")  in  a  VLSI  cell  is  addressed.  Properties  of  average  defect 
clusters  in  silicon  are  calculated  as  a  function  of  neutron  energy,  including 
dimensions,  defect  densities,  and  local  electrical  properties.  The  expected 
effects  of  a  single  cluster  on  M0S  and  bipolar  VLSI  geometries  are  described, 
and  predictions  are  made  of  the  neutron  fluence  required  to  produce  a  cluster 
at  a  critical  location  in  an  integrated  circuit.  Experimental  results  obtained 
for  14-MeV  neutron-irradiated  smal 1 -geometry  bipolar  test  chips  are  described 
and  compared  with  predictions  based  on  cluster  calculations.  Much  of  the  data 
is  accounted  for  satisfactorily  in  terms  of  an  average  cluster  model.  The  pres¬ 
ent  results  suggest  that  hard  errors  may  occur  in  VLSI  structures  at  relatively 
modest  neutron  fluences.  This  situation  may  cause  lower  bounds  to  be  placed  on 
device  sizes  for  applications  requiring  radiation  tolerance. 
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In  a  related  investigation,  an  experimental  study  of  neutron-produced 
disordered  region  properties  using  a  scanning  electron  microscope  was  initi¬ 
ated.  In  this  study,  an  attempt  is  being  made  to  observe  individual  defect 
clusters  through  their  effect  on  carrier  recombination.  The  experimental 
capability  was  established  and  pre-irradiation  measurements  were  made  during 
this  program.  Measurements  on  the  test  sample  following  14-MeV  neutron  bom¬ 
bardment  remain  to  be  performed.  Experimental  techniques  and  results  of  the 
pre-i rradiation  measurements  are  described  in  this  report. 

Thermally  stimulated  current  (TSC)  measurements  have  been  performed  on 
irradiated  S i 0^  films  in  order  to  obtain  basic  mechanisms  information  regard¬ 
ing  the  nature  of  hole  traps  at  the  Sil^-Si  interface.  The  primary  goal  of 
our  work  thus  far  has  been  to  substantiate  that  the  TSC  method,  coupled  with 
C-V  measurements,  is  feasible  for  studying  the  occupancy  of  hole  traps.  This 
has  been  accomplished.  Using  a  Northrop- fabricated  capacitor  as  a  test  vehi¬ 
cle,  a  hole  trap  distribution  ranging  in  energy  from  0.7  to  1.3  eV  was  found, 
with  energy  peaks  at  0.9  and  1.1  eV  being  tentatively  assigned.  The  charge 
contributing  to  the  observed  TSC  accounted  for  most  of  the  observed  radiation- 
induced  flatband  shift.  The  Northrop  sample  also  displayed  a  slow  trapping  or 
negative-bias  instability.  The  electron  traps  responsible  for  this  phenomenon 
were  located  approximately  150  A  from  the  gate  electrode  with  an  energy  dis¬ 
tribution  ranging  from  0.9  to  1.5  eV.  The  effects  produced  by  these  traps 
were  observed  to  disappear  at  •"oom  temperature  within  two  days.  This  insta¬ 
bility  is  inconsistent  with  a  thermally  activated  process.  With  the  feasibil¬ 
ity  of  the  method  proven,  the  program  will  continue  in  1982  with  the  investi¬ 
gation  of  hoK  trapping  in  pedigreed  M0S  oxides  fabricated  by  Sandia.  These 
measurements  should  permit  correlation  of  process  variables  and  hardness  with 
hole  trap  properties.  It  is  expected  that  comparisons  will  be  made  for  the 
following  combinations  of  oxide  properties:  hard  versus  soft,  A1  vs  poly-Si 
gate,  and  wet  vs  dry  oxide  processing.  Of  major  concern  is  whether  there  are 
other  contributions  to  the  TSC  spectra,  such  as  ion  contamination  or  electron 
tunneling.  Our  measurements  will  be  correlated  with  concurrent  chemical  analy¬ 
sis  of  the  interface  to  be  undertaken  on  similar  samples  at  the  Jet  Propulsion 
Laboratory. 
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An  investigation  of  the  effects  of  ionizing  radiation  on  HgCdTe  MIS 
devices  has  been  performed.  Most  of  the  MIS  capacitors  studied  contained  ZnS 
as  the  insulator  and  had  different  native  oxide  layers  between  the  HgCdTe  and 
the  ZnS.  The  net  charge  buildup  measured  in  these  devices  following  77°K  irra¬ 
diation  is  attributed  to  the  charge  trapped  in  the  interfacial  region  (i.e, , 
between  the  HgCdTe  and  the  ZnS).  Either  electron  or  hole  trapping  may  occur 
in  that  region,  depending  on  the  sign  of  the  electric  field  applied  during 
irradiation.  This  finding  implies  that  at  least  a  fraction  of  both  the  elec¬ 
trons  and  holes  produced  by  ionization  are  mobile  in  ZnS  at  77°K.  The  amount 
of  charge  trapped  in  the  interfacial  region  depended  on  the  nature  of  that 
region,  particularly  on  the  thickness  of  the  native  oxide,  and  on  the  magni¬ 
tude  of  the  field  applied  during  irradiation.  Some  studies  of  capacitors  with 
a  CVD  Si 0^  insulator  were  also  performed.  The  results  indicated  that  both 
holes  and  electrons  are  trapped  in  the  oxide  during  77°K  irradiation,  with  the 
resulting  net  trapped  charge  being  relatively  small.  A  room- temperature  anneal 
evidently  removes  at  least  a  portion  of  the  trapped  electrons,  resulting  in 
increased  net  positive  trapped  charge. 


SECTION  2.0 


INVESTIGATION  OF  THE  MECHANISMS  OF  NEUTRON  EFFECTS  ON  VLSI  CIRCUITS 

2.1  INTRODUCTION 

Considerable  emphasis  is  being  placed  at  present  on  the  reduction  of 
integrated  circuit  geometries  to  increase  the  functional  capabilities  of  a 
single  silicon  chip.  The  issue  has  arisen  whether  the  damage  produced  by  a 
single  neutron  incident  on  a  VLSI  circuit  could  cause  a  hard  error,  i.e., 
create  enough  localized  permanent  damage  to  cause  the  failure  of  a  cell. 

Results  of  our  studies  relating  to  that  issue  are  described  in  this  section. 

In  Section  2.2,  results  of  an  analytical  and  experimental  investigation  of 
the  effects  of  14-MeV  neutrons  on  a  small -geometry  bipolar  test  chip  are  pre¬ 
sented.  Section  2.3  discusses  the  status  of  similar  experiments  to  be  per¬ 
formed  on  MOS  transistors  with  VLSI  dimensions.  In  Section  2.4,  progress  is 
presented  for  an  experiment  in  which  we  are  attempting  to  observe  the  struc¬ 
ture  of  individual  neutron- induced  defect  clusters  through  their  effect  on  the 
recombination  of  excess  carriers  produced  by  a  scanning  electron  beam. 

The  investigation  of  the  possibility  of  hard  error  production  by  neu¬ 
tron  bombardment  will  continue  at  our  laboratories  in  1982.  Much  of  the  infor¬ 
mation  presented  here  is  essentially  a  progress  report.  Our  further  findings 
will  be  described  in  the  final  report  for  Contract  DNA001-82-C-0055. 

2.2  EXPERIMENTAL  AND  ANALYTICAL  STUDIES  ON  NEUTRON- IRRADIATED  BIPOLAR 

TRANSISTORS  AND  PINCH  RESISTORS 

2.2.1  Overvi e w 

Detailed  studies  of  the  effects  of  14-MeV  neutron  bombardment  on  a  test 
chip  containing  small-geometry  bipolar  transistors  and  pinch  resistors  have 
been  performed.  Results  of  our  work  through  August,  1981  are  described  in  a 
published  paper  that  is  reproduced  here  in  Appendix  A.  The  interested  reader 
is  referred  to  that  paper  for  a  discussion  of  analytical  and  experimental  find¬ 
ings  relating  to  hard  errors.  Certain  details  were  omitted  from  the  paper  for 
brevity,  which  are  described  below.  Section  2.2.2  elaborates  on  the  experi¬ 
mental  methods  used.  In  Section  2.2.3,  the  approach  utilized  for  statistical 
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analysis  of  the  data  is  presented.  Section  2.2.4  gives  additional  comments 
regarding  the  interpretation  of  data  and  the  calculational  methods  employed. 
Our  experimental  and  analytical  plans  in  this  area  for  1982  are  given  in  Sec¬ 
tion  2.2.5. 


2.2.2  Me  a  ■> 1  jrement  Techniques 

Twenty  bipolar  test  chips  were  studied.  Each  chip  contained  five  pinch 
resistors  and,  depending  upon  yield,  up  to  eleven  transistors.  (Further  details 
regarding  these  devices  are  given  in  Appendix  A.)  The  transistors  and  resis¬ 
tors  were  sequentially  switched  into  test  circuits  under  computer  control.  The 
measurement  system  and  its  capabilities  are  described  here* 


Figure  1  shows  a  simplified  version  of  the  test  circuit  used  to  measure 
transistor  collector  current  (I  ).  A  fixed  base-emitter  voltage  (VRF)  of  720  mV 


was  applied  and  monitored  for  each  I  measurement. 


BE' 

Each  I£  value  was  adjusted 
for  variations  in  V^, and  if  varied  more  than  500  p,V  the  test  was  aborted. 
A  fixed  collector-emitter  voltage  of  1.000  volt  was  applied  and  monitored.  If 
it  varied  more  than  500  p.V  the  test  was  also  aborted.  The  collector  current 
was  converted  to  a  voltage  by  a  Keithley  Model  616  Electrometer  and  monitored 
by  a  Hewlett-Packard  Model  3456A  Digital  Voltmeter  (DVM) .  The  DVM  was  pro¬ 
grammed,  under  computer  control,  to  make  a  series  of  ten  readings  of  the  volt¬ 
age  on  the  616.  Each  of  the  ten  readings  was  averaged  over  one  hundred  power 
line  cycles  to  reduce  any  common-mode  voltages.  The  ten  readings  were  then 
averaged  and  this  average  value,  along  with  its  standard  deviation,  was  stored 
in  the  computer.  (The  percent  standard  deviations  measured  were  typically 
less  than  0.02  with  a  maximum  of  0.1?.)  Finally,  temperature  was  monitored 
for  each  Ic  measurement  and  the  Ic  value  adjusted  accordingly.  If  the  temper¬ 
ature  varied  more  than  0.1°C  from  3QUC  the  test  was  aborted. 


Figure  2  shows  a  simplified  version  of  the  test  circuit  used  to  measure 
the  pinch  resistors.  The  Hewlett-Packard  Model  3456A  DVM  was  used  in  its  re¬ 
sistance-measurement  mode  with  an  excitation  current  of  500  nA.  This  low  cur¬ 
rent  was  used  to  minimize  the  voltage  drop  across  the  pinch  resistors  (typically 
10  to  30  millivolts).  Sets  of  10  readings  were  made  in  the  same  manner  as  for 

*Certain  improvements  were  made  in  the  measurement  system  after  performing  the 
experiments  described  in  Appendix  A  which  enhanced  measurement  accuracy.  The 
description  given  here  is  for  the  final  (improved)  system. 
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Figure  1.  Simplified  diagram  of  test  set-up 
used  to  measure  transistor  collector  current. 


the  transistor  collector  currents.  Temperature  was  also  monitored  and 
corrected  for  in  a  similar  manner.  Figure  3  shows  a  detailed  wiring  diagram 
for  the  circuit  used  in  testing  both  the  transistors  and  the  pinch  resistors. 

A  simplified  version  of  the  computer  test  flow  diagram  is  shown  in  Figure  4. 

2.2.3  Methods  of  Statistical  Analysis  for  Experimental  Results 

In  our  initial  approach  to  statistical  analysis  of  the  experimental 
data  obtained  for  neutron-irradiated  bipolar  transistors  and  pinch  resistors, 
the  binomial  distribution  was  employed.  This  distribution  is  applicable  to  the 
case  in  which  only  two  outcomes  are  possible.  In  the  present  case,  the  two 
possibilities  are  that  a  cluster  is  or  is  not  observed  in  a  particular  device 
after  a  given  trial.  Neglected  in  this  treatment  is  the  possibility  of  two 
clusters  occurring  in  a  single  device  during  one  trial.  To  make  this  assump¬ 
tion  valid,  the  fluence  which  was  taken  to  correspond  to  one  trial  had  to  be 
chosen  sufficiently  small.  We  arbitrarily  chose  that  fluence  which  gave  a 
probability  of  0.1  of  yielding  a  cluster  within  one  device  in  a  single  trial, 
and  referred  to  this  as  unit  fluence.  The  possibility  of  having  two  events 
within  one  device  was  then,  correspondingly,  0.01,  which  was  neglected. 

Clearly,  as  the  unit  fluence  is  reduced  in  magnitude,  this  assumption  becomes 
much  better.  In  the  limit,  the  incidence  of  a  single  neutron  on  a  target  may 
be  considered  a  trial,  in  which  case  only  the  two  stated  outcomes  are  pos¬ 
sible.  For  any  practical  neutron  fluence,  the  number  of  trials  would  be  ex¬ 
ceedingly  large  and  application  of  the  binomial  distribution  would  be  impracti¬ 
cal.  For  this  case,  the  Poisson  distribution  applies,  since  the  probability 
of  a  "success",  i.e.,  the  presence  of  a  cluster,  is  extremely  small,  and  the 
number  of  trials  is  extremely  large.  Strictly  speaking,  then,  the  Poisson  dis¬ 
tribution  applies  rigorously,  but  the  binomial  distribution  may  be  employed 
with  little  error  as  long  as  the  "trial"  is  suitably  chosen.  We  employed  the 
Poisson  distribution  in  most  of  our  work,  but  nevertheless  continued  to  make 
use  of  the  term  unit  fluence  as  a  parameter  which  is  helpful  for  scaling,  and 
visualizing,  the  problem. 


The  Poisson  distribution  is  governed  by  a  single  parameter  A,  which  is 
closely  related  to  the  mean  of  the  distribution.  The  Poisson  probability  dis 
tribution  is  given  by 


P(m)  = 


e~V 

m ! 


(1) 
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Figure  3.  Wiring  diagram  for  the  circuit  used  in  testing  transistors  and  pinch  resistors 


1 .  ENTER  DATE  i  TIME 

2.  ADJUST  VgE  «  0.72  *D.DDD5V 

3.  ADJUST  VCE  -  I.OD  *0.0DD5V 


PROMPT  TEST  TECHNIC  I  AN  TD  INSERT 
CORRECT  OUT  INTO  TEST  SOCKET 


1.  SWITCH  APPROPRIATE  RELAYS  TD 
DESIRED  TRANSISTOR 

2.  SET  UP  DVM  TO  READ  10  SAMPLES 
EACH  DNE  MEASURED  OVER  100 
POWER  LINE  CYCLES;  AVERAGE 
THEM  E,  DETERMINE  THE  STANDARD 
DEVIATION 

3.  SWITCH  APPROPRIATE  RELAYS  TD 
MEASURE  VDr 

DC 

<t.  SWITCH  APPROPRIATE  RELAYS  TD 
MEASURE  OVEN  TEMPERATURE 

5.  CORRECT  lc  DATA  PDR  VgE  & 
TEMPERATURE 


- 1 

1.  PRINT  ERRDR  MESSAGE 

2.  ABORT  TEST 


- 1 _ 

t.  SWITCH  APPROPRIATE  RELAYS  TO 
DESIRED  PINCH  RESISTOR 

2.  SET  UP  DVM  TO  READ  10  SAMPLES 
EACH  MEASURED  DVER  100 
PDWER  LINE  CYCLES;  AVERAGE 
THEM  &  DETERMINE  THE  STANDARD 
DEVIATION 

3.  SWITCH  APPROPRIATE  RELAYS  TD 
MEASURE  OVEN  TEMPERATURE 


h.  CDRRECT  RESISTANCE  READINGS 
FOR  TEMPERATURE 


Figure  4.  Simplified  computer  test  flow  diagram  employed  in 
transistor  and  pinch  resistor  measurements. 
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where  P(m)  is  the  probability  of  having  m  clusters  in  a  given  device.  The 
parameter  X  is  equal  to  np,  where  n  is  the  number  of  trials  (neutrons/cm  in 
this  case),  and  p  is  the  probability  of  success  per  trial.  For  large  X ,  this 
distribution  approaches  the  normal  distribution,  which  is  what  we  would  expect 
to  observe  if  the  data  were  subject  only  to  a  variety  of  random  experimental 
errors. 


The  fit  to  the  pinch  resistor  data  was  the  first  to  be  obtained.  In 
this  case,  the  data  is  so  obviously  skewed  (not  normally  distributed)  that  it 
is  clearly  governed  by  events  of  low  probability  and  of  high  unit  effective¬ 
ness,  i.e.,  by  a  low  value  of  the  parameter  X.  In  order  to  make  the  problem 
tractable,  it  had  to  be  assumed  that  all  clusters  were  alike  in  their  effect 
on  the  observed  resistance.  With  this  simplification,  a  fit  to  the  data  could 
be  obtained  which  yielded  the  parameter  X  and  the  unit  effectiveness  of  the 

cluster.  The  fit  is  shown  in  Figure  8  of  Appendix  A.  The  parameters  of  the 

11  2 

fit  are  a  1  of  1.7,  a  unit  fluence  of  2  x  10  n/cm  ,  and  a  unit  effectiveness 
per  cluster  of  0.35%.  The  fit  was  made  "by  eye",  and  no  rigorous  criterion 
was  employed  to  assure  that  a  best  fit  had  been  obtained.  The  crudeness  of 
the  underlying  assumption  of  a  single  value  for  the  effectiveness  per  cluster 
did  not  justify  a  more  refined  treatment  of  the  data.  The  unit  fluence  deter¬ 
mined  from  the  fit  is  a  factor  of  two  less  than  the  expected  value.  If  this 
were  taken  at  face  value,  it  would  mean  that  the  clusters  were  somewhat  larger 
in  physical  extent  than  we  have  assumed.  However,  the  data  are  not  sufficient 
in  a  statistical  sense,  to  draw  that  conclusion  at  this  point. 

The  unit  fluence  is  determined  as  follows.  By  definition,  at  unit 

fluence  the  probability  of  having  no  clusters  is  0.9.  So  P(0)  =  0.9.  For  n 

times  the  unit  fluence,  the  probability  of  having  no  clusters  is  simply  0.9n. 

The  actually  observed  value  for  P(0),  as  taken  from  the  fit,  is  used  to  deter- 

12  2 

mine  n.  The  actual  fluence  to  which  the  fit  applies  is  3  x  10  n/cm  ,  from 
which  the  unit  fluence  is  then  determined  directly. 

The  transistor  data  present  a  sharp  contrast  to  the  highly  skewed  resis 

tor  data.  As  seen  in  Figure  7  of  Appendix  A,  the  data  exhibit  a  significant 

displacement  from  the  origin,  and  a  nearly  normal  shape.  A  Poisson  fit  to  the 

data  is  also  shown  in  the  figure  (dotted).  The  parameters  of  the  fit  are  a 

10  2 

X  of  26,  a  unit  fluence  of  1.3  x  10  n/cm  ,  and  a  unit  effectiveness  per 
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cluster  of  0.08%.  As  they  stand,  these  results  contradict  the  hypothesis  that 

the  data  are  governed  by  events  occurring  in  the  intrinsic  base,  '.nis  is  so 

because  the  observed  unit  fluence  is  a  factor  of  40  less  than  the  expected 
112 

value  of  5  x  10  n/cm  .  For  consistency,  the  clusters  would  have  to  be  much 
larger  than  we  suppose,  or  the  cluster  density  is  greater.  The  cluster  density 
is  dependent  only  on  the  neutron  cross  section,  which  is  well  known,  and  the 
cluster  size  we  assumed  is  not  likely  to  be  greatly  in  error.  It  was  subse¬ 
quently  realized,  as  rioted  in  Section  2.2.4,  that  the  transistor  data  were 
most  likely  affected  significantly  by  neutron-induced  changes  in  the  extrinsic 
base.  At  low  neutron  fluences,  this  probably  was  the  dominant  effect.  In  this 
context,  the  effective  volume  is  then  much  larger  than  that  of  the  intrinsic 
base,  and  the  unit  effectiveness  per  cluster  is,  on  average,  much  less  than 
that  expected  in  the  intrinsic  base.  The  experimental  data  are,  then,  a  com¬ 
posite  of  effects  occurring  in  the  intrinsic  and  extrinsic  base,  and  the  above 
fit  can  be  taken  as  one  limit,  in  which  the  entire  distribution  is  ascribed  to 
effects  in  the  extrinsic  base. 

An  additional  fit  to  the  transistor  data  was  obtained  under  the  assump¬ 
tion  that  the  distribution  had  been  offset  from  the  origin  as  the  result  of  a 
systematic  error.  In  retrospect,  it  is  likely  that  the  offset  was  due  to 
changes  in  the  extrinsic  base.  With  this  viewpoint,  the  fit  would  then  be 
consistent  with  the  extrinsic  base  being  responsible  simply  for  the  offset 
and  not  contributing  to  the  spread  in  the  data.  This  spread  would  be  ascribed 
entirely  to  effects  in  the  intrinsic  base.  This  assumption,  then,  represents 

the  other  limiting  case.  The  fit  is  shown  in  Figure  7  of  Appendix  A  (dashed). 

11  2 

The  parameters  of  the  fit  are  a  A.  of  3,  a  unit  fluence  of  1.1  x  10  n/cm  , 

and  a  unit  effectiveness  of  0.2%  per  cluster.  The  unit  fluence  is  again  less 

11  2 

than  the  calculated  value  of  5  x  10  n/cm  ,  but  the  difference  in  this  case 
is  not  so  great  as  to  violate  our  hypotheses.  The  indication  is  again  that 
the  clusters  are  somewhat  larger  in  physical  extent  than  assumed.  Since  we 
have  tentatively  ascribed  the  total  spread  in  the  data  to  effects  occurring  in 
the  intrinsic  base,  which  is  an  overestimate,  we  have  in  effect  determined  an 
upper  limit  to  the  mean  effectiveness  of  a  cluster  occurring  there. 

Given  the  fit  to  the  transistor  data,  it  is  possible  to  assign  upper 
limits  to  the  amount  of  damage  to  be  expected  in  a  small  device  due  to  the 
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presence  of  multiple  clusters.  It  is  apparent  from  Figure  7  of  Appendix  A 

12  2 

that  for  a  fluence  of  3  x  10  n/cm  a  significant  probability  occurs  for 
having  more  than  six  clusters  impact  a  single  critical  volume.  If  the  criti¬ 
cal  level  for  a  hard  error  in  a  device  of  such  dimensions  is  1.0%  in  collec¬ 
tor  current  change,  for  instance,  then  the  probability  for  failure  in  a  single 
device  at  such  a  fluence  is  given  by 

CD 

Xrf  P\=3^  ’  ^ 

m=5 


where  P^(m)  is  the  Poisson  probability  for  m  events,  as  determined  with  a  co¬ 
efficient  X.  In  this  particular  case,  that  probability  sums  to  0.185.  The 
parameter  X  simply  scales  with  critical  volume  and  with .fluence,  if  calcula¬ 
tions  are  to  be  made  for  other  cases.  Thus,  if  the  critical  volume  is  1  pm 
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x  1  pm  x  0.2  pm,  and  the  fluence  1  x  10  n/cm  ,  we  obtain  X  =  0.2,  for  which 
the  probability  of  obtaining  three  clusters  within  the  critical  volume  is 

3 

0.001,  which  is  a  non-negl igible  value  if  chips  contain  in  excess  of  10  de¬ 
vices.  On  average,  under  these  conditions  the  probability  of  having  a  cluster 
in  any  one  device  is  approximately  0.2. 

2.2.4  Additional  Comments 

After  performing  the  work  described  in  Appendix  A,  further  study  caused 
us  to  alter  our  viewpoint  on  several  issues.  The  first  issue,  already  raised 
in  Section  2.2.3,  is  that  measurements  of  collector  current  (Ic)  for  irradiated 
transistors  were  most  likely  influenced  not  only  by  radiation-induced  defects 
in  the  intrinsic  base  region  but  also  by  those  present  in  the  extrinsic  base. 
Figure  5  illustrates  the  structure  of  the  bipolar  test  transistors  employed  and 
shows  the  intrinsic  and  extrinsic  base  regions.  In  our  experiment,  collector 
current  was  monitored  for  an  applied  positive  emitter-base  bias.  Our  view  was 
that  Ic  would  degrade  only  due  to  neutron-produced  disordered  regions  in  the 
intrinsic  base  region.  However,  we  later  realized  that  a  portion  of  the  car¬ 
riers  injected  into  the  extrinsic  base  region  can  also  reach  the  collector  and 
contribute  to  Ic-  Defect  clusters  produced  in  the  extrinsic  base  would  then 
also  cause  the  latter  component  to  decrease.  As  discussed  in  Section  2.2.3, 
we  now  feel  that  our  experimental  data  for  irradiated  bipolar  transistors  were 
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Figure  5.  Schematic  illustration  of  the  bipolar  test 
transistor  showing  the  intrinsic  and  extrinsic  base 
regions. 

influenced  by  this  effect.  It  can  be  argued  that  the  extrinsic  component  of  Ic 
should  degrade  nonlinearly  and  reach  saturation  at  high  fluences.  Such  satura¬ 
tion  would  be  expected  when  the  diffusion  length  becomes  small  compared  to  the 

13  -2 

thickness  of  the  extrinsic  base  region.  Fluences  in  excess  of  10  cm  are 
required  to  satisfy  this  condition.  If  devices  were  irradiated  to  such  fluences 
(which  we  may  refer  to  as  preconditioning  fluences),  incremental  (lower-fluence) 
irradiations  performed  subsequently  would  be  expected  to  reveal  effects  occur¬ 
ring  in  the  intrinsic  base.  The  cluster  density  at  a  14-MeV  neutron  fluence  of 
13  2  3 

10  n/cm  is  ~one  per  p,m  ,  so  a  finite  probability  exists  that  subsequent 

incremental  irradiations  would  produce  clusters  which  overlap  with  pre-existing 

ones  within  the  intrinsic  base.  The  effects  of  such  clusters  could  well  be 

different  from  those  which  are  non-overlapping.  Basic  studies*  which  have  been 

made  specifically  to  reveal  the  effects  of  cluster  overlap  show  that  in  '^terial 

of  0.4  ohm-cm  resistivity  (relevant  to  the  present  case)  the  clusters  ir.  effect 

14  2 

act  independently  up  to  fluences  of  10  n/cm  .  Fission  neutrons  were  &  ployed 
in  that  study.  On  the  basis  of  that  work,  it  is  expected  that  the  problem  of 
cluster  overlap  is  not  a  significant  one  even  for  14-MeV  neutrons  at  fluences 
up  to  10*^  n/cm^  at  doping  concentrations  of  5  x  10*^  cm-:*  or  greater.  Thus, 
preconditioning  fluences  of  this  magnitude  are  not  expected  to  affect  results 
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adversely.  It  is  also  likely,  however,  that  useful  information  about  clusters 
in  the  intrinsic  base  can  be  obtained  even  at  fluences  less  than  those  required 
to  saturate  the  extrinsic  base  effect.  This  is  because  the  effect  of  addi¬ 
tional  clusters  in  the  extrinsic  base  asymptotically  approaches  zero  with  in¬ 
creasing  fluence.  The  variation  from  device  to  device  of  the  extrinsic  base 
effect  also  approaches  zero  well  before  complete  saturation  occurs.  Therefore, 
it  is  required  only  to  precondition  the  devices  by  irradiation  to  such  fluences 
that  subsequent  device-to-device  variations  are  less  than  the  effects  expected 
to  occur  in  the  intrinsic  base.  Measurements  performed  after  the  preconditioning 
fluences  yield  the  relevant  statistical  parameters  applicable  to  the  extrinsic 
base  region.  Subsequent  irradiations  then  yield  intrinsic  base  information 
after  the  extrinsic  base  effects  are  subtracted  out. 

Our  previous  calculations  of  neutron  damage  in  silicon  have  been  per¬ 
formed  using  an  "average  cluster"  model,  as  described  in  Appendix  A  and  in  more 
detail  in  Reference  2.  Our  viewpoint  regarding  several  aspects  of  those  calcu¬ 
lations  has  subsequently  been  modified,  and  those  modifications  are  briefly 
discussed  here. 

The  total  range  R  of  primary  recoil  atoms  was  used  to  calculate  the 
dimensions  of  average  clusters  produced  by  neutron  bombardment.  The  procedure 

3 

employed  is  discussed  in  Appendix  A  and  was  based  on  the  work  of  Sigmund  et  al . 
The  range  R  is  the  total  distance  (path  length)  that  the  recoil  travels,  and 
thus  overestimates  the  largest  dimension  of  a  cluster.  We  now  feel  that  it 
is  more  accurate,  and  more  consistent  with  Sigmund's  work,  to  use  the  pro¬ 
jected  range  Rp  instead  of  using  R  in  determining  cluster  dimensions.  Over  the 
energy  range  of  interest,  a  good  approximation  ( i . e . ,  within  about  ±7%)  is 
that  Rp  st  0.75  R.^ 

An  error  was  made  in  determining  the  minor  axis  dimension  in  our  previ¬ 
ous  calculations.  Using  Sigmund's  notation,  the  minor  axis  (b)  should  have 
been  obtained  from 


where  the  bracketed  quantities  are  mean  square  values  of  the  damage  distribution 
in  the  y  and  z  directions,  respectively.  (The  incident  direction  is  assumed  to 
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be  x  here. ) 


Assuming  cylindrical  symmetry,  (y  /p  is  equal  to 


(z2)r 


Thus, 


b  = 


(4) 


Values  for  minor  axis  obtained  using  this  expression  would  be  a  factor  of  V  2 
larger  than  those  calculated  originally  (Appendix  A),  except  that  using  Rp 
instead  of  R  in  such  determinations  will  make  the  increase  only  about  7%. 

(See  example  below.) 


Again  using  Sigmund's  notation,  the  average  damage  depth  in  the  incident 
direction  is  denoted  by  (x)^.  For  an  energetic  recoil  in  silicon,  such  as  that 
produced  by  an  incident  14-MeV  neutron,  most  of  the  damage  is  located  in  the 
last  one-half  of  the  range.  The  major  axis  of  the  region  containing  most  of 
the  damage  is  given  by  2(<Ax^)q)^^,  where 


<Ax2)d  =  (x2)d-  <xd>2  .  (5) 

Previously  (Appendix  A),  we  obtained  the  major  axis  (a)  of  an  average  defect 
cluster  using  the  relation 


a  =  R  (x)D/(x)R  ,  (6) 

where  (x)R  is  the  average  projected  range.  A  more  appropriate  expression  for 
obtaining  the  major  axis  is 

<ax2>d 

To/ 

Values  for  the  two  bracketed  ratios  in  this  equation  are  given  in  Sigmund's 

3  2 

paper.  The  use  of  Rp  instead  of  R  and  the  use  of  (Ax  as  a  measure  of 

damage  extent  causes  values  for  the  major  axis  to  be  somewhat  lower  than  those 

given  in  our  earlier  work. 

An  example  will  illustrate  the  effect  of  the  above  changes  on  average 
cluster  dimensions.  The  case  of  14-MeV  neutron  bombardment  is  considered,  and 
Table  1  lists  relevant  parameters.  Using  Eqs.  (6)  and  (7),  comparison  values 
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Table  1.  Parameters  used  to  calculate  the  dimensions  of  average 
clusters  produced  in  silicon  by  14-MeV  neutron  bombardment. 


R 

R 

P 

<Ax2)d 

<X>D 

<A 

(p.m) 

(lim) 

<*>D 

<X>R 

<X>D 

0.623 

~0.47 

0.213 

0.714 

0.057 

for  the  major  axis  were  calculated  and  are  listed  in  Table  2.  A  reduction  to 
70%  of  the  original  value  is  obtained  using  the  revised  formulation.  To  obtain 
the  minor  axis  b,  we  express  Eq.  (4)  as 


which  yields  0.227  p.m.  This  value  is  7%  larger  than  that  originally  determined. 
Table  2  also  compares  revised  values  for  cluster  volume  V  and  defect  density  N^. 
The  revised  volume  is  ~80%  of  that  originally  determined  and  the  revised  defect 
density  is  ~25%  larger  than  the  original  value. 

Table  2.  Comparison  of  original  and  revised  values  for  calcu¬ 
lated  properties  of  average  clusters  produced  by  14-MeV  neutron 
bombardment. 


a 

b 

V  7 

d  , 

(em) 

(t-m) 

(cm  3) 

(cm  3) 

ORIGINAL 

0.445 

0.212 

-14 

1.05  x  10  1 

1.81  x  1017 

REVISED 

0.310 

0.227 

8.36  x  10"15 

2.27  x  1017 

Although  the  above  modifications  make  our  average  cluster  calculations 
more  accurate,  such  a  model  is  an  oversimplification  of  the  actual  nature  of 
neutron- produced  damage.  For  example,  14-MeV  neutrons  incident  on  silicon 
will  give  rise  to  a  distribution  of  cluster  sizes  due  to  the  energy  distribu¬ 
tion  of  primary  recoils  that  are  produced.  In  subsequent  modeling,  it  may  be 
necessary  to  take  into  account  this  and  other  detailed  aspects  of  neutron  damage 
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in  order  to  successfully  account  for  experimental  data.  This  may  prove  to  be 
an  important  step  in  predicting  whether  hard  errors  will  occur  in  submicron- 
geometry  circuits. 

2.2.5  Plans 

During  1982,  we  plan  to  continue  experimental  and  analytical  studies 
on  14-MeV  neution-irradiated  bipolar  transistors  and  pinch  resistors.  Test 
samples  will  be  of  two  types:  1)  fresh  devices  of  the  same  type  as  those  ex¬ 
amined  in  1981;  2)  smal ler- geometry  bipolar  transistors.  The  latter  devices 
will  be  sought  from  an  appropriate  semiconductor  manufacturer.  It  is  highly 
desirable  that  device  dimensions  be  ^1  pm  to  improve  sensitivity  to  the 
presence  of  defect  clusters  in  critical  regions. 

The  experimental  results  described  in  Appendix  A  were  obtained  based 
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on  a  single  irradiation  to  3  x  10  n/cm  .  In  1982,  we  plan  to  irradiate  a 
device  population  several  times  at  lower  fluences  and  measure  changes  in  elec¬ 
trical  properties  after  each  incremental  irradiation.  This  procedure  will 
allow  the  degradation  process  to  be  observed  in  more  detail.  Considerable 
information  regarding  the  statistics  of  neutron  damage  in  silicon  should  also 
result.  In  addition,  we  plan  to  perform  irradiations  at  fluences  appropriate 
for  separating  extrinsic  and  intrinsic  base  region  effects,  as  described  in 
Section  2.2.4. 

Modeling  of  neutron  damage  in  silicon  will  also  continue.  We  plan  to 
modify  our  previous  calculations  in  1982  by  incorporating  changes  such  as 
those  discussed  in  Section  2.2.4.  More  refined  calculations  may  also  be  de¬ 
sirable,  depending  on  the  results  of  the  experimental  program. 

2.3  INVESTIGATION  OF  HARD  ERRORS  IN  MOS/VISI  CIRCUITS 

The  potential  effects  of  a  single  neutron-produced  disordered  region 
on  the  electrical  properties  of  smal 1- geometry  MOS  transistors  are  described 
in  the  paper  contained  in  Appendix  A  'f  this  report.  Anticipated  effects 
include:  channel  mobility  degradation,  ncreased  drain  leakage  current,  and 
a  shift  in  threshold  voltage.  During  th i 3  program,  we  were  successful  in 
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obtaining  MOS  devices  with  dimensions  small  enough  such  that  these  effects 
should  be  observable. 

A  variety  of  MOS  transistors  contained  on  several  types  of  test  chips 
were  obtained:  n-  and  p-channel  transistors;  enhancement  and  depletion  mode; 
bulk  and  SOS.  Devices  with  the  following  channel  widths  and  lengths  are  on 
hand:  0.5  x  0.5,  1  x  0.5,  1.5  x  0.5,  2  x  0.5,  3  x  0.5,  4  x  0.5,  lxl, 

1.5  x  1,  2  x  1,  4  x  1,  1.5  x  1.5,  and  2.2  x  1.5  p,m.  (Several  larger-sized 
devices  are  also  available  on  the  same  chips  which  are  of  lesser  importance 
to  the  current  study.) 

These  small -geometry  MOS  transistors  were  obtained  late  in  the  present 
program,  so  the  experimental  study  of  the  effects  of  14-MeV  neutron  bombard¬ 
ment  will  be  performed  during  1982.  Results  of  that  investigation  will  be 
reported  in  the  Final  Report  for  Contract  DNA001-82-C-0055. 
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2.4  EXPERIMENTAL  STUDY  OF  DISORDERED  REGION  PROPERTIES  USING  A  SCANNING 

ELECTRON  MICROSCOPE 

2.4.1  Intro  auction 

The  structure  of  neutron-produced  defect  clusters  on  a  microscopic 
scale  is  of  importance  with  regard  to  the  likelihood,  and  the  severity,  of 
hard  errors  in  smal 1- geometry  devices.  The  vast  majority  of  experimental  evi¬ 
dence  relating  to  clusters  is  indirect,  however,  and  reliance  must  be  placed 
on  fairly  simple  models  to  yield  insight  into  the  physical  structure  of  the 
cluster  itself.  Of  principal  concern  in  the  present  context  is  the  physical 
extent  of  the  cluster  and  the  statistical  variation  in  cluster  size  for  sili¬ 
con.  Recent  calculations  of  defect  configurations  using  the  MARLOWE  code  have 
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yielded  pictures  of  defect  clusters  dominated  by  extended  damage  tracks. 

Since  there  is  little  uncertainty  about  the  average  number  of  defects  created 
by  a  silicon  knock-on  of  a  given  energy,  those  calculations  imply  simply  that 
the  given  number  of  defects  is  distributed  over  a  larger  volume  in  the  crystal 
than  would  otherwise  be  the  case,  and  thus  a  lesser  likelihood  prevails  that 
all  of  the  defects  associated  with  a  cluster  would  fall  within  the  active 
region  of  a  device.  Only  a  small  change  in  the  large- angle  scattering  proba¬ 
bility  assumed  in  the  code  would  alter  the  cluster  shape  considerably,  how¬ 
ever,  and  it  is  unclear  just  how  much  credence  to  give  to  the  calculated 
results. 


Given  the  prevailing  uncertainty  regarding  microscopic  cluster  models, 
two  approaches  were  pursued  to  gain  further  insight.  The  first,  discussed  in 
previous  sections  and  in  Appendix  A,  relies  on  statistical  analysis  of  the 
event  rate  unoer  neutron  irradiation  to  yield  information  regarding  the  maxi¬ 
mum  physical  extent  of  the  cluster.  The  second  approach,  treated  below,  is  an 
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attempt  to  observe  clusters  directly.  In  our  previous  work,  such  an  attempt 
was  made  by  examining  the  contribution  of  individual  defect  clusters  to  gener¬ 
ation  current.  High  precision  measurements  of  capacitance  transient  recovery 
from  deep  depletion  were  made  on  MOS  capacitors  and  an.lyzed  to  yield  local 
generation  rate  with  high  spatial  resolution.  The  experimental  sensitivity 
obtained  with  the  available  samples  was  such  that  we  tell  somewhat  short  of 
being  able  to  observe  individual  clusters.  In  view  of  the  manifest  difficulty 
of  that  approach,  an  attempt  was  then  initiated  to  observe  individual  clusters 
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through  their  effect  on  carrier  recombination.  The  experimental  capability 
was  established  and  pre-i rradiation  measurements  were  made  during  this  pro¬ 
gram.  The  test  sample  will  be  irradiated  with  14-MeV  neutrons  along  with  the 
small -geometry  devices  employed  in  the  statistical  study  of  cluster  proper¬ 
ties.  Experimental  techniques  and  results  of  the  pre-irradiation  measurements 
are  described  here. 

2.4.2  Experimental  Method 

The  experimental  approach  employed  in  the  investigation  of  cluster  prop¬ 
erties  was  based  on  the  use  of  a  scanning  electron  microscope  to  yield  a  highly 
localized  source  of  excess  carriers.  The  variation  in  electron-beam-induced 
current  (EBIC)  observed  as  the  beam  is  scanned  was  expected  to  reveal  the  in¬ 
fluence  of  defect  clusters.  The  selected  test  vehicle  was  a  space-quality 
thin  silicon  solar  cell.  The  experimental  configuration  is  shown  in  Figure  6. 

A  low-energy  electron  beam  is  shown  impinging  on  the  back  surface  of  a  3-mil 
soTar  cell.  The  back-surface  field  serves  to  reflect  the  minority  carriers 
diffusing  to  the  rear  surface.  The  fraction  of  carriers  available  for  collec¬ 
tion  at  the  front  surface  is  deemed  to  be  impacted  directly  by  clusters  falling 
within  the  energy  deposition  volume  of  the  beam,  whereas  clusters  distributed 
in  the  bulk  would  yield  only  an  average  reduction  of  collected  current  governed 
by  the  bulk  diffusion  length.  Thus,  if  observation  is  restricted  to  the  spa¬ 
tial  scale  of  the  clusters  (0.1  urn),  then  sensitivity  is  restricted  to  clusters 
at  the  back  surface.  High  spatial  resolution  in  terms  of  depth  is  obtained  by 
slight  variation  of  the  electron  beam  energy. 

The  motivation  for  the  observation  of  clusters  using  the  recombination 
process  rather  than  the  generation  process  was  principally  that  defect  clusters 
are  more  prominent  in  recombination.  This  is  due  to  the  presence  of  the  space 
charge  region  around  the  cluster,  which  enhances  minority  carrier  capture,  and 
to  the  fact  that  essentially  all  of  the  deep  levels  created  by  neutron  damage 
are  active  in  recombination,  whereas  generation  is  sensitive  preferentially  to 
the  subset  of  levels  which  is  very  near  midgap.  Since  we  are  attempting  to 
observe  on  the  order  of  200  defects  per  cluster  in  total,  observing  the  defects 
via  recombination  promised  greater  prospects  of  success.  A  general  motivation 
for  this  undertaking  was  also  the  desire  to  extend  existing  diagnostic  capabil¬ 
ities  to  the  spatial  scale  appropriate  to  current  device  technologies. 
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Figure  6.  Diagram  of  experimental  method  for  examining  cluster 
dimensions.  Rear-surface  beam  injection  into  a  thin,  space- 
quality  Si  solar  cell  is  employed  for  localized  generation  of 
excess  carriers. 

Advantage  was  taken  of  the  fact  that  a  new  experimental  capability 
utilizing  a  scanning  electron  microscope  was  recently  implemented  at  our  labo¬ 
ratory.^  By  use  of  a  modulated  beam  in  conjunction  with  pnase-sensitive  de¬ 
tection,  sensitivity  to  extremely  small  (pA-level)  electron-beam-induced  cur¬ 
rents  can  be  obtained.  With  this  sensitivity,  in  conjunction  with  signal 
averaging  techniques,  it  was  anticipated  that  the  recombination  current  asso¬ 
ciated  with  an  individual  defect  cluster  could  be  resolved. 

The  experimental  technique  yields  spatial  resolution  which  is  a  function 
of  beam  energy.  Energy  deposition  profiles  for  the  low  energies  employed  are 
reproduced  in  Figure  7.  They  are  based  on  the  work  of  Everhart  and  Hoff/ 
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Figure  7.  Energy  deposition  profiles  in  silicon  for 
a  number  of  moderate  electron  beam  energies,  based 
on  the  work  of  Everhart  and  Hoff. 7 

The  lateral  resolution  is  on  the  order  of  one-half  of  the  extrapolated  range 
at  a  given  energy.  Thus,  for  the  case  of  6  keV,  the  extrapolated  range  is 
0.4  pm,  and  the  lateral  resolution  is  0.2  pm.  The  lateral  dose  has  been  calcu 
lated  using  Monte-Carlo  methods  by  several  researchers.  Results  of  such  calcu 
lations  by  Chadsey  for  the  specific  case  of  lateral  dose  distribution  at  a 
depth  of  0.6  pm  for  a  10- keV  beam  are  reproduced  in  Figure  8.  The  dose  at 
distance  r  from  the  axis,  expressed  as  a  ratio  with  respect  to  the  dose  on- 
axis,  is  down  to  \%  at  about  50%  of  the  range  at  10  keV,  which  is  1  pm.  At 
lower  beam  energies,  the  lateral  beam  spread  is  correspondingly  smaller.  This 
criterion  sets  an  upper  limit  on  the  usable  beam  energy.  Depth  resolution  is 
also  enhanced  at  lower  beam  energies.  By  operating  over  a  range  of  energies 
from  4-6  keV,  depth  resolution  of  less  than  0.1  pm  can  be  achieved.  By  use  of 
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Figure  8.  Lateral  dose  profile  at  10  keV  in  silicon  at 
a  depth  of  0.6  mti.  °  Shown  is  the  ratio  of  the  dose  at 
distance  r  from  the  beam  axis  relative  to  the  dose  on 
the  axis. 

phase-sensitive  detection,  beam-induced  current  sensitivities  of  200  fA  have 
been  demonstrated  with  this  experimental  technique. 

Using  the  simplest  assumptions  about  neutron-induced  clusters  which 
were  capable  of  yielding  to  solution,  the  effect  of  an  individual  cluster  on 
recombination  current  was  calculated  to  be  on  the  order  of  0.1/;'  of  the  total 
EBIC  current  for  an  assumed  beam  current  of  2  pA.  It  was  therefore  necessary 
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to  reduce  the  noise  level  of  the  measurement  technique  below  this  level. 
Furthermore,  the  inherent  variation  of  EBIC  current  with  beam  position  also 
needed  to  be  small  compared  to  this  value  on  the  spatial  scale  of  0.1  pm. 

2.4.3  Exper i men  ta 1  Resu 1 ts 

Test  solar  cells  were  fabricated  for  us  by  Applied  Solar  Energy,  Inc. 

The  back-surface  metallization  was  removed  by  etching  over  the  central  area 
of  the  cell  for  beam  access.  Results  of  EBIC  measurements  for  a  lateral  scan 
across  the  etched  area  revealed  large  variations  in  signal,  as  shown  in  Fig¬ 
ure  9a.  The  beam  voltage  used  for  these  measurements  was  11  kV,  and  a  cur¬ 
rent  gain  (EBlC/beam  current)  of  about  1000  was  observed.  The  variation  in 
EBIC  was  about  9%  peak-to-peak.  The  variation  in  signal  amplitude  was  found 
to  be  due  to  the  fact  that  the  silicon  back  surface  had  been  partially  etched 
as  well,  when  the  metal  was  removed,  so  that  the  back-surface  field  was  par¬ 
tially  eliminated  in  certain  areas.  This  was  demonstrated  by  measuring  the 
EBIC  through  the  metallization  at  a  beam  voltage  of  36  kV,  which  is  sufficient 
to  deposit  some  energy  beneath  the  metallization.  The  EBIC  variation  was  less 
than  \%  on  the  same  spatial  scale,  as  shown  in  Figure  9b.  Three  traces  are 
shown  on  the  figure,  showing  the  normalized  response  after  respective  summa¬ 
tions  of  one,  two,  and  four  sweeps  of  the  sample.  The  traces  are  shifted  arbi¬ 
trarily  along  the  vertical  axis  for  convenience  of  comparison.  The  enhancement 
in  sensitivity  afforded  by  signal  averaging  is  apparent,  and  the  prospect  of 
observing  information  on  an  amplitude  scale  of  0.1$  seemed  reasonable.  Unfortu¬ 
nately,  the  same  figure  reveals  a  systematic  drift  of  the  beam  which  was  deter¬ 
mined  to  be  20  nm/min,  which  severely  limits  the  amount  of  averaging  that  can 
be  performed  if  information  on  the  submicron  spatial  scale  is  to  be  observable. 
(One  sweep  requires  50  sec.)  The  obliteration  of  the  signal  is  in  fact  quite 
apparent  in  the  figure,  as  contrast  is  seen  to  be  lost  from  the  top  to  the 
bottom  trace  in  the  feature  being  followed  as  a  reference  point. 

A  second  set  of  solar  cell  specimens  was  fabricated  by  Applied  Solar 

Energy  in  which  the  back  metallization  was  not  put  over  the  entire  cell.  Frag- 
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ments  of  about  1-2  mm  were  then  cut  out  (using  a  diamond  saw)  which  encom¬ 
passed  both  metallized  and  unmetallized  regions.  The  use  of  small  samples  was 
nece-  -’*.ated  by  the  desirability  of  small  device  capacitance,  since  beam  modu¬ 
lation  at  200  kHz  was  employed  and  the  junction  capacitance  sets  the  impedance 
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Figure  9.  a)  EBIC  signal  for  a  lateral  scan  across  a  solar  cell. 
The  beam  voltage  was  11  kV  and  the  current  gain  was  about  1000. 
b)  EBIC  signal  for  a  lateral  scan  across  the  metallized  rear  sur¬ 
face  of  a  solar  cell.  Shown  are  the  normalized  results  for  one, 
two,  and  four  summations  of  repeated  sweeps  (arbitrarily  translated 
along  the  vertical  axis).  Some  lateral  drift  in  beam  position  is 
observed  (20  nm/min).  The  beam  voltage  was  36  kV  and  current  gain 
was  about  20. 


scale.  The  bulk  resistivity  of  these  solar  cells  was  0.1  ohm-cm,  which  was 
selected  on  the  grounds  that  the  space-charge  region  width  surrounding  the 
cluster  should  be  comparatively  small.  Assuming  a  barrier  potential  of  0.1  eV 
for  the  cluster,  the  depletion  region  width  for  this  resistivity  (p-type  mate- 
rial,  5  x  10  cm  doping  concentration)  is  on  the  order  of  200  A,  which  is 
small  compared  to  the  assumed  cluster  dimensions. 

Results  obtained  with  these  samples  before  irradiation  are  as  follows. 
The  amplitude  resolution  which  could  be  obtained  with  summation  of  three  suc¬ 
cessive  traces  was  0.06%.  Conditions  were  therefore  marginal  for  the  observa¬ 
tion  of  clusters  after  irradiation,  since  some  degradation  in  signal-to-noise 
ratio  had  to  be  anticipated.  A  penalty  in  the  use  of  low-resistivity  material 
is  that  the  diffusion  length  degrades  more  readily  than  in  higher  resistivity 
material.  In  addition,  the  pre-irradiation  diffusion  length  is  lower  as  well. 

In  the  course  of  the  EBIC  measurements  on  these  samples,  a  slow  drift 
in  gain  was  observed  in  those  regions  which  received  a  significant  beam  expo¬ 
sure.  Sample  results  are  shown  in  Figure  10,  where  the  beam  was  scanned  past 
a  surface  feature  which  served  as  a  fiducial  mark.  The  EBIC  signal  is  seen  to 
decrease  near  the  left  periphery  (labeled  "a")  which  represents  the  boundary 
of  the  field  of  exposure.  A  second  strong  decrease  is  observed  within  that 
field  (labeled  "b"),  which  resulted  from  a  period  of  exposure  of  the  sample  at 
higher  magnification.  Such  a  scan  was  employed  to  adjust  for  critical  beam 
focus.  The  reduction  in  EBIC  as  a  result  of  beam  exposure  is  seen  to  be  on  the 
order  of  4%  for  typical  exposures  encountered  in  a  test.  No  determination  was 
made  as  to  whether  the  effect  was  reversible  or  permanent.  It  is  conceivable 
that  it  may  be  due  to  charge  buildup  in  the  native  oxide  layer  at  the  back 
surface,  which  may  result  in  surface  inversion  and  thus  parasitic  carrier  col¬ 
lection.  It  would  seem  to  be  quite  difficult,  however,  to  achieve  such  inver¬ 
sion  of  the  heavily  doped  layer  at  the  back  surface.  Alternatively,  we  may  be 
observing  the  effects  of  defect  migration  into  the  bulk  under  the  influence  of 

the  electron  beam.  Defects  have  been  observed  in  germanium  following  50-kV 
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electron  bombardment.  The  permanence  of  the  effect  observed  would  allow  a 
choice  to  be  made  between  these  possibilities. 
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Figure  10.  EBIC  signal  for  a  lateral  scan  past  a  surface 
feature  serving  as  a  fiducial  mark.  The  decrease  in  EBIC 
signal  at  points  (a)  and  (b)  appears  to  be  related  to  cumu¬ 
lative  beam  exposure,  since  the  region  of  decrease  coincides 
with  the  area  scanned.  The  decrease  at  (b)  is  ascribed  to 
previous  scans  at  higher  magnification. 

The  dependence  of  current  gain  on  beam  energy  was  determined  for  our 
samples  over  the  range  of  4  to  22  keV.  Results  are  shown  in  Figure  11.  Using 
a  formalism  described  el sewhere^ ’ these  data  can  be  interpreted  in  terms  of 
a  local  carrier  collection  probability  if  depth-dose  curves  such  as  those  of 
Figure  7  are  assumed.  Using  this  formalism,  a  carrier  collection  probability 
of  21%  at  a  depth  of  70  p,m,  just  beneath  the  back-surface  field  region,  was 
determined.  The  finite  collection  efficiency  was  ascribed  to  the  finite  bulk 
diffusion  length  in  this  low-resistivity  material.  Since  the  collection  effi¬ 
ciency  was  determined  using  ac  modulation,  a  corresponding  analysis  of  the 
diffusion  process  was  required  to  determine  the  diffusion  length. 

A  one-dimensional  analysis  of  transport  is  sufficient  to  the  purpose, 
and  this  has  been  previously  performed.^’ ^  Under  sinusoidal  excitation  con¬ 
ditions,  the  diffusion  length  is  foreshortened,  and  the  apparent  diffusion 
length  1(a)  is  given  by  the  relationship 
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Figure  11.  Dependence  of  current  gain  UrBjC/In  .  ). 
normalized  to  beam  energy,  on  beam  energy.  e  1 


where  t  is  the  carrier  lifetime  and  a  the  angular  modulation  frequency.  The 
experiment  yields  L(a)  directly,  and  the  dc  diffusion  length  L ( 0 )  can  be  deter¬ 
mined  iteratizely  usirg  the  above  relationship,  if  a  bulk  diffusivity  D  is 
assumed,  together  with  the  relationship  L  =v'  Dt.  By  this  method,  an  apparent 
diffusion  length  1(a)  of  44  ^m  was  determined,  with  the  corresponding  true 
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bulk  diffusion  length  being  52  *.m.  The  diffusivity  was  taken  to  be  11.8  cm  / 
sec. 

On  the  basis  of  these  results,  an  appropriate  neutron  fluence  was 
selected.  A  tradeoff  was  required  between  the  desire  on  the  one  hand  to  have 
a  reasonably  high  density  of  clusters  to  examine,  and  the  necessity  on  the 


other  hand  not  to  degrade  the  diffusion  length  excessively  because  of  the 
resulting  loss  in  signal-to-noise  ratio.  For  purposes  of  this  calculation, 

the  lifetime  damage  constant  was  assumed  to  be  1  x  10  n-sec/cm  .  This  value 

12  112 
was  extrapolated  from  higher  resistivity  data.  A  fluence  of  10  n/cm 

was  expected  to  yield  a  reduction  in  signal  strength  of  a  factor  of  four.  The 

results  of  post-irradiation  testing  on  this  sample  will  be  presented  in  the 

final  report  for  Contract  DNA001-82-C-0055. 


SECTION  3.0 


THERMALLY  STIMULATED  CURRENT  MEASUREMENTS  ON  IRRADIATED  MOS  CAPACITORS 

3.1  INTRODUCTION 

The  susceptibility  of  silicon  MOS  devices  to  ionizing  radiation  is  an 
important  problem  that  is  being  addressed  at  various  research  laboratories  and 
semiconductor  fabrication  facilities.  The  purpose  of  the  present  investigation 
is  to  provide  basic  mechanisms  information  regarding  the  nature  of  hole  traps 
which  have  been  found  to  exist  preferentially  at  the  Si-SiO^  interface  in  irra¬ 
diated  MOS  devices  so  that  hardening  procedures  may  be  systematically  improved. 
Thermally  stimulated  current  (TSC)  measurements  on  irradiated  oxides  are  used 
here  to  indicate  the  energy  level  position  of  the  traps  or,  at  a  minimum,  the 
thermal  activation  energies  of  the  trapped  holes.  These  measurements  are  per¬ 
formed  in  conjunction  with  conventional  capacitance- voltage  (C-V)  techniques  to 
track  the  amount  of  charge  residing  in  the  oxide.  The  measurements  performed 
under  this  program  are  intended  to  be  complemented  with  other  basic  mechanisms 
information  being  obtained  using  X-ray  Photoelectron  Spectroscopy  (XPS)  on  the 
same  materials  at  the  Jet  Propulsion  Laboratory.  The  XPS  technique  yields 
detailed  information  concerning  the  physical  chemistry  of  the  Si -Si 0^  inter¬ 
face,  the  distribution  of  network  strain,  and  the  effects  of  the  introduction 
and  removal  of  charge  at  the  interface. 

The  full  set  of  experimental  data  (TSC,  C-V,  and  XPS)  is  ultimately  to 
be  acquired  for  a  matrix  of  specimens  which  cover  the  conventionally  used  gate 
oxide  processes.  The  matrix  includes  aluminum  gate  and  polysilicon  gate  tech¬ 
nologies,  dry  and  wet  (pyrogenic)  oxide  growth  processes,  and  both  radiation- 
hard  and  radiation-soft  varieties.  In  the  latter  regard,  the  soft  specimens 
are  prepared  by  the  usual  formula  for  fabrication  of  hard  oxides,  with  the 
addition  of  a  high-temperature  nitrogen  anneal.  Sample  fabrication  is  being 
performed  at  the  Sandia  National  Laboratories.  The  specimens  are  intended  to 
be  representative  of  the  current  state  of  the  art  in  radiation-hard  technolo¬ 
gies.  They  also  assume  particular  significance  because  of  the  increased  reli¬ 
ance  expected  to  be  placed  on  the  new  Center  for  Radiation-Hardened  Microelec¬ 
tronics  (CRM)  at  Sandia  by  the  Department  of  Defense,  the  Department  of  Energy, 
NASA,  and  their  contractors.  This  study  is  therefore  appropriate  and  timely. 
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During  the  initial  phase  of  the  investigation  reported  on  here,  the 
experimental  technique  was  established  and  demonstrated  using  MOS  capacitors 
previously  fabricated  at  NRTC.  Additionally,  a  mask  set  for  the  fabrication 
of  the  full  matrix  of  samples  was  furnished  by  NRTC  to  Sandia,  and  fabrication 
of  the  matrix  of  specimens  was  underway  at  the  end  of  the  current  reporting 
period.  An  initial  shipment  of  samples  was  received  prior  to  the  end  of  the 
reporting  period,  but  not  in  sufficient  time  for  tests  to  be  completed.  Results 
of  those  measurements  will  be  described  in  the  final  report  for  Contract 
DNA001-82-C-0055. 

3.2  BACKGROUND 

Prior  to  discussing  the  specific  experimental  technique  and  the  initial 
results,  we  present  relevant  background  material.  A  number  of  different  mea¬ 
surement  techniques  have  been  applied  in  the  past  to  the  task  of  identifying 
the  nature  of  hole  traps  in  SiO,,  films.  These  include  annealing  studies, 
thermoluminescence,  photodepopulation,  and  prior  measurements  of  thermally 
stimulated  current  and  capacitance.  Various  experiments  have  also  been  con¬ 
ducted  with  a  view  toward  correlating  the  introduction  rate  of  trapped  holes 
with  oxide  growth  parameters  and  such  characteristics  as  mobile  ion  concen¬ 
trations. 

A  variety  of  annealing  results  hfrve  been  obtained.  Woods  and  Williams^ 

report  the  disappearance  of  trapped  positive  charge  after  anneals  above  150°C. 
14 

Aitken  found  that  80%  of  the  radiation-induced  charge  is  removed  upon  anneals 

to  400UC.  Regarding  the  activation  energy  for  the  release  of  trapped  holes, 

again  a  variety  of  values  is  reported.  Nicollian  quotes  an  activation  energy 

of  0.35  eV.  lindmayer^  found  activation  energies  to  range  from  1.3-1. 5  eV. 

Using  thermally  stimulated  capacitance  measurements,  Manchanda  et  al.^  found 

18 

an  activation  energy  for  anneal  of  1-1.5  eV.  Bakowski  performed  thermally 
stimulated  current  measurements  on  dry  and  wet  oxides  and  found  two  values  of 
activation  energy,  one  in  the  range  of  0.8-1. 3  eV,  and  the  other  1.8-2. 4  eV. 
(The  method  employed  by  Bakowski  to  extract  activation  energies  differed  from 
that  employed  here,  and  if  the  latter  method  were  applied  the  calculated  acti¬ 
vation  energies  would  be  considerably  smaller  (by  as  much  as  40%). 


« 
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Hughes  and  Thomas  also  performed  extensive  TSC  measurements  in  a  study 
which  compared  hard  and  soft,  dry  and  wet  oxides  using  several  diagnostic  tech¬ 
niques.  The  results  were  somewhat  ambiguous  due  to  problems  in  sample  fabrica¬ 
tion.  The  wet  specimens  appeared  to  be  dominated  by  mobile  ion  contributions  to 
the  thermally  stimulated  current.  Dry  oxides,  which  were  not  thus  affected, 
did  not  exhibit  the  structure  in  the  TSC  response  observed  by  Bakowski.  Sample 
behavior  was  also  found  to  be  altered  by  pre-irradiation  bias-temperature 
stressing. 

In  contrast  to  the  above  results,  optical  measurements  of  the  photoioni¬ 
zation  cross  section  of  trapped  holes  have  consistently  yielded  no  photodepopu- 

20 

lation  of  the  trapped  holes.  Harari,  Wang,  and  Royce  were  unable  to  photo¬ 
depopulate  the  traps  in  measurements  performed  at  77°K  with  photons  up  to  5  eV 

21 

in  energy.  Snow,  Grove,  and  Fitzgerald  found  the  threshold  of  emptying  to  be 
at  4.3  eV,  which  is  the  energy  at  which  election  excitation  from  the  silicon 
valence  band  to  the  SiO^  conduction  band  becomes  energetically  allowed  (internal 
photoemission).  It  is  quite  likely  that  it  was  this  process  which  they  observed, 
rather  than  photodepopulation. 

Bringing  these  various  results  into  consonance  is  one  of  the  tasks  being 
addressed  in  the  present  study.  Three  possible  explanations  of  the  difference 
between  optical  and  thermal  values  for  the  apparent  hole  binding  energies  come 
to  mind.  First,  within  the  framework  of  the  smal 1-polaron  hopping  model  of  hole 
transport,  it  is  possible  that  the  reference  level  for  the  thermal  activation 
energies  is  not  the  valence  band  edge,  but  rather  that  energy  at  which  the  hole, 
or  small  polaron,  L  comes  mobile  within  the  Si 0^  bandgap.  The  polaron  binding 
er-rgy  is  on  the  order  of  0.5  eV.  Second,  it  is  possible  that  what  is  being  ob- 
s*_’  /ed  in  the  TSC  measurements  is  an  electrochemical  reaction  which  renders  the 
interface-trapped  holes  mobile.  The  third  is  the  conventional  interpretation  of 
discrepancies  between  optical  and  thermal  energies  in  terms  of  the  Franck-Condon 
principle.  For  the  present  case,  this  principle  implies  that  the  trapped  hole, 
and  its  associated  lattice  distortion,  remains  rigid  during  ar  optical  inter¬ 
action  in  that  the  system  does  not  have  time  to  adjust  on  the  time  scale  of  the 
photoexcitation.  If  we  describe  the  mechanical  freedom  of  the  complex  (hole  + 
l'ttv-e)  in  terms  of  a  configurational  coordinate,  the  dependence  of  system 
energy  on  this  coordinate  has  harmonic  oscillator  properties  as  shown  in 
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Figure  i2  for  the  ground  state  and  a  nearby  configuration  which  involves  a 
different  lattice  adjustment.  The  optical  transition  is  shown  in  the  figure, 
•?.nd  is  seen  to  be  considerably  larger  in  energy  than  the  energy  separation  of 
the  two  inequivalent  sites  which  would  govern  the  thermal  excitation  process. 
The  continuation  of  the  present  work  should  shed  light  on  these  hypotheses. 


Figure  12.  Configuration  diagram  showing  the  dependence 
of  energy  on  the  configurational  coordinate  q  for  two 
inequivalent  sites  of  the  trapped  hole.  A  large  differ¬ 
ence  between  the  optical  and  the  thermal  excitation  ener¬ 
gies  is  apparent  for  transitions  from  one  site  to  another. 


3.3  THEORY 

The  technique  of  thermally  stimulating  the  depopulation  of  traps  has  been 
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used  for  many  years  in  the  study  of  solids.  These  include  thermally  stimu¬ 
lated  luminescence  (thermoluminescence  or  "glow  curves")  in  which  electron 

23 

traps  are  depopulated  in  insulators,  and  more  recently,  the  use  of  thermally 

stimulated  ionic  currents  (TSIC)  in  studies  of  ionic  impurities  in  Si09  and 

24  c 

Al^O^.  Thermoluminescence  is  also  a  standard  technique  for  radiation  dosime¬ 
try.  The  extension  of  the  thermally  stimulated  current  technique  to  the  study 
of  hole  traps  in  SiO^  films  is  a  straightforward  application  of  previous  work. 
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As  in  all  TSC  measurements,  the  fundamental  problem  is  the  interpretation  of 
the  observed  current  signature. 

A  TSC  measurement  is  made  by  applying  a  bias  across  the  insulator  in  an 
MOS  capacitor  and  recording  the  current  as  the  sample  temperature  is  raised 
(Figure  13(a)).  The  measurement  is  repeated  in  order  to  determine  the  back¬ 
ground  current  which  is  unrelated  to  thermal  depopulation  of  traps,  such  as 
that  due  to  carrier  injection.  The  difference  between  the  two  currents  is 
interpreted  as  arising  from  the  thermal  depopulation  of  charged  traps.  The 
form  of  the  temperature  ramp  (linear,  exponential,  hyperbolic)  impacts  the 
conversion  of  the  TSC  spectra  into  a  form  from  which  trap  energies  and  the 
density  of  states  may  be  extracted.  The  fact  that  both  electrons  and  holes 
may  be  present,  that  discrete  and  continuous  trap  levels  may  exist,  and  that 
recombination  and  retrapping  processes  can  occur  complicates  the  analysis.  If 
retrapping  can  be  ignored,  as  is  usually  assumed  in  the  case  of  high-field  TSC, 
many  of  these  problems  can  be  eliminated.  In  the  present  case  of  irradiated 
SiO^,  it  is  believed  that  the  hole  traps  have  much  larger  capture  cross  sec¬ 
tions  than  electron  traps  and  so  the  TSC  spectra  should  reflect  the  hole  traps 
predominantly. 

The  simplified  rate  equation  governing  the  depopulation  of  traps  is 

dn  -E/kT  ( t) 

ST  '  '"tve  ’  (1) 

where  n^  is  the  occupied  trap  density,  v  is  the  "attempt- to-escape"  frequency, 

E  is  the  activation  energy  of  the  detrapping  process,  and  T  is  the  temperature, 

which  is  a  function  of  time.  In  general,  an  analytic  solution  of  this  equation 
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does  not  exist.  However,  as  shown  by  Hickmott  for  TSIC  studies  with  discrete 

energy  levels,  a  hyperbolic  temperature  ramp  allows  an  analytic  integration  of 

the  rate  equation.  Using  high-field  TSC  for  a  continuous  distribution  of  energy 
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levels,  Simmons  and  coworkers  have  shown  that  the  TSC  spectrum  is  a  direct 
image  of  the  energy  distribution  of  occupied  traps.  For  a  linear  temperature 
ramp  with  a  constant  rate  3,  the  sample  temperature  is  given  by 

T  =  to  +  ei  , 
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Figure  13.  Schematic  representation  of  a  TSC 
measurement.  J  is  the  measured  current,  while 
N(E)  is  the  derived  density  of  occupied  states. 

where  Tq  is  the  initial  temperature  and  t  is  time.  The  measured  current  is 
plotted  versus  temperature;  a  representative  curve  is  shown  in  Figure  13(b) 
(curve  1).  A  second  measurement  is  made  to  establish  the  baseline  (Figure 
13(b),  curve  2).  The  TSC  spectrum  is  the  difference  between  the  two  currents 
(Figure  13(c)).  For  a  distribution  of  occupied  hole  traps,  N(E),  where  N(E) 

3 

has  dimensions  of  traps/cm  -eV,  the  current  per  unit  area  I  is  given  by 


I  =  \  qLDN(E)  ,  (3) 

where  q  is  the  magnitude  of  unit  charge,  L  the  thickness  of  the  insulating 
layer,  and 

D  =  1.2  ME  -  Ev)/T  .  (4) 
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The  energy  scale  (eV)  is  linearly  related  to  temperature  (°K)  by 

E-E  =  T  (1.92  x  10"4  log  £  +  3.2  x  10'4)  -  0.0155  .  (5) 

v  p 

Because  the  dependence  of  (E-  E  }  on  v  is  weak,  the  specific  value  of  v  used 

v  10  12  -1 

is  not  critical.  A  value  between  10  and  10  sec  is  usually  taken.  Explicit 
determination  of  the  best  value  is  possible  by  measuring  the  spectra  at  two  dif¬ 
ferent  ramp  rates,  and  f2 ,  and  measuring  the  temperatures,  Tj  and  at 
which  a  prominent  peak  occurs.  The  emission  rate,  v,  is  then  calculated  from 
v  =  10y  where 


To  log  P2  -  T.  log  pj 

y  . - — - 1.66 

i2  ij 


3.4  EXPERIMENTAL  METHOD 


(6) 


A  significant  amount  of  effort  was  expended  during  this  program  to  estab¬ 
lish  the  experimental  method.  Difficulties  were  encountered  in  meeting  the 
requirements  of  sensitivity  to  sub-picoamp-level  currents  in  the  presence  of  a 
noisy  environment  arising  from  the  heater  and  thermal  ramp  control  circuits. 
Difficulties  were  also  encountered  in  eliminating  spurious  currents  from  the 
sample  headers.  The  various  experimental  configurations  are  described  in  this 
section. 

Initially,  the  use  of  a  scanning  electron  microscope  (SEM)  was  considered 
for  irradiation  and  measurement.  The  SEM  produces  large  ionizing  dose  rates, 
which  allow  the  desired  doses  of  104  to  10^  rads(Si)  to  be  accumulated  in  min¬ 
utes.  The  Co-60  source  at  our  disposal,  by  contrast,  is  limited  in  practice 
to  dose  rates  of  about  20  rads/sec.  Further,  the  SEM  is  equipped  with  a  hot 

_5 

stage,  and  offers  a  dark,  vacuum  environment  at  the  10  torr  level.  Some  ques¬ 
tions  can  be  raised  about  the  equivalence  of  low-energy  electron  irradiation 
and  Co-60  irradiation.  It  was  our  intention  to  address  these  concerns  experi¬ 
mentally  by  performing  at  least  some  experiments  at  the  Co-60  source  for  pur¬ 
poses  of  comparison.  The  principal  question  in  this  regard  concerns  the  effec¬ 
tiveness  in  generating  holes  which  has  been  found  to  differ  for  the  two  sources 

of  irradiation.  This  question  has  been  most  recently  addressed  by  Brown  and 
26 

Dozier.  For  our  purposes,  however,  absolute  calibration  is  of  secondary 
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importance.  Rather,  it  is  essential  that  reproducibility  of  the  total  dose 
from  run  to  run  be  assured.  In  this  regard,  it  should  be  noted  that  our  SEM 
has  recently  been  equipped  with  a  unique  capability  of  measuring  the  beam  cur¬ 
rent  impinging  on  a  sample  without  perturbation  or  interception,  thus  making 
precise  beam  current  measurements  during  irradiation  possible. ^  Additional 

considerations  in  the  case  of  low-energy  electron-beam  irradiation  are  the  non¬ 
uniformity  of  dose  with  depth  into  the  sample,  and  the  possibility  of  lateral 

nonuniformities  in  dose  due  to  the  Gaussian  intensity  profile  of  finite  width. 
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These  issues  have  been  addressed  quantitatively  by  Galloway  and  Roitman.  It 
was  decided  to  irradiate  with  the  SEM  during  the  feasibility  study,  to  be  fol¬ 
lowed  by  Co-60  irradiations  in  subsequent  work. 

When  the  SEM  was  employed  for  TSC  measurements,  it  was  found  that  the 
noise  coupling  from  the  heater  circuit  to  the  current  measurement  circuit  was 
excessive.  Hence,  the  SEM  was  used  only  for  irradiation  purposes,  and  TSC 
measurements  made  in  a  separate  fixture.  Accordingly,  a  fixture  was  fabri¬ 
cated  in  which  specimens  mounted  in  TO-5  headers  could  be  heated  under  vacuum 
in  a  manner  which  minimized  the  coupling  between  the  two  circuits.  A  schema¬ 
tic  of  the  apparatus  is  shown  in  Figure  14.  Semi-rigid  coax  was  used  for  all 
electrical  feedthroughs.  The  sample  was  mounted  within  a  copper  block  of 
small  thermal  mass  in  order  to  minimize  the  thermal  time  constant  and  ensure 
a  uniform  temperature  environment  for  the  specimen. 

TSC  measurements  made  with  this  fixture  yielded  excessive  currents  at 
high  temperatures.  These  currents  were  unrelated  to  trapping  of  charges  within 
the  oxides,  and  were  traced  to  the  headers.  The  currents  were  found  to  be  a 
function  of  temperature,  to  be  unrelated  to  the  biases  applied  to  the  header 
pins,  and  to  prevail  even  under  conditions  of  uniform  temperature,  thereby 
tending  to  rule  out  the  effects  of  thermal  gradients.  We  were,  in  effect, 
doing  TSC  measurements  on  the  headers  themselves.  A  number  of  different  device 
packages  were  investigated,  and  none  were  found  which  yielded  satisfactory 
levels  of  leakage  currents  at  elevated  temperatures.  This  included  packages 
employing  a  sapphire  dielectric.  Typical  results  of  TSC  measurements  per¬ 
formed  on  a  TO-33  header  with  no  MOS  sample  are  shown  in  Figure  15. 

Since  the  current  levels  observed  with  headers  would  obscure  the  detrap¬ 
ping  of  holes  in  the  SiC^.  a  guarded  configuration  using  a  sapphire  dielectric 
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Figure  14.  Schematic  diagram  of  fixture  for  performance  of 
thermally  stimulated  current  measurements  on  specimens 
mounted  on  TO-5  headers. 


was  tried.  This  is  illustrated  in  Figure  16(a).  The  TSC  currents  arising 
within  the  sapphire  substrate  were  again  found  to  be  prohibitively  large.  Mani¬ 
festly,  all  current  return  paths  to  the  hot  substrate  needed  to  be  eliminated. 
The  first  valid  TSC  measurements  were  made  in  the  above  fixture  using  a  flying 
lead,  as  illustrated  in  Figure  16(b).  A  2-mil  gold  bonding  wire  was  connected 
to  both  the  sample  and  the  substrate  for  stability.  After  mounting  of  the 
sample,  and  soldering  to  the  output  lead,  the  connection  to  the  substrate  was 
cut.  Measurements  were  made  in  this  way  while  another  fixture  was  fabricated 
which  provided  for  probe  contacts  to  the  sample.  Consideration  was  given  to 
the  use  of  an  apparatus  employed  for  readout  of  TLDs,  since  such  equipment 
provides  well-controlled  selectable,  linear  thermal  ramp  rates.  This  approach 
was  abandoned  due  to  the  fact  that  these  instruments  commonly  employ  noisy  SCR 
control  circuits. 
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Figure  15.  Example  of  thermally  stimulated  currents 
observed  for  a  TO-33  header.  An  18-mV  thermocouple 
reading  corresponds  to  330°C. 

The  circuits  employed  for  linear  sample  heating  are  shown  in  Figure  17. 
The  derivative  of  the  temperature  is  obtained  by  use  of  the  high-pass  filter 
on  the  PAR  113  preamp.  The  amplified  signal  is  summed  with  the  manually  ad¬ 
justable  reference  voltage  at  the  summing  junction  of  a  Kepco  power  operational 
amplifier  which  supplies  power  to  the  heater.  The  measurement  thermocouple  was 
mounted  on  an  alumina  substrate  for  isolation.  The  reference  voltage  was  ad¬ 
justed  as  necessary  to  maintain  constant  temperature  derivative,  as  monitored 
by  means  of  an  X-Y  plotter. 

The  experimental  apparatus  using  probe  contacts  to  the  sample  is  shown  in 
Figure  18.  This  apparatus  was  used  for  both  TSC  and  C-V  measurements.  Changes 
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Figure  16.  Sample  mount  using  metallized  alumina  substrate  in  a 
guarded  configuration  (a),  and  using  a  "flying  lead"  connection 
technique  (b)  which  eliminates  return  to  heated  substrate. 


with  respect  to  the  previous  fixture  consisted  mainly  in  the  use  of  two  inde¬ 
pendent  heater  circuits,  one  of  which  was  connected  to  a  dc  supply,  and  the 
other  of  which  was  connected  to  the  servo  circuit.  Additionally,  independent 
thermocouples  were  provided  for  measurement  and  control.  To  assure  thermal 
integrity  the  measurement  thermocouple  was  thermally  lagged  to  the  sample  block, 
and  for  maximum  thermal  response  frequency  the  control  thermocouple  was  imbed¬ 
ded  in  the  sample  block.  Samples  were  die-attached  to  gold-chrome  metallized 


POWER  SUPPLY 


TC  INPUT 


£ 


(0-< 

1 

< 

< 

< 

50V)* 

) 

>62  kn 
> 

,  IN 

KFPCO 

2.4  k ft 

3PERATI0NAI 

AMPLIFIER 

TSC  CHAMBER 


OUT 


OUT 


20  ••C/Inch 


KEITHLEY  616 


’MANUALLY  AD  JUSTED  TO  MAINTAIN  CONSTANT  TEMPERATURE  TIME  DERIVATIVE 


Figure  17.  Diagram  of  circuits  used  for  linear  sample  heating 
for  data  acquisition. 
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Figure  18.  Experimental  TSC 
apparatus  incorporating  probe 
contacts  to  MOS  test  samples. 


25-mil  alumina  substrates,  which  were  in  turn  clamped  to  the  copper  surface 
of  the  sample  block  by  means  of  lugs  fabricated  from  machineable  ceramic.  The 
measurement  thermocouple  was  clamped  to  the  substrate  using  one  of  these  lugs, 
and  the  metallization  was  scribed  using  a  diamond  saw  in  such  a  way  as  to 
maintain  electrical  isolation  of  the  thermocouple  from  the  current  measure¬ 
ment  circuit.  The  bias  supply  wire  was  then  clamped  under  the  other  lug. 

The  associated  electronics  were  modified  as  well,  as  illustrated  in 
Figure  19.  Thermocouple  signals  were  conditioned  using  isolation  amplifiers 
with  integral  reference  compensation.  The  temperature  ramp  procedure  was  the 
following:  1)  establish  the  appropriate  setpoint  for  the  closed-loop  control 
system,  and  2)  increment  the  bias  heater  supply  as  necessary  in  order  to  main¬ 
tain  reasonably  stable  output  on  the  Kepco  operational  amplifier.  A  TSC  run  re¬ 
quires  about  15  minutes  under  the  chosen  conditions  of  ramp  rate  (0.3°C/sec). 


SETPOINT 


Figure  19.  Diagram  of  the  final  TSC  system. 
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An  Alessi  Model  CW-100  probe  tip  was  initially  used  to  contact  to  the 
bonding  pad  of  the  capacitor.  The  probe  tip  consisted  of  a  5-mil  tungsten  "S" 
spring  etched  to  a  point  of  0.07-p.m  radius  which  is  welded  to  a  Dumet  wire. 

The  wire  was  mounted  in  an  alumina  tube  which  was  mounted  in  a  joystick  type 
probe  manipulator.  A  fine  wire,  welded  to  the  Dumet  wire,  was  connected  to 
rigid  coax  which  was  brought  out  of  the  vacuum  chamber  by  a  BNC  vacuum  con¬ 
nector.  Even  with  the  aid  of  a  microscope  and  a  C-V  plotter  as  a  sensor,  it 
was  very  difficult  to  contact  the  bonding  pad  with  the  "S"  probe  without  caus¬ 
ing  damage  to  both  probe  tip  and  pad.  The  "S"  probe  was  replaced  with  an 
Everett/Charles  spring  contact  with  a  rounded  tip.  This  w:s  found  to  produce 
a  good  electrical  contact  without  damaging  either  probe  or  sample. 

The  current  measurement  was  made  using  a  Keithley  616  digital  electrome¬ 
ter  operating  in  the  fast  (feedback)  mode.  The  capacitor  substrate  was  biased 
to  +3.0V  during  TSC  measurements.  The  output  of  the  electrometer  was  fed  to 
the  Y  input  of  one  of  the  X-Y  recorders.  The  C-V  and  G-V  measurements  were 
made  using  a  PAR  Model  410  C-V  plotter  operating  at  1  MHz  with  15  mV  modula¬ 
tion  signal.  A  typical  voltage  ramp  rate  used  was  0.2  V/sec. 

As  previously  mentioned,  initial  plans  called  for  irradiation  of  the  MOS 
capacitors  using  a  Co-60  source  for  the  feasbility  tests.  However,  because  of 
logistics  (the  source  was  located  off  the  premises)  the  feasibility  tests  were 
conducted  by  irradiating  with  the  SEM.  As  previously  described,  the  disadvan¬ 
tages  of  using  the  SEM  were:  1)  complexity  in  estimating  radiation  dose  because 
of  the  appreciable  effect  of  the  metallization,  and  2)  lateral  nonuniformity  of 
dose  due  to  finite  beam  spread.  The  SEM  was  operated  at  a  5-kV  accelerating  po¬ 
tential  and  currents  of  5  and  20  picoamps.  These  currents  had  uncertainties  of 
at  least  25°'.  In  order  to  estimate  the  energy  deposited  in  the  Si 0^ >  an  inverse 
energy  dependence  of  -dE/dx  was  assumed.  Integrating  the  energy  loss  equation 
incorporating  losses  in  the  A1  gate  electrode,  and  correcting  for  the  nonuni- 
formity  of  the  scan  gave  an  approximate  total  dose  of  0.75  Mrad ( S i 0^ )  received 
by  the  sample  used  in  the  feasibility  test.  Since  the  irradiation  required 
that  the  substrate  be  biased  at  a  negative  potential,  a  small  fixture  was  built 
for  mounting  the  sample  and  electrically  grounding  the  gate  electrode  without 
appreciably  masking  the  capacitor. 
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The  procedure  for  irradiation  in  the  SEM  was  to  position  the  sample, 
take  a  scan  at  a  minimum  power  level,  adjust  the  position  and  fully  irradiate 
using  multiple  scans.  C-V  measurements  were  made  pre  and  post  irradiation. 

If  the  flatband  shift  was  insufficient,  the  procedure  was  repeated. 

3.5  EXPERIMENTAL  RESULTS 

The  experimental  approach  was  governed  by  a  general  concern  that  the  use 
of  probes  and  thermocompression  bonded  wires  directly  on  the  thin  gate  oxides 
might  give  rise  to  spurious  effects.  Since  the  initially  envisioned  experi¬ 
mental  method  involved  the  use  of  wire  bonds,  a  structure  consisting  of  both 
field  and  gate  oxides  was  called  for.  To  maintain  compatibility  with  samples 
furnished  to  JPL  for  XPS  measurements,  the  field  oxide  was  to  be  grown  on  all 
samples,  but  would  be  completely  removed  from  the  XPS  specimens  prior  to  growth 
of  the  gate  oxide.  Masks  for  the  fabrication  of  1  mm  x  1  mm  field  plates,  with 
a  4  mil  x  4  mil  bonding  pad  on  the  field  oxide,  were  furnished  by  js  to  Sandia 
Laboratories.  The  matrix  of  specimens  was  prepared  at  Sandia  over  the  course 
of  the  program.  While  waiting  for  these  samples,  the  measurement  method  was 
tested  by  means  of  MOS  capacitors  previously  fabricated  at  Northrop.  Results 
of  these  measurements  are  presented  below. 

The  Northrop  MOS  capacitor  design  is  shown  schematically  in  Figure  20. 

A  silicon  dioxide  film  of  ~1000  A  thickness  was  grown  in  dry  oxygen  at  1100°C 
on  polished  (111)  surfaces  of  3  Q-cm  n-type  silicon  wafers.  An  aluminum  field 
plate  1  mm  in  diameter  and  600  A  thick  was  evaporated  through  a  mask  over  the 
oxide.  An  n+  layer  about  1.5  nm  thick  was  diffused  into  the  back  side  of  the 
wafer  and  a  gold  layer  deposited  to  make  good  electrical  contact. 

The  sequence  of  measurements  for  the  feasibility  sample  was  as  follows. 
Before  and  after  each  TSC  run  a  C-V  measurement  was  made.  For  the  TSC  data 
run,  the  temperature  was  ramped  with  a  negative  gate  bias  applied.  After 
cooling,  a  C-V  measurement  was  made,  then  the  TSC  baseline  was  run.  Another 
C-V  measurement  was  made  immediately  before  irradiation.  During  room- tempera¬ 
ture  irradiation  a  positive  bias  was  applied  to  the  gate.  After  irradiation 
a  C-V  measurement  was  made.  Then  the  TSC  data  was  taken,  another  C-V  measure¬ 
ment  made,  the  TSC  baseline  measured,  and  finally  C-V  was  measured  again.  Usu¬ 
ally  there  was  no  difference  between  C-V  measurements  before  and  after  the  TSC 


48 


'  "  i 


25  jam 


At  (0.06  pm) 


Figure  20.  Schematic  diagram  showing  construction 
of  the  MOS  capacitor  used  as  a  test  vehicle. 

baseline  after  irradiation.  This  indicated  that  no  hole  traps  were  left  occu¬ 
pied  in  the  energy  range  accessible  to  thermal  detrapping  for  the  temperatures 
used  in  these  measurements. 

Figure  21  shows  the  results  of  C-V  measurements  prior  to  irradiation. 
Note  that  there  was  a  negative  flatband  shift  each  time  a  TSC  measurement  was 
made.  The  TSC  measurement  is  equivalent  to  a  negative-bias  stressing  of  the 
MOS  oxide.  The  phenomenon  responsible  for  negative  flatband  shifts  observed 

under  negative-bias  stressing  is  called  slow  trapping  or  the  negative-bias 

28 

instability.  The  negative  flatband  shift  implies  that  there  is  an  increase 
in  the  effective  positive  charge  in  the  oxide.  Holes  and  ions  cannot  be  re¬ 
sponsible  for  this  increase  since  they  move  to  the  gate  under  negative  bias 
and  would  not  contribute  to  the  flatband  shift.  Since  TSC  cannot  distinguish 
between  positive  charge  flowing  to  the  gate  or  negative  charge  flowing  to  the 
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Figure  21.  Pre- irradiation  C-V  measurements: 
(a)  pre-TSC,  (b)  post-TSC  data,  (c)  post-TSC 
baseline. 


Si-SiO^  interface,  a  plausible  explanation  for  the  negative-bias  instability 
is  that  it  is  caused  by  electron  emission  from  pre-existing  neutral  traps  in 
the  oxide.  If  these  traps  are  located  at  a  mean  distance  x  from  the  interface, 
where  L  is  the  oxide  thickness,  then  the  flatband  shift  A  VCD  and  integrated 

r  D 

current  Qc  as  determined  from  the  TSC  measurement  would  be  given  by 
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where  Qq  is  the  charge  left  in  the  oxide  after  electron  emission  and  CQx  is 
the  oxide  capacitance.  Solving  for  x  gives 
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In  Figure  22  the  TSC  spectra  are  plotted.  These  curves  are  the  differ¬ 
ence  between  the  TSC  data  and  the  TSC  baseline  measurements.  The  method  of  con¬ 
verting  the  temperature  scale  to  an  energy  scale  was  described  in  Section  3.3. 
The  energy  distribution  of  these  states  covers  a  broad  band  from  0.9  to  1.45  eV 
without  any  distinct  peaks.  Integrating  the  TSC  data  gives  a  Qc  of  512  pC. 
Using  a  AVpg  of  0.25  V  places  the  position  of  the  electron  traps  approxi¬ 
mately  150  I  from  the  gate  in  a  1000  A  thick  oxide. 


042  E  -  Ev  (eV) 

Figure  22.  Results  of  TSC  measurements  on  the  Northrop  M0S 
capacitor.  Uncertainty  in  the  current  measurements  is  indicated 
by  the  error  bar. 
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It  was  observed  that  after  two  days  at  room  temperature  without  bias, 
the  flatband  voltage  returned  to  its  pre-TSC  value.  This  observation  was 
somewhat  surprising  since  the  energy  level  distribution  of  the  trapped  holes 
ranged  from  0.9  to  1.45  eV,  with  an  average  at  1.15  eV.  The  Qq  of  ~600  pC 
corresponds  to  an  average  hole  emission  rate  (p)  of  1.8  x  10^  holes/sec.  If 
this  were  a  thermally  activated  process  it  should  be  governed  by 

-E/kT 

p  =  ve  ,  (10) 

where  E  is  the  trap  activation  energy  and  v  is  the  "attempt  frequency".  Sub- 

23  -1 

stituting  1.25  eV  for  E  and  solving  for  v  gives  v  =  3  x  10  sec  .  After 
comparing  this  value  to  the  usually  accepted  range  of  10^-10^  sec-'*'  it  was 
concluded  that  this  room  temperature  instability  was  not  consistent  with  a 
thermally  activated  process. 

The  post-irradiation  C-V  results  are  shown  in  Figure  23.  For  comparison 
the  initial  pre-irradiation  C-V  measurement  is  also  displayed  (curve  a ).  Curve  b 
was  taken  immediately  after  the  SEM  electron  beam  irradiation,  while  curve  c 
was  recorded  after  overnight  room- temperature  annealing.  Again  a  room  tempera¬ 
ture  instability  is  observed  for  hole  traps;  however,  these  traps  are  located 
near  the  interface  as  opposed  to  those  associated  with  the  negative-bias  insta¬ 
bility,  which  are  located  approximately  150  A  from  the  gate.  After  the  TSC  mea¬ 
surements,  curve  d  was  obtained.  There  were  no  differences  between  C-V  measure¬ 
ments  obtained  after  the  TSC  data  and  TSC  baseline  runs.  Curve  d  demonstrates 
that  most  of  the  charge  introduced  by  the  irradiation  was  released  during  the 
TSC  measurement.  The  small  residual  negative  flatband  shift  may  be  caused  by 
holes  in  deeper  traps  or  by  the  negative-bias  instability. 

Observe  in  Figure  23  that  the  post-irradiation  curves  show  not  only  a 
shallower  slope  as  compared  to  the  pre-irradiation  curve,  but  also  three 
shoulders.  It  is  reasonable  to  interpret  this  behavior  as  being  due  to  par¬ 
allel  capacitance  caused  by  the' lateral  nonuniform  irradiation  induced  by  the 
raster  scan  of  the  SEM  or  by  shadowing  caused  by  the  fixture  grounding  the 
gate  electrode.  Ignoring  the  very  small  shoulder  located  at  approximately 
-14  V,  it  is  possible  to  estimate  the  fraction  of  area  of  the  M0S  capacitor 
associated  with  each  parallel  capacitor.  This  is  done  by  assuming  that  of  the 
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Figure  23.  Post-irradiation  C-V  measurements:  (a)  pre¬ 
irradiation  (shown  for  reference),  (b)  post-irradiation, 

(c)  post-irradiation  after  overnight  room- temperature 
annealing,  (d)  post-irradiation,  post-TSC. 

three  regions  in  the  C-V  plot,  one  has  both  capacitors  in  accumulation,  another 
both  in  depletion,  while  the  center  region  has  one  in  each.  Using  this  ap¬ 
proach,  it  is  estimated  that  19«  of  the  MOS  capacitor  was  irradiated  to  give 
a  flatband  voltage  of  -3.25  V  while  the  rest  of  the  capacitor  had  a  flatband 
voltage  of  -10  V.  Calculating  the  stored  charge  induced  by  the  irradiation  by 
weighting  the  contribution  from  each  capacitor  gives  1750  pC. 

The  post-irradiation  TSC  spectra  is  displayed  in  Figure  22.  The  trap 
distribution  ranges  in  energy  from  0.7  to  1.3  eV.  It  is  possible  that  there 
are  two  peaks,  one  at  0.9  eV  and  the  other  at  1.1  eV.  Because  of  the  uncer¬ 
tainty  in  the  current  measurement  in  this  feasibility  study,  this  assignment 
must  be  considered  tentative.  The  total  charge  obtained  by  integrating  the 
TSC  spectra  is  1680  pC.  This  is  in  excellent  agreement  with  the  equivalent 
total  charge  obtained  from  the  flatband  shifts  (1750  pC). 
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3.6  CONCLUSIONS 


The  primary  goal  of  this  year's  work  was  to  substantiate  that  the  TSC 
method,  coupled  with  C-V  measurements,  is  feasible  for  studying  the  occupancy 
of  hole  traps  at  the  Si/Si^  interface.  This  has  been  accomplished.  Using  a 
Northrop- fabricated  capacitor  as  a  test  vehicle,  a  hole  trap  distribution 
ranging  in  energy  from  0.7  to  1.3  eV  was  found,  with  energy  peaks  at  0.9  and 
1.1  eV  being  tentatively  assigned.  The  charge  contributing  to  the  observed 
TSC  accounted  for  most  of  the  observed  radiation-induced  flatband  shift. 

The  Northrop  sample  also  displayed  a  slow  trapping  or  negative-bias  in¬ 
stability.  The  electron  traps  responsible  for  this  phenomenon  were  located 
approximately  150  A  from  the  gate  of  the  1000  A  oxide  with  an  energy  distribu¬ 
tion  ranging  from  0.9  to  1.5  eV.  The  effects  produced  by  these  traps  were  ob¬ 
served  to  disappear  at  room  temperature  within  two  days.  This  instability  was 
inconsistent  with  a  thermally  activated  process. 

With  the  feasibility  of  the  method  proven,  the  program  will  continue  in 
1982  with  the  investigation  of  hole  trapping  in  pedigreed  M0S  oxides  fabri¬ 
cated  by  Sandia.  These  measurements  should  permit  correlation  of  process  vari¬ 
ables  and  hardness  with  hole  trap  properties.  It  is  expected  that  comparisons 
will  be  made  for  the  following  combinations  of  oxide  properties:  hard  versus 
soft,  A1  vs  poly-Si  gate,  and  wet  vs  dry  oxide  processing.  Of  major  concern 
is  whether  there  are  other  contributions  to  the  TSC  spectra,  such  as  ion  con¬ 
tamination  or  electron  tunneling.  Our  measurements  will  be  correlated  with 
concurrent  chemical  analysis  of  the  interface  undertaken  on  similar  samples 
by  F.  Grunthaner  at  JPL. 
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SECTION  4.0 


IONIZING  RADIATION  EFFECTS  ON  HgCdTe  MIS  DEVICES 

A  study  of  the  effects  of  ionizing  radiation  on  HgCdTe  MIS  devices  was 
performed  during  this  program.  The  results  obtained  were  presented  at  two 
conferences  during  1981:  the  IRIS  Infrared  Detector  Specialty  Group  Meeting 
and  the  IEEE  Conference  on  Nuclear  and  Space  Radiation  Effects.  Those  papers 
were  published  and  are  reproduced  here  in  Appendices  B  and  C. 
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Abstract 

Results  of  analytical  and  experimental  studies  of 
the  mechanisms  of  neutron  effects  on  VLSI  structures  are 
described  in  this  paper.  The  issue  of  whether  a  single 
neutron-produced  defect  cluster  can  produce  a  signifi¬ 
cant  amount  of  permanent  damage  ("hard  error")  in  a  VLSI 
cell  is  addressed.  Properties  of  average  defect  clusters 
in  silicon  are  calculated  as  a  function  of  neutron  en¬ 
ergy,  including  dimensions,  defect  densities,  and  local 
electrical  properties.  The  expected  effects  of  a  single 
cluster  on  MOS  and  bipolar  VLSI  geometries  are  described, 
and  predictions  are  made  of  the  neutron  fluence  required 
to  produce  a  cluster  at  a  critical  location  in  an  inte- 
:  iiuuil  lR  wHzVEdfl  taa'tl  I  # 

neutron-irradiated  small -geometry  bipolar  test  chips  are 
described  and  compared  with  predictions  based  on  cluster 
calculations.  Much  of  the  present  data  is  accounted  for 
satisfactorily 'in  terms  of  an  average  cluster  model.  The 
present  results  suggest  that  hard  errors  will  occur  in 
VLSI  structures  at  relatively  modest  neutron  fluerces. 
This  situation  may  cause  lower  bounds  to  be  placed  on 
device  sizes  for  applications  requiring  radiation  toler¬ 
ance. 

Introduction 

Single  particle  effects,  such  as  alpha-particle-in¬ 
duced  soft  errors  in  MOS  memories, 1  are  of  considerable 
practical  importance  at  present.  Transient  upset  is  pro¬ 
duced  through  generation  of  electron-hole  pairs  by  the 
incident  particle,  and  this  effect  becomes  increasingly 
important  as  integrated  circuit  geometries  are  reduced. 
The  issue  of  whether  a  single  particle  can  produce  a 
significant  amount  of  permanent  damage  ("hard  error")  in 
a  VLSI  cell  has  also  arisen.  In  particular,  we  have  con¬ 
sidered  the  case  of  defect  clusters,  or  disordered  re¬ 
gions,  produced  in  silicon  by  fast  neutrons.  For  levels 
of  integration  through  LSI,  circuit  dimensions  typically 
are  much  larger  than  the  largest  dimension  of  a  defect 
cluster.  However,  this  situation  does  not  necessarily 
hold  for  VLSI.  Thus,  it  is  important  that  the  impact  of 
disordered  regions  on  smal 1 -geometry  circuits  be  evalu¬ 
ated  through  analysis  and  experiment. 

In  this  paper,  results  of  analytical  and  experimen¬ 
tal  studies  of  the  mechanisms  of  neutron  effects  on  VLSI 
structures  are  described.  We  first  present  calculated 
properties  of  average  disordered  regions  in  silicon  as  a 
function  of  neutron  energy.  Next,  the  expected  effects  of 
a  single  defect  cluster  on  MOS  and  bipolar  VLSI  geome¬ 
tries  are  described,  followed  by  a  prediction  of  the 
neutron  fluence  required  to  produce  a  cluster  at  a  crit¬ 
ical  location  in  an  integrated  circuit.  Experimental  re¬ 
sults  on  a  small -geometry  bipolar  test  chip  are  then 
presented  and  compared  with  predictions. 

Disordered  Region  Properties 

In  this  section  we  present  results  of  calculations  of 
the  properties  of  neutron-produced  disordered  regions  in 
silicon  as  a  function  of  neutron  energy.  Included  are  di¬ 
mensions  of  the  damaged  zone  and  the  surrounding  space- 
charge  region,  local  defect  densities,  the  effects  of 
short-term  annealing,  and  determination  of  local  elec¬ 
trical  properties  within  a  disordered  region. 

C 1 uster  S izes,  Defect  De nsities,  and  Space-Charge - 
Region  t  imers  ions' 

? 

Sigmund  developed  an  expression  which  gives  an 
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upper  limit  for  the  number  (N)  of  atoms  displaced  by  an 
energetic  recoil: 

N  ■  •  (!) 
n 

In  this  equation,  Ercj  is  that  portion  of  the  primary  re¬ 
coil  energy  Er  which  goes  into  displacement  processes 
and  Ea  is  an  effective  displacement  energy,  which  is 
~25ev  for  most  materials.^  That  value  of  Ej  is  employed 
here  for  silicon,  and  is  about  twice  the  minimum  recoil 
energy  necessary  before  displacement  is  observed. 4>6 
Equation  (1)  then  reduces  to  N  =  17Er(j  ( Ercj  in  keV). 

Stein6  calculated  average  values  for  Er  and  Erd  over  the 
fwhj*.  1  |-!I  lUHtn  a 

into  account.  His  values  were  used  here  to  determine  the 
number  of  displacements  as  a  function  of  neutron  energy. 

The  total  range  R  of  primary  recoil  atoms  in  silicon 
was  used  in  obtaining  the  dimensions  of  an  individual 
neutron-produced  defect  cluster.  To  determine  R  as  a 
function  of  incident  neutron  energy,  the  model  of  Lind- 
hard  et  al.7,8  was  employed.  Sigmund  and  coworkers  (Ref.  9 
and  citations  therein)  performed  calculations  of  mean 
damage  distributions  produced  by  ion  implantation,  and 
their  work  is  used  here  to  estimate  the  dimensions  of 
defect  clusters.  Based  on  their  calculations,  a  prolate 
spheroidal  shape  is  assumed  for  the  damage  distribution. 
For  energetic  ions  incident  on  a  semiconductor  surface, 
the  ions  will  be  distributed  slightly  deeper  into  the 
material  than  the  concomitant  damage  produced  during  bom¬ 
bardment.8^  Sigmund  et  al.9  determined  the  ratio  of  the 
average  ion  range  distribution  (x)r  to  the  averaqe  dam¬ 
age  distribution  <x)[>.  To  obtain  the  major  axis  (a)  of 
the  prolate  spheroid,  we  divided  primary  recoil  range 
values  by  appropriate  values  for  <.x>r/ <x>(j .  To  obtain 
values  for  th'  minor  axis  (b)  as  a  function  of  neutron 
energy,  we  used  the  tApression  b  =  BUy2^)1  >  where  y 

is  the  transverse  dimension  measured  from  the  x  axis, 
and  \y^ >o  's  the  n,ean  square  value  of  that  dimension. 
Values  for  )D  were  obtained  from  transverse  strag¬ 
gling  calculations  of  Sigmund  et  al.9  The  volume  of  a 
prolate  spheroid  is  given  by  n  ab’/6.  Volumes  are  given 
in  Figure  1,  along  with  the  energy  dependence  of  the  de¬ 
fect  density  per  disordered  region.  Defect  densities  were 
determined  by  dividing  the  number  of  displacements  N  by 
the  volume. 

As  examples  of  the  results  of  the  present  calcula¬ 
tions,  figure  2  shows  average  defect  clusters  produced 
by  incident  neutrons  of  three  different  energies.  As  the 
neutron  energy  increases,  the  clusters  become  more  elon¬ 
gated  and  depart  significantly  from  a  spherical  shape. 

In  the  original  model  for  neutron-produced  disordered  re¬ 
gions,  proposed  by  Gossick,'61  a  spherical  shape  was  as¬ 
sumed  for  simplicity.  (Curtis,  11, l2  Gregory, '■*  and 
Holmes'4  have  qiven  further  consideration  to  the  details 
of  that  model.)  The  damaged  region  is  represented  in  the 
Gossick  model  by  a  sphere  of  radius  Rj.  The  interior  of 
the  damaged  region,  represented  by  a  sphere  of  radius  R0, 
is  taken  to  be  compensated  intrinsic.  The  damaged  region 
between  R0  and  Rj  is  charged  by  majority  carriers  trapped 
at  defect  centers.  Surrounding  the  damaged  region  is  a 
depleted  space-charge  region  of  radius  R^. (Disordered  re¬ 
gions  present  a  potential  barrier  to  majority  carriers 
and  a  potential  well  to  minority  carriers.)  In  the  present 
work,  the  equations  developed  by  Gossick  for  spherical 
symmetry  were  employed  and  the  results  then  adapted  to  a 
prolate  spheroidal  shape.  Figure  3(a)  presents  results 
obtained  for  donor  concentrations  (Nj)  of  10'6  and  ID'' 
cm*3  at  three  neutron  energies. 
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Figure  1.  Oisordered  region  volume  and  defect  density 
per  disordered  region  vs  neutron  energy  for  neutron- 
irradiated  silicon. 


Figure  2.  Illustration  of  damaged  regions  produced  in 
silicon  by  impinging  neutrons  of  three  energies. 
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The  above  calculations  of  disordered  region  proper¬ 
ties  apply  to  the  damage  that  exists  immediately  after 
neutron  irradiation.  This  damage  is  unstable  in  that  sub¬ 
stantial  annealing  occurs  at  room  temperature  on  a  time 
scale  of  seconds  following  pulsed  bombardment. 15  We  mod¬ 
ified  our  results  to  take  short-term  annealing  into  ac¬ 
count.  Based  on  previous  experiments ,15.16  we  use  10  as 
an  approximate  factor  by  which  N,  the  initial  number  of 
defects,  is  reduced  during  the  annealing  period.  We  fur¬ 
ther  assume  that  the  defect  cluster  volume  (i.e.,  dam¬ 
aged  portion)  does  not  change  during  annealing  so  that 
the  defect  concentration  Nt  is  reduced  by  a  factor  of  10. 
R?  and  R0  were  then  recalculated  for  this  situation, 
thereby  yielding  altered  space- charge- region  dimensions 
following  annealing.  Figure  3  compares  disordered  region 
dimensions  for  two  cases,  assuming  a  doping  concentra¬ 
tion  of  1015  cm*3;  1)  before  sho't-term  annealing; 

2)  after  short-term  annealing.  Reduction  in  the  dimen¬ 
sions  of  the  surrounding  space-charge  regions  is  evi¬ 
dent  after  annealing. 

Local  Electrical  Properties 

We  now  estimate  local  electrical  properties  within 
defect  clusters  produced  in  silicon  by  1-MeV  and  14-MeV 
neutron  bombardment.  The  probability  of  neutron  scatter¬ 
ing  per  centimeter  of  material,  denoted  by  P' ,  is  given 
by  Not_for  material  thicknesses  small  compared  to  1  cm, 
where  N  is  the  atomic  density  and  ct  is  total  neutron 
cross  section.  (Values  of  oj  have  been  compiled  as  a 
function  of  neutron  energy  by  Garber  and  Kinsey. 1')  The 
number  of  disordered  regions  produced  per  unit  volume  at 
a  specific  fluence  ;  is  then  given  by  P' Denoting  dis¬ 
ordered  region  volume  (damaged  portion)  as  V^,.,  the  frac¬ 
tional  volume  taken  up  by  disordered  regions  is  then 
given  by  P'iVtfr-  The  approach  used  here  to  determine 
local  electrical  properties  is  to  employ  known  macro¬ 
scopic  damage  coefficients  and  convert  them  to  local  co¬ 
efficients  by  using  this  fractional  volume  expression. 

Carrier  generation  rate  is  considered  first.  Neutron 
damage  in  depleted  regions  of  silicon  devices  causes  an 
in>  r  a  e  in  1^  y  rw  rsiilm  R  ?  1  Ht-avbj  P*- paiXJ k  iJl 

an  Increase  in  dark  current.  The  macroscopic  neutron-in¬ 
duced  increase  in  thermal  generation  rate  can  be  ex¬ 
pressed  as18 

AG  =  n.p/2Kg  ,  (2) 


where  n^  is  the  intrinsic  carrier  concentration  and  Kg 
is  the  generation-lifetime  damage  coefficient.  For  fis¬ 
sion  neutron  irradiation  (—  IMeV),  Kq  was  determined^8 
to  be  7.0  x  108  n-sec/cm?  after  short-term  annealing,  so 
that  AG  =  1  x  10l 5  pairs/cm’-sec  at  room  temperature  fora 
fluence  of  101?  n/cm?.  Using  the  same  experimental  ap¬ 
proach  as  in  Ref.  18,  a  value  for  Kg  of  3.1  x  108  n-sec/ 
cm'  was  determined  for  14-MeV  neutron  bombardment  in  the 
present  investigation.  This  yields  aG  =  2.3  x  10^5  pairs/ 
crrw-sec  as  the  macroscopic  neutron- induced  generation 
rate.  The  local  generation  rate  is  given  by 

iG  ( local  )  =  AG/P';Vdr  =  ni/2KgP'Vdr  .  (3) 

19  3 

Equation  (3)  yields  AG(local )  ^  7x  10  pairs/cm  -sec 
for  1-MeV  neutrons.  For  14-MeV  neutrons,  we  obtain  a 
local  generation  rate  of  -2.5  x  id8  pairs/cm8-sec, 
which  is  a  factor  of  ~30  less  than  that  for  the  1-MeV 
case.  This  difference  occurs  because  the  calculated  vol¬ 
ume  per  disordered  region  is  a  stronger  function  of  neu¬ 
tron  energy  than  the  number  of  defects  per  disordered 
region,  and  hence  the  defect  density  decreases  with  in¬ 
creasing  energy  (Figure  l). 

In  the  above  calculations,  the  volume  used  was  that 
for  the  damaged  portion  of  a  disordered  region  and  did 
not  include  the  surrounding  space-charge-region  volume. 

This  is  reasonable  because  carrier  generation  is  of  im¬ 
portance  when  a  defect  cluster  is  in  a  device  depletion  j 

region,  and  in  that  zone  the  space-charge  region  sur¬ 
rounding  the  cluster  loses  its  significance. 

We  now  estimate  the  minority-carrier  recombination  I 

rate  within  a  single  disordered  region  in  neutral  mate¬ 
rial.  The  recombination  rate  per  carrier  (R')  is  given  by 
the  reciprocal  of  the  minority-carrier  lifetime  (t).  The 
macroscopic  neutron- induced  increase  in  R'  is  given  by  j 

AR'  =  i(  1/t)  =  p/Kr  .  «) 

where  Kj.  is  the  recombination  lifetime  damage  coeffi- 
i'i  F  fiui:  ■  ox .  ti  .n  bombardments  K,.  for  n-iyrt-  ma¬ 
terial  at  low  injection  levels  is  ~  108  n-sec/cm'.iy 
Using  Eq.  (4),  we  obtain  a  value  for  aR*  of  ~  10'  sec"1 
at  a  fluence  of  10 18  n/cm^.  The  local  recombination 
rate  can  be  expressed  as  I 

AR'(local)  =  AR'/P';Vdr  '  l/P'VdrKr  .  (5) 
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_  io15  cm"3  Table  1  summarizes  the  macroscopic  and 

local  electrical  properties  determined  for 
disordered  regions.  (The  estimated  macroscopic 
removal  rate  for  14-MeV  neutrons  is  partially 
based  on  data  of  Cleland  et  al .  )  These  de¬ 

terminations  are  meaningful  only  in  the  con¬ 
text  of  the  average  cluster  model  assumed 
here.  In  fact,  the  damage  is  not  distributed 
uniformly  on  a  microscopic  scale,  but  is  typ¬ 
ically  described  in  terms  of  narrow  damage 
tracks  and  terminal  clusters. “ 


Table  1.  Local  electrical  properties  in  a 
disordered  region.  (Macroscopic  values  of 
generation  rate  and  recombination  rate  are 
for  a  fluence  of  lO'2  n/c.n2.  Recombination 
rates  are  for  low-resistivity  n-type  material 
at  low  injection  levels.) 
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1  MeV 

Microscopic  loot 

11  MeV 

Microscopic  tool 

1  x  1015  ~7xl019 

2.3xl015  -2.5x10** 
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Recombination  Rate 

Per  Carrier  (sec’S 
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- 1  x  107  -7x10** 
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1-0  “ 
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Figure  3.  (a)  Illustration  of  the  spice-charge  region  surround¬ 
ing  neutron-produced  defect  clusters  in  silicon,  (b)  Illustra¬ 
tion  of  the  reduction  in  space-charqe-region  dimensions  follow¬ 


ing  short-term  annealing. 

For  1-MeV  neutron  bombardment,  we  obtain  aR' (local)-  7 
x  10'1  sec*!.  14-MeV  neutrons  are  a  factor  of  2  to  3 
more  effective  than  fission  neutrons  in  degrading  the  re¬ 
combination  lifetime.20  Using  a  factor  of  2.5,  this  re¬ 
duces  K  to  ~4  x  104  n-sec/cm2.  Equation  (5)  then  yields 
iR* (local)  -  2.5  x  10l°  sec’*  for  14-MeV  neutrons.  (In 
these  calculations,  the  volume  used  was  that  of  the  dam¬ 
aged  portion  of  a  disordered  region  since  that  is  where 
recombination  actually  occurs.) 

We  now  consider  neutron-induced  carrier  removal.  The 
rate  at  which  the  majori ty-carrier  concentration  (n)  de¬ 
creases  in  neutron- irradiated  silicon  is  referred  to  as 
the  carrier  removal  rate.  At  a  given  fluence,  the  removed 
carrier  concentration  is  given  by  (An/AA);.  A  disordered 
region  can  be  viewed  as  a  void  containing  no  free  major¬ 
ity  carriers,  and  this  void  has  a  volume  V*.  (where  V*jr 
is  the  totaj  volume,  including  the  space-charge  regipn/ 
since  the  space-charge  region  also  contains  a  negligible 
concentration  of  free  majority  carriers.  Thus,  by  assump¬ 
tion,  the  local  carrier  removal  rate  in  a  disordered  re¬ 
gion  is  infinite,  which  leads  to  a  method  for  checking 
the  consistency  of  the  present  calculations  with  experi¬ 
mental  values  for  removal  rate. 


Figure  4.  Fraction  of  carrier  concentration  remaining  vs 
neutron  fluence.  Results  based  on  experimental  findings 
of  Stein  and  Gereth  as  well  as  calculated  results  are 
shown  for  three  dopant  concentrations. 


Stein  and  Gereth  measured  carrier  removal  rates  as  a 
function  of  dopant  concentration.2!  For  fission-neutron- 
irradiated  n-type  silicon,  they  obtained  ~4,  ~6,  and -9 
neutron*!  cm*  1  for  Nj  values  of  10l4,  1015,  4n<j  iol&,  re¬ 
spectively.  The  fraction  f*  of  the  initial  carrier  con¬ 
centration  n  remaining  after  irradiation  can  be  ex¬ 
pressed  as  ji  -  (:.n/A,  )»•  j/n.  This  expression  was  used 
with  Stein  and  Gereth 's  removal  rates  to  obtain  the 
solid  curves  shown  in  Figure  4.  The  fractional  volume  oc¬ 
cupied  by  disordered  regions  is  given  by  and  thus 

one  minus  this  quantity  is  the  fractional  remaining  car¬ 
rier  concentration.  Calculated  results  for  1-MeV  neutron 
bortiardment  are  shown  dashed  in  Figure  4.  Considering 
the  assumptions  made,  the  agreement  with  experimental  re¬ 
sults  is  reasonable. 


Potential  Effects  of  a  Single  Disordered  Region  on 
Small -Geometry  MQS  and  STpolar  Circuits 

Figure  3  indicates  that  sizes  of  average  disordered 
regions  are  comparable  to  anticipated  dimensions  of  VLSI 
circuits,  particularly  in  the  case  of  14-MeV  neutron 
bombardment.  This  suggests  that  a  single  disordered  re¬ 
gion  could  produce  significant  effects  in  such  a  circuit. 
We  consider  here  the  potential  effects  of  a  single  de¬ 
fect  cluster  on  the  electrical  properties  of  example  MOS 
and  bipolar  cells  in  a  VLSI  array.  The  importance  of 
these  effects  will  depend  on  the  size  of  a  VLSI  cell 
compared  to  the  dimensions  of  a  disordered  region. 
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Figure  5  illustrates  a  single  disordered  region  in 
an  MOS  transistor  placed  at  two  potentially  critical 
locations:  a)  linking  the  source  and  drain  electrodes; 
b)  penetrating  the  qate  oxide.  In  case  a),  the  space- 
charge  region  surrounding  thp  defect  cluster  will  con¬ 
nect  the  junction  space-charge  regions  associated  with 
the  source  and  drain,  giving  rise  to  punch-through.  This 
condition  will  resu’t  in  increased  drain  leakage  currert. 
A  second  expected  effect  is  reduction  of  the  channel  mo¬ 
bility  due  to  introduction  of  new  scattering  centers  in 
the  silicon  near  the  Si-SiC^  interface.  A  third  effect  is 
a  shift  in  threshold  voltage,  and  this  can  occur  via 
three  mechanisms  First,  if  defects  associated  with  the 
disordered  region  intersect  the  Si -  Si  0^  interface,  the 
Fermi  level  at  the  interface  will  be  pinned,  resulting 
in  a  larger  gate  voltage  being  required  to  turn  the  de¬ 
vice  on.  Second,  if  the  interface  is  intersected  only  by 
the  space-charge  region  surrounding  the  cluster,  there 
should  also  be  a  threshold  shift  since  the  surface  will 
now  be  easier  to  invert.  This  shift  will  be  in  a  direc¬ 
tion  opposite  to  that  associated  with  Fermi-level  pin¬ 
ning.  Third,  both  electron  and  hole  traps  are  likely  to 
be  introduced  in  the  oxide  if  the  defect  cluster  pene¬ 
trates  into  it  (such  as  in  Fig.  5(b)).  These  traps  could 
give  rise  to  threshold  voltage  shifts  if  the  device  were 
subsequently  bombarded  with  ionizing  radiation.  Introduc¬ 
tion  of  new  traps  all  the  way  through  the  oxide  is  also 
expected  to  give  rise  to  increased  gate  leakage  current 
due  to  enhanced  trap-assisted  tunneling  at  the  boundaries 
of  the  oxide.  Associated  with  this  leakage  current  in¬ 
crease  is  an  expected  reduction  in  oxide  breakdown  volt¬ 
age. 

p 

Figure  6  presents  an  example  bipolar  device  (II 
cell)  containing  a  disordered  region  in  a  critical  loca¬ 
tion.  If,  as  depicted,  the  disordered  region  penetrates 
the  narrow  base  region  (~0. 1-0.2  urn  in  state-of-the-art 
lzl  circuits),  two  effects  are  expected.  Leakage  current 
should  increase  due  to  punch- through  and  current  gain 
should  degrade  due  to  the  enhanced  local  recombination 
rate. 

Prediction  of  Neutron  Fluences  Required  to  Produce 
D’tsordere3~fteq'ibns  at  Critical  locations 
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Figure  5.  Single  disordered  regions  in  MOS  transistors: 
a)  linking  source  and  drain;  b)  penetrating  the  gate 
oxide. 
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The  neutron  fluence  required  to  produce  a  single  de¬ 
fect  cluster  at  a  critical  location  in  a  smal 1 -geometry 
Integrated  circuit  is  now  calculated.  For  simplicity,  we 
treat  the  case  where  the  orientation  of  the  cluster  with 
respect  to  the  circuit  geometry  is  unimportant.  An  exam¬ 
ple  of  this  situation  is  shown  in  Figure  6.  The  disord¬ 
ered  region  is  depicted  as  completely  penetrating  the 
base  region  of  the  vertical  npn  transistor.  If,  for  ex¬ 
ample,  a  14-MeV  neutron-produced  cluster  is  created  in 
the  center  of  the  base  region,  emitter  and  collector  will 
be  linked  by  the  cluster  regardless  of  its  orientation. 


Consider  a  silicon  chip  of  area  A  and  thickness  t. 
There  will  be  certain  regions  on  the  chip  that  can  be 
considered  as  critical  locations  in  terms  of  the  ex¬ 
pected  effect  of  a  single  disordered  region.  For  the  ex¬ 
ample  of  Figure  6,  the  collector  area  Is  the  critical 
area  and  the  intrinsic  basewidth  is  the  critical  thick¬ 
ness.  We  denote  the  critical  area  for  one  device  on  the 
chip  as  Ac,  the  critical  thickness  as  tc,  and  the  number 
of  devices  per  chip  with  critical  regions  as  D.  The  frac¬ 
tional  critical  area,  f>,  can  then  be  expressed  as  ACD/A, 
and  the  fractional  critical  thickness,  ft,  as  t,/t.  If 
one  cluster  is  produced  in  a  single  chip,  the  cluster 
density  is  diver,  by  1/At  clusters/cm3.  The  fluence  ;j 
required  .ircduce  one  cluster  per  chip  is  then  given 
by  1/AtP'.  The  critical  fluence  required  to  produce 
one  cluster  per  chip  at  a  critical  location  can  then  be 
expressed  as 
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Figure  6.  Single  disordered  region  penetrating  the  In¬ 
trinsic  base  and  linking  emitter  and  collector  in  a  bi¬ 
polar  transistor  (I^l  cell). 

When  cluster  dimensions  exceed  the  basewidth  In  a 
bipolar  transistor,  then  the  appropriate  critical  thick¬ 
ness  to  use  with  Eq.  (6)  is  larger  than  the  actual  base- 
width  since  the  cluster  may  be  centered  outside  the  base 
region  but  still  penetrate  It  completely.  As  an  example, 
consider  14-MeV  neutron  bombardment  (P*  *  0.09)  of  a 
chip  containing  10^  devices  with  critical  dimensions  of 
lxl  x  0.1  pm.  Assuming  the  cluster  orientation  of  Fig¬ 
ure  6,  the  effective  value  of  tc  Is  ~0.35  pm.  Equation 
(6)  then  yields  a  critical  fluence  of  ~3  x  10$  n/cm‘. 

Experimental  Studies 

A  bipolar  test  chip  was  designed  and  fabricated  at 
our  laboratories  and  electrical  measurements  were  then 
performed  on  test  devices  before  and  after  14-MeV  neu¬ 
tron  bombardment.  The  results  of  those  initial  measure¬ 
ments  are  described  in  this  section. 

Experimental  Details 

2 

A  six-mask,  "washed-emitter"  l  l  process  was  used  to 
fabricate  the  bipolar  test  chip  which  contained  devices 
with  ~2-pm  minimum  geometry.  Included  were  bipolar 
transistors  with  the  same  structure  as  the  vertical  npn 
device  in  Figure  6.  Collector  dimensions  were  measured 
for  45  devices  and  the  result  was  2.3  ±  0.3  pm.  Thus;  the 
actual  collector  geometry  Is  typically  2.3  x  2.3  pm,  and 
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the  basewidth  (W)  is  ~0.2  pm.  The  test  chip  also  con¬ 
tained  pinch  resistors  with  approximate  dimensions  of 
2.3  x  3.0  x  0.2  pm.  The  intrinsic  base  doping  concentra¬ 
tion  was  ~5  x  10^6  cnr3,  which  was  also  the  concentra¬ 
tion  for  the  p-type  pinch  resistors. 

For  transistors,  measurements  of  collector  current 
(Ic)  were  made  for  three  values  of  forward  bias  applied 
to  the  base-emitter  junction:  0.66,  0.69,  and  0.72  V.  By 
measuring  Ic,  we  obtain  information  solely  concerning  re¬ 
combination  in  the  narrow  base  region  (the  intrinsic 
base),  as  opposed  to  recombination  in  the  wider  base  re¬ 
gion  adjoining  it  (the  extrinsic  base).  Recombination  in 
the  critical  intrinsic  base  region  was  of  primary  con¬ 
cern  here,  as  depicted  in  Figure  6.*  Transistor  leakage 
current  was  also  measured  for  a  collector  bias  of  1.6  V 
with  the  base  and  emitter  grounded.  For  pinch  resistors, 
resistance  was  determined  by  measuring  the  current  flow 
foi  three  values  ut  ffppttw  bias,  fl,  10u,  and  150  m7. 
These  voltages  were  selected  to  be  sufficiently  low  to 
avoid  pn  junction  current  flow  or  resistor  pinch-off. 
Device  temperatures  were  maintained  at  31°C  during  mea¬ 
surements,  with  temperature  control  to  within  0.1®C. 
Device  voltages  were  maintained  within  0.1  mV  of  desired 
values. 

To  improve  precision,  five  measurements  were  made  on 
each  device  both  before  and  after  irradiation.  The  mean, 
standard  deviation,  and  percent  standard  deviation  (i.e., 
expressed  as  a  percentage  of  the  mean  value)  were  then 
calculated  for  every  set  of  readings.  Unirradiated  pinch 
ftirjtjrt  acj  tee* Msfc,r»  w<rrt  enipluyed  ftr  re-k^wnet 
purposes  to  verify  the  integrity  of  our  experimental 
technique.  As  discussed  below,  we  expected  neutron-in¬ 
duced  changes  in  device  properties  to  be  relatively 
small.  Thus,  it  was  important  that  the  percent  standard 
deviation  for  a  given  set  of  five  measurements  be  small 
compared  to  device  changes.  We  were  able  to  achieve  such 
conditions.  For  measurements  of  Ic  in  transistors,  at 
0.72  V  applied  the  percent  standard  deviation  averaged 
over  the  132  devices  that  were  irradiated  was  0.33%  be¬ 
fore  Irradiation  and  0.30%  after  bombardment.  The  cor¬ 
responding  quantities  for  six  reference  devices  were 
0.33%  and  0.25%,  respectively.  For  current  measurements 
in  pinch  resistors  at  150  mV.  the  percent  standard  devi¬ 
ation  averaged  over  90  irradiated  devices  was  0.27%  be¬ 
fore  irradiation  and  0,28%  afterwards.  For  10  reference 
devices,  the  average  percent  standard  deviation  was 
0.24%  before  irradiation  and  0.25%  when  measured  after 
irradiation  of  the  other  90  devices. 

Irradiations  were  performed  at  room  temperature  at 
the  Lawrence  Livermore  Laboratory  RTNS-II  14-MeV  neutron 
source.  Post-irradiation  measurements  were  performed 
typically  within  one  to  two  weeks  after  bombardment. 

Expected  Changes  in  Device  Properties 

We  now  consider  the  expected  effects  of  a  single  14- 
MeV  neutron-produced  defect  duster  on  the  electrical 
properties  of  the  present  pinch  resistors  and  bipolar 
transistors.  The  damaged  portion  of  an  average  duster 
is  calculated  to  have  a  major  axis  of  0.45  i*m  and  a 
minor  axis  of  0.21  pm.  For  a  doping  concentration  of 
5  x  10*6,  ^e  surrounding  space-charge-region  thickness 
was  determined  to  be  negligible.  We  then  assumed  that 
the  cluster  is  a  non-conducting  void  (an  overestimate) 
and  calculated  the  change  in  resistance  that  occurs  when 
one  cluster  is  added.  For  the  cluster  orientation  of 
Figure  6.  which  is  the  most  likely  situation  to  occur  in 
our  experiments,  the  calculated  increase  in  resistance 
is  ~0.5%. 

IT  n  ,n_r 

Current  gain  measurements  were  not  made  here  because 
recombination  in  regions  other  than  the  intrinsic 
tutu  1  met  u  ut  turiwr  «*rta*4e. 


For  bipolar  transistors,  we  again  assume  the  cluster 
orientation  shown  in  Figure  6.  After  short-term  anneal¬ 
ing,  the  defect  concentration  calculated  here  for  a  14- 
MeV  neutron-induced  cluster  (1.8  x  lO1”  cm*2)  is  less 
than  the  dopant  concentration.  For  this  situation,  the 
damaged  zone  is  fully  charged,  the  surrounding  space- 
charge-region  width  is  negligible,  and  the  potential 
barrier  is  reduced  to  near  zero.  The  charged  defects  in 
the  cluster  will  still  be  effective  for  recombination, 
but  not  as  effective  as  when  a  potential  well  for  minor¬ 
ity  carriers  is  presented  by  the  disordered  region.  To 
calculate  the  amount  of  recombination,  we  first  estimate 
the  lifetime  within  the  damaged  region.  Using  the  ap¬ 
proach  given  above  for  determining  local  recombination 
rate,  a  lifetime  of  0.14  nsec  was  determined  for  an  as¬ 
sumed  damage  coefficient  (appropriate  for  p-type  mate¬ 
rial  at  relatively  high  injection  levels)  of  1.5  x  105 
n-sec/cm2.  To  determine  the  corresponding  diffusion 
length  (L),  a  value  for  dlffusivity  (0)  in  the  cluster 
is  needed.  Since  electron  mobilities  in  polycrystalline 
silicon  are  as  high  as  100,  as  an  estimate  we  employ 
100  cm2/V-sec  as  the  local  mobility  in  a  defect  cluster. 
This  assumption  leads  to  D  »  2.6  cm2/sec  and  L  s=0.2  pm. 
Thus,  the  calculated  diffusion  length  is  comparable  to 
the  basewidth  for  the  present  devices. 

A  carrier  collection  volume  assumed  equal  to  the  vol¬ 
ume  within  one-half  the  basewidth  of  the  cluster  boundary 
was  employed.  Appropriate  device  equations  were  used  to 
determine  the  fraction  of  carriers  entering  the  cluster, 
either  from  the  emitter  or  from  the  surrounding  collec¬ 
tion  volume,  that  reounbUie.  omd  on  this  caleutelfor. 
and  volume  considerations,  the  expected  reduction  in 
collector  current  for  the  present  geometry  is  ~0.8%. 
Associated  with  a  small  change  in  Ic  is  a  relatively 
large  change  in  current  gain.  For  the  intrinsic  base  re¬ 
gion  of  the  present  devices,  if  the  pre-irradiation  gain 
is  100, the  post-irradiation  gain  is  calculated  to  be 
~50.  A  larger  effect  is  expected  in  smaller  devices. 

In  general,  for  W  fixed,  as  the  collector  area  is  re¬ 
duced  the  effects  of  a  single  cluster  are  enhanced. 

These  effects  can  be  reduced  by  reducing  the  basewidth. 
The  important  point  is  that  a  single  duster  produced 
by  14-MeV  neutrons  is  predicted  to  have  a  significant 
effect  on  current  gain  in  smal 1 -geometry  transistors. 

Experimental  Results 

Figure  7  is  a  histogram  showing  the  number  of  bi¬ 
polar  transistors  exhibiting  a  percentage  change  in  the 
mean  value  of  Ic  within  a  specific  range  (0.2%  bins) 
after  irradiation  to  3.1  x  1012  n/cm2.  The  six  unirradi¬ 
ated  reference  devices  exhibited  a  slight  decrease  in 
current,  so  the  distribution  was  shifted  slightly  to¬ 
ward  the  origin  to  account  for  this. 

The  distribution  of  percentage  change  in  resistance 
for  pinch  resistors  is  shown  as  a  histogram  in  Figure  8. 
Percent  change  in  resistance  was  observed  to  correlate 
with  the  magnitude  of  the  current,  presumably  due  to 
geometry  variations.  Using  regression  analysis,  an  ap¬ 
propriate  correction  factor  was  obtained  which  “uncor- 
related"  the  data,  thereby  assuring  comparability  of  all 
devices. 

Leakage  current  measurements  on  transistors  did  not 
reveal  any  significant  radiation-induced  increases, 
which  is  consistent  with  our  calculations  of  a  negli¬ 
gible  space-charge-region  thickness  for  a  defect  cluster 
in  the  intrinsic  base  region.  Thus,  an  appreciable 
punch-through  current  is  not  expected  for  the  present 
devices. 

Discussion 

For  the  present  2.3  x  2.3  x  0.2  iun  transistors, 
with  ten  of  these  critical  structures  per  chip,  a 
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Figure  7.  Histogram  showing  percent  change  in  mean  value 
of  collector  current  for  bipolar  transistors  at  a  base- 
to-emitter  voltage  of  0.72  V.  Analytical  results  are  also 
shown.  The  numbers  shown  correspond  to  the  number  of 
clusters  present  in  a  given  device  for  the  dashed  fit. 


evidence  that  relatively  rare  events  of  high  specific 
effectiveness  are  responsible  for  the  observed  changes. 

A  statistical  fit  to  the  data  was  made  to  determine  the 
magnitude  of  the  resistance  change  per  event  and  the 
event  probability.  A  fit  using  the  Poisson  distribution 
is  shown  in  Fig.  8.  If  unit  fluence  is  taken  to  be  that 
value  for  which  the  probability  of  obtaining  at  least 
one  cluster  in  a  critical  location  is  0.1,  then  the 
best  fit,  obtained  by  visual  estimation,  corresponds  to 
a  fluence  equal  to  16  times  the  unit  fluence.  For  this 
case,  the  probability  of  obtaining  no  clusters  in  a  par¬ 
ticular  critical  region  is  0.19;  for  one  cluster,  0.31; 
for  two  clusters,  0.26;  for  three  clusters,  0.15;  and 
1j*  4  c/at  tVatAers,  v.jt.  ttk 

was  3.1  x  10*2  n/cm?,  the  fit  to  the  data  yields  a  unit 
fluence  of  about  2  x  10“  n/cm2  as  compared  to  the  cal¬ 
culated  value  of  4  x  10^  n/cm?.  An  additional  implica¬ 
tion  of  the  analytical  fit  is  that  an  average  cluster  is 
responsible  for  a  resistance  change  of  about  0.35%,  as 
compared  to  the  calculated  value  of  -0.5%. 

From  Table  1,  an  upper  limit  for  removal  rate  is 
-30  n*l  cm-1.  For  the  present  experimental  fluence,  a 
macroscopic  resistance  change  of  -0.2%  is  expected  on 
this  basis  alone,  neglecting  mobility  changes.  For  com¬ 
parison,  the  mean  change  observed  experimentally  (Fig¬ 
ure  8)  is  0.5%.  The  fact  that  these  quantities  are  com¬ 
parable  tends  to  validate  the  analysis. 
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Figure  8.  Histogram  showing  percent  change  in  mean  value 
of  resistance  for  pinch  resistors  at  an  applied  voltage 
of  150  mV.  Analytical  results  are  also  shown.  The  num¬ 
bers  shown  correspond  to  the  number  of  clusters  present 
in  a  given  device  for  the  analytical  fit. 

critical  14-MeV  neutron  fluence  of  8.4  x  10**  n/cm^  is 
obtainea  using  Eq,  (6).  Thus,  the  actual  fluence  of  3.1 
x  1017  n/cm?  is  predicted  to  produce  a  single  cluster  at 
-4  critical  locations  per  chip.  For  pinch  resistors,  we 
obtain  a  critical  fluence  of  6.4  x  10“  n/cm?  to  produce 
one  cluster  per  chip  (containing  ten  resistors  each)  at 
a  critical  location.  Thus,  the  predicted  probability  that 
a  given  resistor  in  the  population  will  contain  one  clus¬ 
ter  is  -0.5,  Perhaps  more  realistically,  one  can  argue 
that  if  a  cluster  penetrates  only  through  a  portion,  such 
as  one-half,  of  the  critical  region  thickness,  then  sig¬ 
nificant  effects  will  be  observed.  For  this  example,  fur¬ 
ther  modification  of  tc  in  Eq.  (6)  yields  the  following 
critical  fluences:  1)  transistors:  5  x  loll  n/cm-';  2)  re¬ 
sistors:  4  x  10“  n/cm?.  These  fluences,  then,  are  our 
best  estimates  of  the  unit  critical  fluences  applicable 
to  the  present  devices. 

The  distribution  of  resistance  changes  (Fig.  8)  is 
found  to  be  sufficiently  broad  to  constitute  persuasive 

64 


As  demonstrated  by  the  analytical  fit  in  Figure  8, 
the  pinch  resistor  data  are  consistent  with  the  conven¬ 
tional  interpretation  of  neutron-produced  disordered 
regions.  On  the  other  hand,  the  bipolar  transistor  data 
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tion  of  fractional  changes  in  collector  current  is  con¬ 
siderably  tighter  than  the  distribution  for  the  resis¬ 
tors.  A  fit  to  the  data  with  the  model  which  postulates 
discrete,  damaging  events  is  shown  in  Figure  7  (dotted), 
where  the  Poisson  distribution  is  again  employed.  The 
parameters  of  the  fit  are  a  change  per  event  of  0.08% 
and  a  unit  fluence  of  1.3  x  1010  n/cm?,  or  40  times  less 
than  the  expected  unit  fluence.  Thus,  the  transistor 
data,  as  they  stand,  appear  to  be  inconsistent  with  the 
present  model  and  with  the  resistor  data. 

Since  the  experimental  distribution  for  transistors 
is  rather  narrow,  we  explored  the  issue  of  whether  the 
observed  changes  could  be  explained  cn  a  macroscopic 
basis.  An  estimate  was  made  of  the  expected  change  in 
collector  current  at  3.1  x  10l?  n/cm?  by  using  a  macro¬ 
scopic  damage  coefficient  to  determine  the  amount  of 
diffusion  length  degradation  in  the  intrinsic  base  re¬ 
gion.  The  value  obtained  was  0.051,  which  is  well  below 
the  mean  of  the  experimental  data  (abojt  21).  For  a  ma¬ 
croscopic  analysis  to  apply,  agreement  between  the  two 
values  should  be  within  the  uncertainty  of  the  mean  of 
the  experimental  data. 


The  transistor  data  do  not  appear  to  be  satisfacto¬ 
rily  described  by  the  average  cluster  model,  in  contrast 
to  the  situation  for  resistors.  This  observation  and  the 
result  of  the  macroscopic  calculation  led  us  to  hypothe¬ 
size  that  the  transistor  results  were  subject  to  a  sys¬ 
tematic  error  that  gave  rise  to  a  relatively  large  off¬ 
set  in  the  distribution.  Assuming  that  this  offset  is  a 
constant  and  does  not  contribute  to  the  statistical  dis¬ 


tribution  of  data  points,  then  the  data  can  be  fit  in  a 
manner  more  consistent  with  the  average  cluster  model 
than  the  previous  fit  (dotted  in  Fig.  7).  Using  the 


Poisson  distribution,  the  resulting  fit  for  the  case  of 
1.41  origin  offset  is  shown  dashed  in  the  figure.  The 
parameters  of  the  fit  are  a  percentage  change  in  collec¬ 
tor  current  per  cluster  of  0.2*  and  a  unit  fluence  of 
1.1  x  10“  n/cm?.  (The  calculated  values  for  these  pa¬ 
rameters  are  0.8*.  and  5  x  loll  n/rm?,  as  noted  above.) 
Experimental  confirmation  is  required,  however,  before 
we  can  confidently  describe  the  data  in  this  manner. 
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Mueller  and  Guenzer^4  performed  calculations  of 
displacement  damage  cascades  in  silicon  using  the 
MARLOWE  simulation  code.  Their  calculations  indicate 
that  the  structure  of  "defect  clusters"  varies  sub¬ 
stantially  from  that  of  the  average  cluster  calculated 
here.  They  found  that  many  "clusters"  take  the  form  of 
relatively  long  strings  of  defects  with  few  branches, 
and  thus  few  terminal  clusters.  If  their  simulations 
are  appropriate,  then  one  would  expect  the  electrical 
properties  of  small -geometry  devices  to  be  affected 
somewhat  differently.  Qualitatively,  we  expect  that 
Mueller  and  Guenzer's  defect  structures  would  produce 
significantly  less  degradation  at  a  given  fluence  and 
the  critical  fluence  would  be  lower.  The  fact  that  the 
present  resistor  data  can  be  accounted  for  satisfacto¬ 
rily  in  terms  of  average  cluster  calculations  provides 
support  for  the  average  cluster  model.  Further  support 
is  provided  by  the  electron  microscopy  work  of  8erto- 
1 o 1 1 i 2 5  ana  Sargent. 26  (in  the  latter  reference, 
craters'  observed  in  14-MeV  neutron-irradiated  sili¬ 
con  were  comparable  in  size  to  the  average  cluster  di¬ 
mensions  determined  here  (including  the  surrounding 
space-charge  region).)  Although  we  have  no  compelling 
reason  at  present  to  question  the  conventional  picture 
of  neutron- induced  disordered  regions  in  silicon,  fur¬ 
ther  work  will  be  required  to  establish  the  microscopic 
structure  of  defect  clusters. 

We  now  use  the  present  analytical  and  experimental 
results  to  predict  for  an  example  case  the  critical  flu- 
wee  ^  Ciecni  Mitfc  of  frtPrt!  efrvrt  Tit  T15T.  foAsiflE.  a 
chip  with  1C4  devices  of  1  x  1  x 0.2  ^m  critical  area  that 
is  irradiated  with  14-MeV  neutrons.  Using  Eq.  (6)  with 
an  effective  tc  of  0.25  ^m  yields  a  critical  fluence  of 
-4.5  x  109  n/cm2.  If  the  initial  current  gain  is  100, 
we  calculate  that  the  post-irradiation  gain  for  a  de¬ 
vice  containing  one  cluster  will  be  -20.  The  analyti¬ 
cal  fits  to  the  transistor  and  resistor  data  can  also 
be  used  to  estimate  critical  fluences  for  this  example. 
By  applying  appropriate  scaling  lor  yeometry  and  for 
the  number  of  devices  per  chip,  the  dashed  fit  of  Fig¬ 
ure  7  yields  a  critical  fluence  of  ~6  x  10®  n/cm?  and 
a  degraded  gain  of  -50.  Scaling  the  fit  to  the  resis¬ 
tor  data  (Figure  8)  yields  a  critical  fluence  of  ~1.5 
x  10g  n/cm2.  Thus,  for  this  example,  critical  fluences 
ranging  from  -6  x  108  to  ~4.5  x  10g  n/cm2  are  esti¬ 
mated,  depending  on  the  approach  used.  Estimated  post¬ 
irradiation  gain  values  range  from  -20  tn  -50.  De¬ 
pending  on  the  application,  such  degradation  could 
constitute  failure  of  the  device  and,  as  a  result,  the 
entire  array.  This  would  then  be  a  manifestation  of  a 
hard  error.  There  is  some  uncertainty  regarding  both 
the  fluence  at  which  this  effect  would  occur  and  the 
magnitude  of  the  9ain  change.  However,  the  important 
point  is  that  a  single  cluster  is  expected  to  have  a 
significant  effect  at  a  relatively  low  'luence. 

The  present  study  has  dealt  with  the  effects  of 
stable  neutron- induced  damage.  At  early  times  following 
a  neutron  burst,  the  effect  of  a  single  cluster  in  a 
VLSI  cell  may  be  up  to  an  order-of-magnitude  more  se¬ 
vere  than  at  long  times.  One  can  also  envision  a  tempo¬ 
rary  hard  error,  with  the  recovery  time  being  on  the 
order  of  seconds. 

Although  further  experiments  are  needed  to  re- 
k in  i+giarv  "  twCwup  frtlUfir  w3 

transistor  data,  many  of  the  present  experimental  re¬ 
sults  are  accounted  for  reasonably  well  using  average 
cluster  calculations.  Further,  our  findings  suggest 
that  hard  errors  will  occur  in  VLSI  structures  at  rela¬ 
tively  modest  neutron  fluences.  This  situation  may 
cause  lower  bounds  to  be  placed  on  device  sizes  for 
applications  requiring  radiation  tolerance. 
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ABSTRACT 

A  common  feature  of  HgCdTe  array  structures  is  the  presence  of 
insulating  layers,  either  as  gate  insulators,  passivation  layers,  or 
field  insulators.  When  most  insulators  are  exposed  to  ionizing  radi¬ 
ation,  a  charge  buildup  occurs  which  can  affect  device  operation.  This 
program  examined  ionization-induced  charge  buildup  in  HgCdTe  MIS  capac¬ 
itors  which  contained  ZnS  as  the  insulator  and  had  different  native- 
oxide  layers  between  the  HgCdTe  and  the  ZnS.  Charge  buildup  was  found 
to  occur  in  the  interfacial  region  between  the  ZnS  and  the  HgCdTe  and 
was  attributed  to  either  electron  or  hole  trapping,  depending  on  the 
operating  conditions.  The  amount  of  charge  buildup  depended  on  the 
nature  of  the  interfacial  layer,  particularly  on  the  thickness  of 
the  native-oxide  layer,  and  on  the  field  applied  during  irradiation. 

Charge  buildup  can  have  several  effects  on  array  characteristics. 

It  can  change  the  HgCdTe  surface  condition  (either  to  accumulation  or 
depletion)  which  usually  has  a  deleterious  effect  on  device  operation. 

For  example,  an  increased  effective  collection  area  for  pn  junction 
detectors,  which  could  lead  to  increased  cross-talk  in  arrays,  has 
been  observed.  Further,  charge  buildup  in  gate  insulators  will  cause 
a  shift  in  device  operating  voltage. 

INTRODUCTION 

A  common  feature  of  HgCdTe  array  structures  is  the  presence  of  insulating 

layers,  either  as  gate  insulators,  passivation  layers,  or  field  insulators.  When 

most  insulators  are  exposed  to  ionizing  radiation,  a  charge  buildup  occurs  which 

1  2 

can  affect  device  operation.  Results  of  recent  studies  '  indicate  that  the  re¬ 
sponse  of  HgCdTe  photodetectors  to  ionizing  radiation  is  dominated  by  the  insulator 
response.  This  appears  to  be  the  case  in  pn  junction  devices,  where  the  insulator 
is  the  passivation  layer  around  the  active  junction  area,  and  in  MIS  dev-'ces,  where 
the  gate  insulator  is  important.  Thus,  all  HgCdTe  array  structures  appear  to  be 
vulnerable  to  ionization-induced  charge  buildup  in  the  insulator,  and  the  damage 
process  in  these  insulators  should  be  studied  in  greater  detail. 

This  program  was  undertaken  to  examine  charge  buildup  in  insulator  structures 
commonly  used  in  HgCdTe  devices.  The  goals  of  the  program  were  to  obtain  the  infor¬ 
mation  required  to  explain  radiation-induced  degradation  of  HgCdTe  photodetector 
array  structures  and  to  develop  procedures  that  could  increase  the  hardness  of 
array  structures. 

J - 

Work  supported  by  the  Defense  Nuclear  Agency  under  Contract  DNA001-80-C-0146 . 

^F . A.  Junga  et  al.,  IEEE  Trans.  Nucl.  Sci.  2_5,  1274  (1978). 

2A.H.  Raima  and  R.A.  Cesena,  IEEE  Trans.  Nucl.  Sci.  26,  4833  (1979). 
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EXPERIMENTAL  METHODS 


Useful  test  structures  for  studying  insulator  and  interface  properties  are  MIS 
capacitors  where  the  fixed  charge  in  the  insulator  and  the  interface  state  density 
can  be  obtained  using  conventional  analysis  of  C-V  measurements.  The  devices  used 
in  this  effort  were  supplied  by  Santa  Barbara  Research  Center  and  were  fabricated 
on  HgCdTe  substrate  material  that  had  a  bandgap  of  0.2  eV  or  greater  to  minimize 
tunneling  which  can  affect  C-V  measurements.  The  standard  insulator  used  for 
HgCdTe  devices  is  ZnS.  For  the  test  devices,  it  was  prepared  by  thermal  evapora¬ 
tion  to  a  nominal  thickness  of  ~400  nm,  although  the  actual  thicknesses  varied  from 
~  280  nm  to  600  nm.  Three  different  HgCdTe  surface  treatments  were  employed  before 
deposition  of  the  ZnS.  In  one  set  of  devices,  the  HgCdTe  was  etched  and  no  inten¬ 
tional  native  oxide  was  grown  on  the  surface.  (However,  a  thin  (~  1  nm)  native- 
oxide  layer  does  form  thermally.)  In  a  second  set  of  devices,  a  70-nm  native-oxide 
layer  was  grown  anodically  using  a  KOH  solution  in  water  and  glycerine.  In  a  third 
set  of  devices,  a  6-nm  native-oxide  layer  was  grown  in  an  oxygen  plasma. J  Details 
of  the  devices  are  given  in  Table  1. 

Irradiation  exposures  were  performed  using  a  Co^°  source  with  devices  biased 
during  exposure.  The  devices  were  held  at  77°k  during  irradiation  exposure  and 
post-irradiation  measurement,  with  no  sample  heating  between  exposure  and  measure¬ 
ment.  The  effect  of  bias  alone  was  examined  by  applying  a  bias-stress  to  the 
samples  for  the  same  time  as  the  irradiation  exposure  time  (10  min) .  Charge  build¬ 
up  in  insulating  layers  was  determined  by  determining  the  flatband  voltage  shift 
from  C-V  curves. 

TABLE  1.  PROPERTIES  OF  HgCdTe  MIS  CAPACITORS 


Insulator 

ZnS  Thickness 
(nm) 

Native  Oxide 

Type 

Native  Oxide 
Thickness  (nm) 

ZnS 

590 

- 

- 

ZnS  + 

native  oxide 

500 

plasma  oxide 

6 

ZnS  + 

native  oxide 

280 

anodic  oxide 

70 

EXPERIMENTAL  RESULTS 

Irradiaticn-induced  flatband  shifts  indicative  of  charge  buildup  in  the  insula¬ 
tor  were  produced  for  all  three  of  the  insulator  structures  upon  exposure  to  ion¬ 
izing  radiation.  In  all  devices  tested,  the  charge  buildup  was  linear  with  total 

4  * 

dose,  up  to  the  maximum  test  dose  of  3.6  x  10  rad(ZnS).  Therefore,  comparisons 
between  devices  can  be  made  at  any  (lower)  dose. 

^D.R.  Rhiger  and  R.E.  Kvaas,  paper  presented  at  the  Electronic  Materials  Confer¬ 
ence,  June  24-27,  1980. 

* 

Ionizing  doses  are  stated  in  rad (ZnS)  because  the  dominant  effects  are  produced  in 
the  ZnS.  For  conversion  purposes,  1  rad(ZnS)  =  0.96  rad  (si)  for  Co®°  gamma  rays 
(E  a:  1.25  MeV)  . 
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The  dependence  of  flatband  shift  on  applied  field*  is  shown  in  Figures  1 
through  3  for  the  three  device  types.  (The  sign  convention  used  shows  the  gate-to- 
substrate  bias.)  The  bias-stress  results  (labeled  no  irradiation)  are  also  shown 
on  the  same  figures.  For  devices  with  an  intentional  native-oxide  layer  (Figures  2 
and  3),  bias-induced  flatband  shifts  were  so  small  that  no  allowance  for  them  was 
necessary.  For  devices  with  a  ZnS-only  insulator  (Figure  1),  bias-induced  flatband 
shifts  were  large  enough  that  a  correction  to  the  irradiation  data  had  to  be  made. 
The  correction  made  was  to  subtract  the  bias-induced  flatband  shifts  from  the  radi¬ 
ation-induced  flatband  shifts.  The  measured  flatband  shifts  for  these  devices  fell 
into  three  slightly  different  classes,  as  can  be  seen  in  Figure  1.  However,  the 
correction  procedure  resulted  in  the  dashed  line  shown  in  Figure  1  for  all  three 
classes  of  devices.  This  dashed  line  is  considered  to  be  the  field  dependence  of 
the  irradiation-induced  flatband  shift  for  devices  with  ZnS-only  insulators  and  is 
used  in  all  subsequent  analysis. 

To  facilitate  comparison  between  the  different  insulator  types,  irradiation- 
induced  flatband  shifts  for  the  three  insulator  structures  are  shown  in  Figure  4. 

The  sign  of  the  radiation-induced  flatband  shifts  indicates  that,  in  all  cases 
(except  at  very  low  applied  bias),  a  net  positive  charge  remains  in  the  insulator 
when  devices  are  irradiated  under  positive  bias,  while  a  net  negative  charge  re¬ 
mains  when  devices  are  irradiated  under  negative  bias.  The  amount  of  flatband 
shift  differs  for  the  different  device  types.  Because  the  ZnS  is  much  thicker  than 
the  native  oxide  in  all  device  types,  we  assume  that  all  of  the  measured  flatband 
shifts  are  produced  by  charge  that  was  initially  created  in  the  ZnS.  Assuming  that 
the  structure  of  the  ZnS  is  the  same  in  all  device  types,  the  physical  processes 
occurring  in  the  ZnS  should  be  qualitatively  the  same.  Quantitative  differences 
could  exist  because  of  the  different  ZnS  thicknesses.  However,  devices  with  the 
anodic-oxide  insulator  had  the  thinnest  ZnS  layer,  yet  showed  the  largest  flatband 
shifts,  while  devices  with  the  ZnS-only  insulator  had  the  thickest  ZnS  layer,  yet 
exhibited  the  smallest  flatband  shifts.  This  indicates  that  the  thickness  of  the 
1  .S  layer  is  not  the  cause  of  the  observed  differences,  and  suggests  that  the  inter¬ 
facial  region  (including  the  native  oxide)  between  the  ZnS  and  the  HgCdTe  is  the 
cause.  Further,  the  observation  that  differences  exist  among  the  three  device  types 
for  both  signs  of  bias  in  Figure  4,  together  with  the  indication  that  the  interfa¬ 
cial  region  is  the  cause  of  the  differences,  suggest  that  both  types  of  charge  must 
reach  the  interfacial  region  and  are  affected  by  it.  Thus,  at  least  a  portion  of 
both  types  of  charge  evidently  are  mobile  in  the  ZnS. 

The  increasing  flatband  shift  with  bias  is  probably  the  result  of  an  increasing 
escape  probability  as  the  field  is  increased.  This  escape  could  be  from  recombina¬ 
tion  (geminate,  columnar,  etc.)  or  trapping  of  the  radiation-induced  charges  in  the 
insulator.  There  is  no  indication  of  saturation  up  to  fields  of  ~0.75  MV/cm;  at 
higher  fields,  the  devices  burned  out.  The  exception  to  the  lack  of  saturation  is 
for  devices  with  the  ZnS-only  insulator  irradiated  under  negative  bias  (Figure  1 
or  4).  Here,  the  flatband  shift  is  essentially  bias  independent.  It  is  likely  that 

*The  applied  field  shown  in  these  figures  is  the  applied  bias  voltage  divided  by  the 
total  insulator  thickness.  For  dual-insulator  devices,  we  used  the  sum  of  the  two 
thicknesses.  This  assumes  that  the  field  is  the  same  in  both  insulators. 
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FIGURE  1.  RADIATION-INDUCED  FLATBAND 
SHIFT  IN  HgCdTe  MIS  CAPACITORS  WITH 
ZnS  INSULATOR. 


FIGURE  3.  RADIATION-INDUCED  FLATBAND 
SHIFT  IN  HgCdTe  MIS  CAPACITORS  WITH 
ZnS -PLUS- PLASMA-OX IDE  INSULATOR. 


FIGURE  2.  RADIATION-INDUCED  FLATBAND 
SHIFT  IN  HgCdTe  MIS  CAPACITORS  WITH 
Z nS - PLUS -ANODIC-OX IDE  INSULATOR. 


FIGURE  4.  COMPARISON  OF  RADIATION- 
INDUCED  FLATBAND  SHIFTS  IN  HgCdTe 
MIS  CAPACITORS. 
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a  different  mechanism  is  responsible  for  such  behavior,  because  the  escape 
probability  in  ZnS  should  be  independent  of  the  nature  or  method  of  formation  of 
the  native-oxide  layer  in  the  interfacial  region,  and  also  of  the  sign  of  the  bias. 

Based  on  the  above  findings,  we  can  begin  to  develop  a  model  to  explain  the 
observed  radiation-induced  effects.  Ionizing  radiation  produces  electron-hole 
pairs  in  the  ZnS,  and  at  least  some  of  these  pairs  are  separated  by  the  applied 

field  and  escape  the  bulk  region  of  the  ZnS  with  a  probability  that  depends  on  the 

magnitude  of  the  field  (but  which  is  independent  of  the  interfacial  region).  At 
least  a  fraction  of  both  the  electrons  and  holes  are  mobile  in  the  ZnS,  or  the  dif¬ 
ferences  observed  among  device  types  for  both  signs  of  applied  field  cannot  be  ex¬ 
plained.  (This  finding  is  qualitatively  different  from  that  for  SiOj  where  all  of 

the  holes  are  effectively  immobile  at  77°K.*)  The  mobile  charge  of  one  sign  is 

swept  toward  the  gate  from  which  it  exits  the  device.  The  mobile  charge  of  the 
other  sign  is  swept  toward  the  semiconductor  where  at  least  a  fraction  of  it  is 
trapped  in  the  interfacial  region.  A  trapping  probability  which  depends  on  the 
nature  of  the  interfacial  region  is  required  to  explain  the  different  flatband 
shifts  in  the  different  device  types.  When  a  native-oxide  layer  is  present,  the 
radiation-induced  holes  driven  toward  the  interface  under  positive  bias  evidently 
are  trapped  in  pre-existing  hole  traps  in  the  interfacial  layer.  The  anodic  oxide, 
having  more  of  these  traps  perhaps  because  of  its  greater  thickness,  exhibits  the 
largest  negative  flatband  shift,  even  though  its  ZnS  layer  is  much  thinner  and  thus 
would  generate  fewer  holes.  The  plasma  oxide  traps  only  a  slightly  larger  fraction 
of  the  radiation-induced  holes  than  does  the  interfacial  region  with  no  intentional 
native  oxide.  The  radiation-induced  electrons  driven  toward  the  interface  under 
negative  bias  are  trapped.  Some  of  these  trapping  sites  could  be  those  containing 
the  trapped  positive  charge  present  in  the  native  oxide  prior  to  irradiation  to 
which  the  electrons  would  be  coulombically  attracted.  The  anodic  oxide,  with  its 
greater  number  of  trapped  holes,  exhibits  only  slightly  larger  positive  flatband 
shifts  than  does  the  plasma  oxide.  Because  of  the  thinner  ZnS  layer  in  devices 
with  the  anodic-oxide  layer,  fewer  electrons  would  be  generated,  indicating  that 
the  anodic  oxide  traps  a  larger  fraction  of  electrons  as  well  as  holes  than  does 
the  plasma  oxide.  The  interfacial  region  with  no  intentional  native-oxide  layer 
appears  to  be  transparent  to  electrons.  The  flatband  shifts  produced  in  such  de¬ 
vices  are  likely  to  be  the  result  of  some  net  trapped  negative  charge  in  the  bulk 
of  the  ZnS.  If  so,  then  not  all  of  the  radiation-induced  charge  recombines  or  es¬ 
capes  the  ZnS.  The  trapping  probability  in  the  ZnS  for  this  case  appears  to  be 
bias  independent  because  the  flatband  shift  is  bias  independent  for  negative  bias 
greater  than  -0.2s  MV/cm. 

By  estimating  the  carrier-pair-creation  energy  in  ZnS  and  assuming  that  the 
charge  producing  the  flatband  shift  is  trapped  in  the  interfacial  region,  we  can 
calculate  the  fraction  of  the  radiation-induced  charge  trapped  at  that  location. 

The  carrier-pair-creation  energy  is  approximately  three  times  the  bandgap  (E^)  in 
most  materials,®  although  it  is  lower  than  this  in  a  few  materials.®  The  bandgap 

^J.R.  Srour  end  K.Y.  Chiu,  IEEE  Trans.  Nucl.  Sci.  24^  2140  (1977). 

SC.A.  Klein,  J.  Appl.  Phys.  39,  2029  (1968). 

®A.  Rothwarf,  J.  Appl.  Phys.  4j[,  752  (1973). 
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of  ZnS  is  3.6  eV,  so  the  maximum  pair-creation  energy  is  probably  ~ 11  eV/pair. 
Using  this  value,  we  calculate  that  a  flatband  shift  of  30  to  150  V  (depending  on 
the  ZnS-layer  thickness)  would  be  produced  by  a  total  dose  of  6.8  x  10^  rad(ZnS) 
if  all  of  the  radiation-induced  charge  of  one  sign  were  trapped  in  the  interfacial 
region  while  all  of  the  charge  of  the  opposite  sign  were  swept  out  of  the  ZnS. 

Thus,  the  observed  charge  trapped  in  the  interfacial  region  is  less  than  10*  of  the 
total  radiation-induced  charge.  The  fate  of  the  rest  of  the  charge  cannot  be  de¬ 
termined  from  the  present  measurements.  The  results  for  devices  with  no  inten¬ 
tional  native  oxide  indicate  that  at  least  some  of  the  charge  is  trapped  in  the 
ZnS  bulk,  but  whether  this  is  a  large  (i.e.,  large,  but  nearly  equal  numbers  of 
electrons  and  holes)  or  a  small  fraction  of  the  radiation-induced  charge  is  not 
known. 

IMPLICATIONS  FOR  ARRAY  STRUCTURES 

Irradiation-induced  charge  buildup  in  insulating  layers  can  have  a  number  of 
harmful  effects  on  the  properties  of  devices  in  HgCdTe  array  structures.  The  most 
straightforward  effect  occurs  when  the  insulator  is  a  gate  insulator.  In  this 
case,  the  flatband  voltage  shift  is  analogous  to  a  threshold  or  operating  voltage 
shift.  If  the  applied  device  biases  cannot  adapt  to  these  shifts,  the  device  oper¬ 
ating  conditions  will  no  longer  be  those  desired,  and  the  device  operation  could  be 
degraded.  For  example,  the  active  region  of  an  optical  charge  collection  device 
collects  and  integrates  the  charge  while  in  depletion.  If  irradiation  produced  a 
negative  flatband  shift  such  that  the  device  were  actually  in  accumulation  rather 
than  depletion  at  the  applied  bias,  it  would  no  longer  collect  the  charge. 

Another  effect  that  charge  buildup  can  produce  is  an  altered  surface  condition 
in  the  semiconductor.  Many  HgCdTe  optical  devices  exhibit  optimal  operation  (low¬ 
est  leakage  current  and  noise)  with  the  surface  at  approximately  flatband.  Irradi¬ 
ation-induced  charge  buildup  could  cause  the  HgCdTe  surface  to  accumulate  or  de¬ 
plete,  and  this  usually  has  a  deleterious  effect  on  device  operation. 

The  changed  HgCdTe  surface  condition  could  be  particularly  troublesome  when  it 
occurs  under  the  passivation  oxide  in  an  array  structure.  Conducting  paths  could 
be  induced  where  none  existed  prior  to  irradiation,  and  these  could  increase  the 
effective  areas  of  the  active  regions.  A  schematic  of  this  effect  for  an  n-on-p 
mesa- junction  structure  with  a  passivation  layer  between  the  junctions  is  shown  in 
Figure  5.  Trapped  positive  charge  would  invert  the  HgCdTe  surface  and  increase  the 
collection  area  of  the  devices.  This  was  the  mechanism  postulated  to  be  respon¬ 
sible  for  the  racuation-induced  changes  in  a  photodetector  array  with  this  struc- 

2 

ture.  One  difficulty  with  the  postulated  mechanism  is  that  the  array  passivation 
layer  was  ZnS  with  no  intentional  native  oxide.  In  the  present  study,  MIS  capaci¬ 
tors  with  this  insulator  exhibited  a  net  trapped  negative  charge  (Figure  1  or  4) 
rather  than  the  trapped  positive  charge  required  by  the  postulated  mechanism.  Per¬ 
haps  some  structural  difference  (e.g.,  the  capacitors  were  MIS  structures  while  the 
array  passivation  layer  had  no  metal  gate),  test  condition  difference  (e.g.,  the 
capacitor  gate  was  connected  to  the  substrate  for  2ero-bias  irradiation  while  the 
array  passivation  layer  was  floating),  or  nonuniform  charge  trapping  in  the  non- 
planar  array  structure  can  account  for  this  apparent  discrepancy.  Additional  work 
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FIGURE  5.  SCHEMATIC  REPRESENTATION  OF  SURFACE  INVERSION  IN 
A  HgCdTe  n-ON-p  PHOTOVOLTAIC  DETECTOR  ARRAY  STRUCTURE  AS  A 
RESULT  OF  TOTAL  DOSE  EXPOSURE. 


is  required  to  resolve  this  mechanistic  issue.  The  observed  effect  is  an  example 
of  a  problem  that  is  unique  to  array  structures  and  cannot  be  evaluated  by  examina¬ 
tion  of  discrete  devices. 

The  results  of  this  study  have  shown  that  exposure  to  ionizing  radiation 
causes  charge  buildup  in  the  common  insulator  structures  used  in  HgCdTe  arrays. 

Even  though  only  a  small  fraction  of  the  ionization-induced  charge  is  observed,  the 
resulting  flatband  shifts  are  significant,  with  measurable  shifts  produced  by  a 
total  dose  of  103  rad(ZnS)  or  less.  Most  of  the  observed  charge  buildup  appears  to 
occur  in  the  interfacial  region  where  it  would  be  particularly  effective  in  alter¬ 
ing  the  surface  condition  of  the  HgCdTe.  All  HgCdTe  array  structures  are  likely  to 
be  affected  in  one  way  or  another  by  charge  buildup  in  insulating  layers. 

The  amount  of  charge  buildup  differs  among  the  different  insulator  structures, 
and  this  could  lead  to  potential  device  hardening  techniques.  For  example,  the  in¬ 
sulator  structure  with  the  plasma  oxide  exhibits  the  least  flatband  shift  at  zero 
bias,  and  thus  would  be  the  logical  candidate  for  the  formation  of  radiation-toler¬ 
ant  passivation  layers.  On  the  other  hand,  for  a  radiation-tolerant  gate  insulator 
operated  at  negative  bias  (greater  than  -0.25  MV/cm),  the  candidate  insulator  is 
the  one  with  no  intentional  native  oxide.  In  addition  to  achieving  hardness  by 
insulator  selection,  additional  studies  to  determine  techniques  to  decrease  charge 
trapping  in  the  interf  tial  region  could  lead  to  more  radiation-tolerant  insulator 
structures  and  thus  to  harder  HgCdTe  array  structures. 
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Abstract 

Ionization-induced  charge  buildup  in  HgCdTe  MIS 
capacitors  has  been  examined.  Most  of  the  capacitors 
studied  contained  ZnS  as  the  insulator  and  had  different 
native  oxide  layers  between  the  HgCdTe  and  the  ZnS.  The 
net  charge  buildup  measured  in  these  devices  following 
77°K  irradiation  is  attributed  to  the  charge  trapped  in 
the  interfacial  region  (i.e.,  between  the  HgCdTe  and  the 
ZnS).  Either  electron  or  hole  trapping  may  occur  in  that 
region,  depending  on  the  sign  of  the  electric  field  ap¬ 
plied  during  irradiation.  This  finding  implies  that  at 
least  a  fraction  of  both  the  electrons  and  holes  pro¬ 
duced  by  ionization  are  mobile  in  ZnS  at  77°K.  The 
amount  of  charge  trapped  in  the  interfacial  region  de¬ 
pended  on  the  nature  of  that  region,  particularly  on  the 
thickness  of  the  native  oxide,  and  on  the  magnitude  of 
the  field  applied  during  irradiation.  Some  studies  of 
capacitors  with  a  CVD  Si02  insulator  were  also  per¬ 
formed.  The  results  indicated  that  both  holes  and  elec¬ 
trons  are  trapped  in  the  oxide  during  77°K  irradiation, 
witn  tne  resulting  net  trapped  cnarge  being  relatively 
small.  A  room-temperature  anneal  evidently  removes  at 
least  a  portion  of  the  trapped  electrons,  resulting  in 
increased  net  positive  trapped  charge. 


Introduction 

A  common  feature  of  HgCdTe  arrays  used  for  infrared 
detection  and  signal  processing  is  the  presence  of  in¬ 
sulating  layers,  either  as  gate  insulators,  passivation 
layers,  or  field  insulators.  When  most  insulators  are 
are  exposed  to  ionizing  radiation,  charge  buildup  oc¬ 
curs  which  can  affect  device  operation.  Results  of  re¬ 
cent  studies1'2  indicate  that  the  response  of  HgCdTe 
photodetectors  to  ionizing  radiation  is  dominated  by 
the  insulator  response.  This  appears  to  be  the  case  in 
pn-junction  devices,  where  the  insulator  is  the  passi¬ 
vation  layer  around  the  active  junction  area,  and  in  MIS 
devices,  where  the  gate  insulator  is  important.  Thus, 
all  HgCdTe  array  structures  appear  to  be  vulnerable  to 
ionization-induced  charge  buildup  in  insulating  layers, 
but  the  damage  processes  in  these  insulators  are  not  yet 
fully  understood. 

We  have  performed  an  investigation  of  the  mechanisms 
responsible  for  ionization-induced  damage  in  the  insu¬ 
lators  used  for  HgCdTe  devices.  The  principal  mechan¬ 
ism  investigated  was  charge  buildup  in  insulators.  The 
primary  goals  were  to  obtain  the  information  required 
to  account  for  radiation-induced  degradation  of  HgCdTe 
photodetector  array  structures  and  to  postulate  proce¬ 
dures  for  increasing  the  radiation  tolerance  of  such 
arrays. 


Experimental  Methods 

Useful  test  structures  for  studvinq  insulator  and 
interface  properties  are  MIS  capacitors  where  the  fixed 
charge  in  the  insulator  and  the  Interface  state  density 
can  be  obtained  using  conventional  analysis  of  C-V  mea¬ 
surements.  The  devices  used  in  this  study  were  supolied 
by  Santa  Barbara  Research  Center  and  were  fabricated  on 
HgCdTe  substrate  material  having  a  bandgap  of  0.2  eV  or 


greater  to  minimize  tunneling,  which  can  affect  C-V  mea¬ 
surements.  The  standard  insulator  used  for  HgCdTe  de¬ 
vices  is  ZnS.  For  the  present  devices,  ZnS  layers  of 
thicknesses  from  ~280  to  600  nm  were  prepared  by  ther¬ 
mal  evaporation.  Three  different  HgCdTe  surface  treat¬ 
ments  were  employed  before  deposition  of  the  ZnS.  In  one 
set  of  devices,  the  HgCdTe  was  etched  and  no  intentional 
native  oxide  was  grown  on  the  surface.  (However,  a  thin 
(  -1  nm)  native-oxide  layer  does  form  thermally.)  In  a 
second  set  of  devices,  a  70-nm  native-oxide  layer  was 
grown  anodically  using  a  KOH  solution  in  water  and  glyc¬ 
erine  using  a  process  similar  to  that  described  in  Ref¬ 
erence  3.  In  a  third  set  of  devices,  a  6-nm  native-oxide 
layer  was  grown  in  an  oxygen  plasma. 4  In  addition  to  de¬ 
vices  containing  a  ZnS  layer,  a  fourth  set  of  devices 
containing  a  110-nm  CVD  SiO^  layer  was  fabricated  and 
tested.  The  area  of  all  capacitors  tested  was  0.2  mm2. 

A  schematic  diagram  of  the  structure  of  the  four  types 
of  capacitors  is  shown  in  Figure  1. 
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Figure  1.  Schematic  of  the  structure  of  the  HgCdTe  MIS 

capacitors. 

Irradiations  were  performed  using  a  Co60  source  with 
devices  biased  during  exposure.  Devices  were  held  at 
77°K  during  irradiation  exposure  and  post-irradiation 
measurements,  with  no  sample  heating  between  exposure 
and  measurement.  The  effect  of  bias  alone  was  examined 
by  applying  a  bias  stress  to  the  samples  for  the  same 
time  as  the  irradiation  time  (10  minutes).  Charge  build¬ 
up  in  insulating  layers  was  monitored  by  determining  the 
the  flatband  voltage  shift  from  C-V  curves. 


Experimental  Results 

1) laJiitiOh  p+uOuttU  a  VlilttfnJ  -SI, Hi  it  *if  trf  the 
devices  tested.  An  example  of  the  results  is  shown  in 
Figure  2  for  a  device  with  a  plasma-oxide  layer.  In  all 
cases,  the  flatband  shift  was  linear  with  total  dose  up 
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Figure  2.  Effect  of  gamma  irradiation  on  a  HgCdTe  MIS 
capacitor  with  ZnS-plus-plasma-oxide  insulator. 

to  the  maximum  test  dose  of  3.6  x  10  rad(ZnS).  This 
result  is  illustrated  in  Figure  3  for  a  device  with  a 
n  -iittljf  irtullifr^  Th“rrf  r  .  ■  ji^-ri1  jn  hr'w  ;n  »- 
vices  can  be  made  at  any  (lower)  dose. 


Figure  3.  Flatband  shift  as  a  function  of  total  dose 
in  a  HgCdTe  HIS  capacitor  with  a  ZnS  Insulator. 


•Ionizing  doses  are  stated  in  rad(ZnS)  because  the  dom¬ 
inant  effects  are  produced  in  the  ZnS.  For  conversion 
purposes,  1  rad(ZnS)  *  0.96  rad(Si)  for  Co®0  gamma 
rays  (E  &  1.25  MeV). 


The  dependence  of  flatband  shift  on  applied  field1 
is  shown  in  Figures  4  through  6  for  the  three  device 
types  containing  ZnS.  (The  sign  convention  used  shows 
the  gate-to-substrate  bias.)  The  bias-stress  results 
(labeled  "no  irradiation")  are  also  shown  on  the  same 
figures.  For  devices  with  an  intentional  native-oxide 
layer  (Figures  5  and  6),  bias-induced  flatband  shifts 
were  so  small  that  no  allowance  for  them  was  necessary. 
For  devices  with  a  ZnS-only  insulator  (Figure  4),  bias- 
induced  flatband  shifts  were  not  negligible  in  compari¬ 
son  to  the  radiation-induced  flatband  shifts.  The  sense 
of  the  bias-induced  shifts  indicates  that  the  mechanism 
responsible  is  field-induced  injection  of  charge  from 
the  HgCdTe  into  the  ZnS,  followed  by  its  trapping  in  the 
ZnS.  This  process  occurs  for  either  sign  of  charge  de¬ 
pending  on  the  sign  of  the  bias.  To  allow  for  the  bias- 
induced  shifts,  the  bias-alone  data  were  subtracted  from 
the  irradiation  data  to  give  the  dashed  line  shown  in 
Figure  4,  which  assumes  that  the  two  effects  are  inde¬ 
pendent. 

To  facilitate  comparison  between  the  different  insu¬ 
lator  types,  radiation-induced  flatband  shifts  for  the 
three  insulator  structures  are  shown  in  Figure  7.  The 
sign  of  the  radiation-induced  shifts  indicates  that,  in 
all  cases  (except  at  very  low  applied  bias),  a  net  posi¬ 
tive  charge  remains  in  the  insulator  when  devices  are 
irradiated  under  positive  bias,  while  a  net  negative 
charge  remains  when  devices  are  Irradiated  under  nega¬ 
tive  bias.  The  amount  of  flatband  shift  differs  for  the 
n  ffv.Ufrt  ty,  «i ,  'm*  su**  tf#  l  J  l*  mw*  tj-ktee 

than  the  native  oxide  in  all  cases,  we  assume  that  all 
of  the  measured  flatband  shifts  are  produced  by  charge 
that  was  initially  created  in  the  ZnS.  Assuming  that  the 
structure  of  the  ZnS  is  the  same  in  all  device  types, 


Figure  4.  Radiation-induced  flatband  shift  in  HgCdTe 
MIS  capacitors  with  a  ZnS-only  insulator. 


^The  applied  field  shown  in  these  figures  is  the  applied 
bias  voltage  divided  by  the  total  Insulator  thickness. 
For  dual-insulator  devices,  we  used  the  sum  of  the 
two  thicknesses.  This  assumes  that  the  field  is  the 
same  in  both  Insulators. 
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Figure  5.  Radiation-induced  flatband  shift  in  HgCdTe 
MIS  capacitors  with  a  ZnS-plus-anodic-oxide 
insulator. 


Figure  6.  Radiation-induced  flatband  shift  in  HgCdTe 
MIS  capacitors  with  a  ZnS-plus-plasma-oxide 
insulator. 

the  physical  processes  occuring  in  the  ZnS  should  be 
qualitatively  the  same.  Quantitative  differences  could 
exist  because  of  the  different  ZnS  thicknesses.  How¬ 
ever,  devices  with  the  anodic-oxide  insulator  had  the 
thinnest  ZnS  layer,  yet  showed  the  largest  flatband 
shifts,  while  devices  with  the  ZnS-only  insulator  had 
the  thickest  ZnS  layer,  yet  exhibited  the  smallest  flat- 
band  shifts.  Tnis  behavior  indicates  that  the  thickness 
of  the  ZnS  layer  is  not  the  cause  of  the  observed  dif¬ 
ferences,  and  suggests  that  the  interfaciil  region  (in¬ 
cluding  the  native  oxide)  between  the  ZnS  and  the  HgCdTe 


Figure  7.  Comparison  of  radiation-induced  flatband 

shift  in  HgCdTe  MIS  capacitors. 

is  the  cause.  Further,  the  observation  that  differences 
exist  among  the  three  device  types  for  both  signs  of  bias 
in  Figure  7,  together  with  the  indication  that  the  inter¬ 
facial  region  is  the  cause  of  the  differences,  suggests 
that  both  electrons  and  holes  must  reach  the  interfacial 
region  and  are  affected  by  it.  Thus,  at  least  a  portion 
of  both  types  of  charge  evidently  are  mobile  in  the. ZnS. 

No  time -dependent  changes  were  found  following  irradia¬ 
tion,  which  means  that  all  of  the  charge  motion  has  oc¬ 
curred  in  a  time  less  than  the  time  frame  (  ~10  minutes) 
of  the  experiment. 

The  increasing  flatband  shift  with  bias  is  probably 
the  result  of  an  increasing  escape  probability  as  the 
field  is  increased.  This  escape  could  be  from  recombi¬ 
nation  (geminate,  columnar,  etc.)  or  trapping  of  the 
radiation-induced  charges  in  the  insulator.  There  is  no 
indication  of  saturation  up  to  fields  of  -0.75  MV/cm; 
at  higher  fields,  the  devices  burned  out.  The  exception 
to  the  lack  of  saturation  is  for  devices  with  the  ZnS- 
only  insulator  irradiated  under  negative  bias  (Figure  4 
or  7).  Here,  the  flatband  shift  is  essentially  bias  in¬ 
dependent.  It  is  likely  that  a  different  mechanism  Is 
responsible  for  such  behavior,  because  the  escape  proba¬ 
bility  in  the  ZnS  should  be  independent  of  the  nature  or 
method  of  formation  of  the  native-oxide  layer  in  the  in¬ 
terfacial  region,  and  also  of  the  bias  polarity. 

The  experiments  performed  on  devices  with  an  Si02 
gate  insulator  were  not  as  extensive  as  those  performed 
on  devices  with  a  ZnS  layer.  The  only  tests  performed 
were  with  devices  biased  at  zero  and  -0.65  MV/cm  during 
irradiation.  A  negative  flatband  shift  was  observed  in 
both  cases  'which  is  indicative  of  trapped  positive  charge. 
Conventional  SIO2  on  Si,  which  is  thermally  grown  (TG) 
rather  than  prepared  by  CVO  as  was  the  Si02  layer  on 
HgCdTe,  also  shows  a  negative  flatband  shift  when  irrad¬ 
iated  at  77*K.5  A  comparison  of  the  present  CVD  SiO? 
results  with  results  for  TG  Si02  on  Si  based  on  previous 
studies  is  shown  in  Figure  8.  In  both  cases,  the  insu¬ 
lator  thickness  was  1)0  nm.  The  flatband  shift  for  the 
CVD  Si O2  was  much  smaller  and  appears  to  be  bias  indepen¬ 
dent  when  compared  to  that  of  TG  SiO?.  The  mechanism 
responsible  for  the  negative  flatband  shift  in  TG  SiOj 
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Figure  8.  Radiation-induced  flatband  shift  in  capaci¬ 
tors  with  SiOg  insulators. 

is  that  all  of  the  positive  charge  is  trapped  in  place 
(except  at  high  applied  fields)  and  the  negative  charge 
is  mobile  at  low  temperatures.  The  most  probable  mech¬ 
anism  to  explain  the  CVD  Si02  results  is  that  much  of 
the  negative  charge  is  trapped  as  well,  resulting  in  a 
smaller  net  trapped  positive  charge.  Room-temperature 
annealing  caused  an  additional  negative  flatband  shift 
in  the  device  with  the  CVD  Si O2  insulator.  These  data 
are  also  shown  in  Figure  8.  Evidently,  at  least  a  por¬ 
tion  of  the  trapped  electrons  become  mobile  and  are  re¬ 
moved  from  the  oxide  at  room  temperature,  leaving  behind 
a  larger  net  positive  trapped  charge. 


discussion 

Based  on  the  above  findings,  we  suggest  a  qualita¬ 
tive  model  to  explain  the  observed  radiation-induced  ef¬ 
fects  in  devices  containing  ZnS.  Ionizing  radiation 
produces  electron-hole  pairs  in  the  ZnS,  and  at  least 
some  of  these  pairs  are  separated  by  the  applied  field 
and  escape  the  bulk  region  of  the  ZnS  with  a  probability 
that  depends  on  the  magnitude  of  the  field  (but  which  is 
Independent  of  the  interfacial  region).  At  least  a  frac¬ 
tion  of  both  the  electrons  and  holes  are  mobile  in  the 
ZnS,  or  the  differences  observed  among  device  types  for 
both  signs  of  applied  field  cannot  be  explained.  (This 
finding  is  qualitatively  different  from  that  for  Si O2 
where  all  of  the  holes  are  effectively  immobile  at 
77“K.->)  Mobile  charge  of  one  sign  is  swept  toward  the 
gate  electrode  and  is  removed  from  the  device.  Mobile 
charge  of  the  other  sign  is  swept  toward  the  semiconduc¬ 
tor  where  at  least  a  fraction  of  it  is  trapped  in  the 
interfacial  region.  A  trapping  orobability  which  depends 
on  the  nature  of  the  interfacial  region  is  required  to 
explain  the  lifferent  flatband  shifts  in  the  different 
device  type  .  When  a  native-oxide  layer  is  present,  the 
radiation-induced  holes  driven  toward  the  interface 
under  positive  bias  evidently  are  trapped  in  preexisting 
hole  traps  in  the  interfacial  layer.  The  anodic  oxide, 
having  more  of  these  traps  perhaps  beause  of  its  greater 


thickness,  exhibits  the  largest  negative  flatband  shifts, 
even  though  its  ZnS  layer  is  much  thinner  and  thus  would 
result  in  fewer  holes  being  generated.  The  plasma  oxide 
traps  only  a  slightly  larger  fraction  of  the  radiation- 
induced  holes  than  does  the  interfacial  region  with  no 
intentional  native  oxide.  The  radiation-induced  elec¬ 
trons  driven  toward  the  interface  under  negative  bias 
are  trapped.  Some  of  these  trapping  sites  could  be  those 
containing  the  trapped  positive  charge  present  in  the 
native  oxide  prior  to  irradiation  to  which  the  electrons 
would  be  coulombically  attracted.  The  anodic  oxide,  with 
its  greater  number  of  trapped  holes,  exhibits  only 
slightly  larger  positive  flatband  shifts  than  does  the 
plasma  oxide.  Because  of  the  thinner  ZnS  layer  in  de¬ 
vices  with  the  anodic-oxide  layer,  fewer  electrons  would 
be  generated,  indicating  that  the  anodic  oxide  traps  a 
larger  fraction  of  electrons  as  well  as  holes  than  does 
the  plasma  oxide.  The  interfacial  region  with  no  inten¬ 
tional  native-oxide  layer  appears  to  be  transparent  to 
electrons.  The  flatband  shifts  produced  in  such  devices 
are  likely  to  be  the  result  of  some  net  trapped  negative 
charge  in  the  bulk  of  the  ZnS.  If  so,  then  not  all  of 
the  radiation-induced  charge  recombines  or  escapes  the 
ZnS.  The  trapping  probability  in  the  ZnS  for  this  case 
fc  T*  "da.  >*c*ui.*  eh*  fl'tband  shift 

is  bias  independent  for  negative  bias  greater  than  -0.25 
MV/cm. 

During  our  experiments,  flatband  voltages  never  ex¬ 
ceeded  ±  10V.  (Flatband  voltage  shifts  were  limited  to 
less  than  ±5V.)  At  no  time  was  there  any  indication 
that  all  charge  traps  in  the  interfacial  regions  were 
full.  From  this  information,  we  calculate  that  greater 
than  3  x  103  charges  can  be  trapped  in  the  interfacial 
regions.  This  equates  to  a  concentration  of  ~2  x  10'' 
charges/cm-*  in  the  anodic  oxide  and  -2  x  10lB  charges/ 
cm'  in  the  plasma  oxide.  Thus,  these  are  the  minimum 
trap  densities  in  the  oxides.  At  the  ZnS-HgCdTe  inter¬ 
face  ( i . e. ,  devices  with  no  intentional  native  oxide), 
the  hole  trap  density  is  at  least  10^2  traps/cmz,  while 
the  electron  trap  density  is  apparently  zero. 

By  estimating  the  carrier-pair-creation  energy  in 
ZnS  and  assuming  that  the  charge  producing  the  flatband 
shift  is  trapped  in  the  interfacial  region,  we  can  cal¬ 
culate  the  fraction  of  the  radiation-induced  charge  trap¬ 
ped  at  that  location.  The  carrier-pal r-creat ion  energy 
is  approximately  three  times  the  bandgap  (Eg)  in  most 
materials,®  although  it  is  lower  than  this  in  some  mate¬ 
rials.  The  bandgap  of  ZnS  is  3.6  eV,  so  the  maximum 
pair-creation  energy  is  probably  ~11  eV/pair.  Using 
this  value,  we  calculate  that  a  flatband  shift  of  30  to 
150V  (depending  on  the  ZnS-Iayer  thickness)  would  be 
produced  by  a  total  dose  of  6.B  x  103  rad(ZnS)  if  all 
of  the  radiation-induced  charge  of  one  sign  were  trapped 
in  the  interfacial  region  while  all  of  the  charge  of  the 
opposite  sign  were  swept  out  of  the  ZnS.  Thus,  the  ob¬ 
served  charge  trapped  in  the  interfacial  region  is  less 
than  lot;  of  the  total  radiation-induced  charge,  The  fate 
of  the  rest  of  the  charge  cannot  be  determined  from  the 
present  measurements.  The  results  for  devices  with  no 
intentional  native  oxide  indicate  that  at  least  some  of 
the  charge  is  trapped  in  the  ZnS  bulk,  but  whether  this 
is  a  large  (i.e.,  large,  but  nearly  equal  numbers  of 
electrons  and  holes)  or  a  small  fraction  of  the  radiation- 
induced  charge  is  not  known. 

One  feature  of  the  present  total-dose  results  was 
that  the  flatband  shift  was  not  zero  for  zero  applied 
field  (i.e.,  the  curves  in  Figure  7  do  not  pass  through 
the  origin).  The  probaule  cause  of  this  offset  is  that  a 
built-in  field  exists  in  the  ZnS  at  zero  bias.  This  field 
could  be  due  to  work  function  differences  at  the  inter¬ 
faces.  In  devices  with  no  native  oxide,  the  internal  field 
evidently  is  negative  since  a  positive  applied  field  is 
required  to  overcome  it.  Adding  the  plasma  oxide  appears 
to  decrease  the  internal  negative  field,  and  adding  the 
anodic  oxide  changes  the  sign  of  the  internal  field. 
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Although  the  changes  in  internal  field  as  a  function  of 
native  oxide  thickness  appear  to  be  systematic  {thicker 
native  oxide  produces  a  more  positive  field),  the  cause 
is  not  understood.  The  observed  effects  are  not  the  re¬ 
sult  of  the  positive  charge  trapped  in  the  native  ox¬ 
ide,  because  this  would  produce  negative  fields,  which 
is  opposite  to  what  is  observed. 


Concluding  Remarks 

The  results  of  this  study  have  shown  that  a  major 
damage  mechanism  in  HgCdTe  MIS  devices  bombarded  with 
ionizing  radiation  is  the  production  of  charge  in  the 
ZnS  insulator  and  its  subsequent  capture  in  the  inter¬ 
facial  region  between  the  ZnS  and  the  HgCdTe.  Although 
only  a  small  fraction  of  the  radiation-induced  charge 
is  actually  observed,  the  effects  of  this  charge  can  be 
significant  even  at  a  fairly  low  total  dose.  The  trap¬ 
ped  charge  could  clearly  affect  devices  with  an  MIS 
structure  by  shifting  the  threshold  voltage.  It  is 
quite  likely  that  the  effects  observed  in  pn-junction 
devices*  are  also  the  result  of  charge  production  in 
the  ZnS  passivation  layer  because  flatband  shifts  were 
observed  here  even  with  no  applied  field.  The  mechan¬ 
ism  responsible  may  be  a  change  in  the  surface  condi¬ 
tion  (i.e.,  accumulation  or  inversion)  of  the  HgCdTe. 

The  effects  observed  in  this  investigation  depended 
on  the  insulato"-  field  applied  during  irradiation  and 
on  the  nature  of  the  interfacial  region  between  the  ZnS 
and  the  HgCdTe.  Thus,  the  effect  of  radiation  on  HgCdTe 
devices  depends  on  insulator  fabrication  procedures 
even  at  zero  applied  field,  which  would  be  the  condi¬ 
tion  under  which  insulators  on  many  devices  would  oper¬ 
ate.  Thus,  the  damage  mechanism  observed  here  is  ex¬ 
pected  to  dominate  the  radiation  response  and  be  a 


first-order  failure  mechanism  in  all  present-day  HgCdTe 
array  structures. 

The  amount  of  charge  buildup  differs  among  the  dif¬ 
ferent  insulator  structures,  and  this  could  lead  to 
potential  device  hardening  techniques.  Tor  example,  the 
insulator  structure  with  the  plasma  oxide  exhibits  the 
least  flatband  shift  at  zero  bias,  and  thus  would  be 
the  logical  candidate  for  the  formation  of  radiation- 
tolerant  passivation  layers.  On  the  other  hand,  for  a 
radiation-tolerant  gate  insulator  operated  at  negative 
u4as  VgA  atrf  tr.ai  tW.TTi),  Ihc  carfiSett  WrWd 

tor  is  the  one  with  no  intentional  native  oxide.  In  ad¬ 
dition  to  achieving  hardness  by  insulator  selection, 
additional  studies  to  determine  techniques  to  decrease 
charge  trapping  in  the  interfacial  region  could  lead  to 
more  radiation  tolerant  insulator  structures  and  thus 
to  harder  HgCdTe  array  structures. 
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M.  Frank 

ATTN 

P.  Lathrop 

ATTN 

D.  Toomb 

Bendix  Corp 

ATTN 

SV/B7 11/70 

ATTN 

E.  Meeder 

ATTN 

0.  Huffman 

Boeing  Aerospace  Co 

Aerospace  Corp 

ATTN 

MS-B1-36,  W.  Doherty 

ATTN 

R.  Crolius 

ATTN 

MS-2R-D0,  E.  Smith 

ATTN 

J.  Stoll 

ATTN 

MS-2R-00,  C.  Rosenberg 

ATTN 

H.  Phillips 

ATTN 

0.  Mu  1  key 

ATTN 

J.  Hiesner 

ATTN 

MS-2R-00,  I.  Arimura 

ATTN 

J.  Reinheimer 

ATTN 

C.  Dixon 

ATTN 

A.  Carlan 

ATTN 

MS-B1-36,  P.  Blakely 

ATTN 

W.  Kolasinski ,  MS/259 

ATTN 

MS-2R-00,  A.  Johnston 

ATTN 

V.  Josephson,  MS-4-933 

ATTN 

R.  Slaughter 

Boeing  Co 

ATTN 

D.  Fresh 

ATTN 

8K-3B 

ATTN 

I.  Garfunkel 

ATTN 

H.  Wicklein 

ATTN 

S.  Bower 

ATTN 

R.  Caldwell 

ATTN 

C.  Huang 

ATTN 

D.  Egelkrout 

ATTN 

W.  Crane,  A2/1083 

ATTN 

P,  Buchman 

Booz,  Allen  and  Hamilton,  Inc 

ATTN 

G.  Gilley 

ATTN 

R.  Chrisner 

ATTN 

D.  Schmuck 

ATTN 

B.  Blake 

Burroughs  Corp 

ATTN 

Product  Evaluation  Lab 

Aerospace  Industries  Assoc  of  America,  Inc 

ATTN 

S.  Siegel 

Cal  ifornia 

Institute  of  Technology 

ATTN 

D.  Nichols,  T-11BD 

Ampex  Corp 

ATTN 

F.  Grunthaner 

ATTN 

J.  Smith 

ATTN 

J.  Bryden 

ATTN 

D.  Knutson 

ATTN 

W.  Price,  MS-B3-122 

ATTN 

R.  Covey 

Analytic  Svcs,  Inc 

ATTN 

P.  Robinson 

ATTN 

P.  Szymanski 

J.  O'Sullivan 

ATTN 

A.  Shumka 

ATTN 

ATTN 

K.  Martin 

ATTN 

A.  Shostak 

ATTN 

W.  Scott 

AVCO  Systems  Oiv 

Charles  Stark  Draper  Lab,  Inc 

ATTN 

W.  Eroding 

ATTN 

A.  Schutz 

ATTN 

D.  Fann 

ATTN 

W.  Callender 

ATTN 

C.  Davis 

ATTN 

D.  Gold 

ATTN 

D.  Shrader 

ATTN 

J.  Boyle 

ATTN 

A.  Freeman 

Battelle  Memorial  Institute 

ATTN 

N.  Tibbetts 

ATTN 

R.  Thatcher 

ATTN 

R.  Ledger 

ATTN 

R.  Haltmaler 

BDM  Corp 

ATTN 

Tech  Library 

ATTN 

S.  Meth 

ATTN 

P.  Greiff 

ATTN 

C.  Stickley 

ATTN 

R.  Bedingfield 

BDM  Corp 

Cincinnati 

Electronics  Corp 

ATTN 

R.  Antinone 

ATTN 

L.  Hammond 

ATTN 

Marketing 

ATTN 

C.  Stump 

ATTN 

0.  Wunsch 
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DEPARTMENT  OF  OEFENSE  CONTRACTORS  (Continued) 


Computer  Sciences  Corp 
ATTN:  A.  Schiff 

Control  Data  Corp 

ATTN:  0.  Newberry,  BRR  142 
ATTN;  T.  Frey 

University  of  Denver 

ATTN:  Sec  Officer  for  F.  Venditti 

Develco,  Inc 

ATTN:  G.  Hoffman 


General  Electric  Co 

ATTN:  G.  Gati,  MD-E184 

General  Electric  Co 

ATTN:  C.  Hewison 
ATTN:  J.  Gibson 
ATTN:  D.  Cole 

General  Electric  Co 
ATTN:  D.  Pepin 

General  Research  Corp 


ATTN: 

E.  Steele 

Oi kewood 

ATTN: 

R.  Hill 

ATTN:  Tech  Library  for  L.  Davis 

ATTN: 

Tech  Info  Ofc 

ATTN: 

A.  Hunt 

E-Systems,  Inc 

ATTN:  K.  Reis 

E-Systems,  Inc 

ATTN:  Oiv  Library 

Eaton  Corp 

ATTN:  A.  Anthony 
ATTN:  R.  Bryant 

University  of  Florida 
ATTN:  H.  Sisler 

Ford  Aerospace  5  Communications  Corp 
ATTN:  H.  Linder 
ATTN:  J.  Davison 
ATTN:  Tech  Info  Svcs 

Ford  Aerospace  J  Communications  Corp 
ATTN:  E.  Hahn 

Franklin  Institute 

ATTN:  R.  Thompson 

Garrett  Corp 

ATTN:  H.  Weil 

General  Dynamics  Corp 

ATTN:  R.  Fields,  MZ  2839 
ATTN:  0.  Wood 


Goodyear  Aerospace  Corp 

ATTN:  Sec  Con  Station 

Grumman  Aerospace  Corp 
ATTN;  J.  Rogers 


Harris  Corp 
ATTN: 
ATTN: 
ATTN: 

Harris  Corp 
ATTN: 
ATTN: 
ATTN: 
ATTN: 
ATTN: 
ATTN: 
ATTN: 
ATTN: 


W.  Abare 
C.  Davis 
E.  Yost 


Mgr  Linear  Engrg 
J.  Cornell 

C.  Anderson 
J.  Schroeder 

B.  Gingerich,  MS-51-120 

D.  Williams,  MS-51-75 
T.  Sanders,  MS-51-121 
Mgr  Bipolar  Digital  Eng 


Hateltine  Corp 

ATTN:  C.  Meinen 
ATTN:  J.  Okrent 


Honeywell,  Inc 
ATTN:  R 
ATTN:  F 
ATTN:  D 
ATTN:  J 


General  Electric  Co 
ATTN:  R.  Ben 
ATTN:  Tech  L 
ATTN:  J.  Pal 
ATTN:  Tech  I 
ATTN:  W.  Pat 
ATTN:  J.  And 
ATTN:  J.  Ped 
ATTN:  R.  Cas 
ATTN:  0.  Tas 


R.  Gumm 
F.  Hampton 

D.  Neilsen,  MN  14-3015 
J.  Moylan 


R.  Benedict 

Tech  Library 

J.  Palchefsky,  Jr 

Tech  Info  Ctr  for  L.  Chasen 

W.  Patterson 

J .  Andrews 

J.  Peden 

R.  Casey 

0.  Tasc> 


Honeywell,  Inc 
ATTN:  J 
AT1N:  R 
ATTN:  C 
ATTN:  J 
ATTN:  H 
ATTN:  M 


J.  Zawacki 
R.  Reinecke 
C.  Cerulli 
J.  Schafer 
H.  Noble 
MS  725-5 


General  Electric  Co 


Honeywell,  Inc 

ATTN:  Tech  Library 


ATTN: 

L.  Hauge 

Honeywell,  Inc 

ATTN: 

B.  Flaherty 

Attention  Sec  MN11-2D6D 

ATTN: 

G.  Bender 

ATTN:  L.  Lavoie 

ATTN: 

J.  Reidl 
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Honeywell,  Inc 

ATTN:  D.  Herold,  MS-MN,  17-2334 
ATTN:  D.  Lamb,  MS-MN,  17-2334 
ATTN:  R.  Belt,  MS-MN,  17-2334 

Hughes  Aircraft  Co 

ATTN:  K.  Walker 
ATTN:  D.  Binder 
ATTN:  R.  McGowan 
ATTN:  CTDC,  6/ El 10 

Hughes  Aircraft  Co 

ATTN:  E.  Smith,  MS  V347 
ATTN:  E.  kubo 
ATTN:  W.  Scott,  S32/C332 
ATTN:  A.  Narevsky,  S32/C332 
ATTN:  D.  Shumake 

Hughes  Aircraft  Co 

ATTN:  R.  Henderson 

Hughes  Aircraft  Co 

ATTN:  MS-A2408,  J.  Ha'l 
ATTN:  P.  Coppen 

IBM  Corp 

ATTN:  Mono  Memory  Systems 
ATTN:  Electromagnetic  Compatibility 
ATTN:  H.  Mathers 
ATTN:  T.  Martin 

IBM  Corp 

ATTN:  A.  Edenfeld 

ATTN:  N.  Haddad 

ATTN:  W.  Doughten 

ATTN:  H.  Kotecha 

ATTN:  W.  Henley 

ATTN:  MS  110-036,  F.  Tietze 

ATTN:  L.  Rockett,  MS  110-020 

ATTN:  0.  Spencer 

ATTN:  S.  Saretto 

1 1 T  Research  Institute 
ATTN:  I.  Mindel 
ATTN:  R.  Sutkowski 

Institute  for  Defense  Analyses 
ATTN:  Tech  Info  Svcs 

International  Tel  S  Telegraph  Corp 
ATTN:  Dept  608 
ATTN:  A.  Richardson 

IRT  Corp 

ATTN;  J.  Harrity 
ATTN:  R.  Mertz 
ATTN:  N.  Rudie 
ATTN:  R.  Judge 
ATTN:  M.  Rose 
ATTN:  Systems  Effects  Oiv 
ATTN:  Physics  Div 
ATTN:  MDC 


JAYCOR 

ATTN;  R.  Stahl 
ATTN:  L.  Scott 
ATTN:  T.  Flanagan 
ATTN:  J.  Azarewicz 
ATTN:  M.  Treadaway 
ATTN:  R.  Berger 

JAYCOR 

ATTN:  R.  Sullivan 
ATTN:  E.  Alcaraz 

JAYCOR 

ATTN:  R.  Poll 

Johns  Hopkins  University 
ATTN:  R.  Maurer 
ATTN;  P.  Partridge 

Kaman  Sciences  Corp 

ATTN:  Dir  Science  &  Technology  Div 

ATTN:  C.  Baker 

ATTN:  W.  Rich 

ATTN:  N.  Beauchamp 

ATTN:  J.  Erskine 

kaman  Tempo 

ATTN:  W.  McNamara 
ATTN:  DAS  I AC 
ATTN:  R.  Rutherford 

kaman  Tempo 

ATTN:  W.  Alfonte 
ATTN:  DAS  I AC 

Litton  Systems,  Inc 

ATTN:  F.  Motter 
ATTN:  G.  Maddox 
ATTN:  J.  Retzler 

Lockheed  Missiles  &  Space  Co,  Inc 
ATTN:  F.  Junga,  S2/54-202 
ATTN:  J.  Smith 
ATTN:  Reports  Library 
ATTN:  J.  Crowley 

Lockheed  Missiles  &  Space  Co,  Inc 
ATTN:  k.  Greenough 
ATTN:  B.  ktmura 
ATTN:  L.  Rossi 

ATTN:  S.  Tiamuty,  Dept  81-74/154 
ATTN:  D.  Wolfhard 
ATTN:  J.  Cayot,  Dept  81-63 
ATTN:  J.  Lee 

ATTN:  Dr  G.  Lum,  Dept  81-63 
ATTN:  P.  Bene 
ATTN:  E.  Hessee 

M . I . T .  Lincoln  Lab 

ATTN:  P.  McKenzie 
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Magnavox  Advanced  Products  &  Sys  Co 
ATTN:  W.  Hagemeier 

Magnavox  Govt  A  Indus  Electronics  Co 
ATTN:  W.  Richeson 

Martin  Marietta  Corp 

ATTN:  S.  Bennett 
ATTN:  H.  Cates 
ATTN:  W.  Janocko 
ATTN:  W.  Brockett 
ATTN:  T1C/MP-3D 
ATTN:  R.  Yokomoto 
ATTN:  J.  Ward 
ATTN:  J.  Tanke 
ATTN:  R.  Gaynor 
ATTN:  MP-163,  W.  Bruce 
ATTN:  P.  Fender 
ATTN:  MP-163,  N.  Redmond 

Martin  Marietta  Denver  Aerospace 
ATTN:  M.  Shumaker 
ATTN:  D-6074,  G.  Freyer 
ATTN:  E.  Carter 
ATTN:  P.  Kase 
ATTN:  Research  Library 
ATTN:  MS-06074,  M.  Polzella 

McDonnell  Oouglas  Corp 

ATTN:  M.  Stitch,  Dept  E003 

ATTN:  Library 

ATTN:  R.  Kloster,  Oept  E451 

ATTN:  T.  Ender,  33/6/618 

ATTN:  0.  Dohm 

rtiTN:  A.  Munie 

McDonnell  Douglas  Corp 
ATTN:  P.  Bretch 
ATTN:  0.  Fitzgerald 
ATTN:  M.  Onoda 
ATTN:  P.  Albrecht 
ATTN:  M.  Ralsten 
ATTN:  J.  Holmgrem 
ATTN:  J.  Imai 
ATTN.  R.  Lothringer 

McDonnell  Douglas  Corp 

ATTN:  Tech  Library 

Mission  Research  Corp 


ATTN: 

C. 

Longmire 

ATTN: 

M. 

Van  Blaric 

ion  Researc 

Corp 

ATTN: 

R. 

Pease 

ATTN: 

R. 

Turfler 

ATTN: 

D. 

A1 exander 

ATTN: 

0. 

Merewether 

Mission  Research  Corp 

ATTN:  B.  Passenheim 
ATTN:  J.  Raymond 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

Mission  Research  Corp 
ATTN:  R.  Curry 
ATTN:  W.  Ware 
ATTN:  J.  Lubell 

Mitre  Corp 

ATTN:  M.  Fitzgerald 

Motorola,  Inc 

ATTN:  A.  Christensen 

Motorola,  Inc 

ATTN:  C.  Lund 

National  Academy  of  Sciences 

ATTN:  National  Materials  Advisory  Board 

National  Semiconductor  Corp 
ATTN:  A.  London 
ATTN:  F.  Jones 
ATTN:  J.  Martin 

New  Technology,  Inc 
ATTN:  0.  Divis 

Norden  Systems,  Inc 
ATTN:  0.  Longo 
ATTN:  Tech  Library 

Northrop  Corp 

ATTN:  P.  Eisenberg 

2  cy  ATTN:  A.  Kalma 

Z  cy  ATTN;  Z.  Shanfield 

Z  cy  ATTN:  A.  Bahraman 
Z  cy  ATTN:  J.  Srour 
2  cy  ATTN:  S.  Othmer 
Z  cy  ATTN:  R.  Hartman 
2  cy  ATTN:  M,  Hopkins 

Northrop  Corp 

ATTN:  L.  Apodaca 
ATTN:  P.  Besser 
ATTN:  0.  Strobel 
ATTN:  P.  Eisenberg 
ATTN:  T.  Jackson 
ATTN:  P.  Gardner 
ATTN:  S.  Stewart 
ATTN:  E.  King,  C3323/WC 

Pac ific-Sierra  Rsch  Corp 

ATTN:  H.  Brode,  Chairman  SAGE 

Palisades  Inst  Rsch  Svcs,  Inc 
ATTN:  Secretary 

Physics  International  Co 
ATTN:  Div  6000 
ATTN:  J.  Shea 

ATTN:  J.  Huntington 

R  S  0  Associates 

ATTN:  W.  Karzas 
ATTN:  C.  Rogers 
ATTN:  P.  Haas 
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DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 


Rand  Corp 

ATTN:  C.  Crain 


Raytheon  Co 


G.  Joshi 
T.  Wein 
J.  Ciccio 


Raytheon  Co 

ATTN:  H.  Flescher 
ATTN:  A.  Van  Doren 


RCA  Corp 
ATTN: 

RCA  Corp 

ATTN: 

ATTN: 

ATTN: 

ATTN: 

ATTN: 

RCA  Corp 

ATTN: 

RCA  Corp 

ATTN: 

ATTN: 

ATTN: 

RCA  Corp 

ATTN: 

ATTN: 

ATTN: 

ATTN: 


V.  Mancino 


L.  Minich 
0.  O'Connor 
G.  Hughes 
R.  Smeltzer 
L.  Napoli 


R.  Killion 


W.  Allen 
L.  Oebacker 
E.  Schmitt 


J.  Sault2 
E.  Van  Keuren 
W.  Heagerty 
R.  Magyarics 


Rensselaer  Polytechnic  Institute 
ATTN:  R.  Ryan 

Research  Triangle  Institute 
ATTN:  M.  Simons 

Rockwell  International  Corp 
ATTN:  J.  Bell 

ATTN:  GA50,  TIC/L,  G.  Green 

ATTN:  A.  Rovell 

ATTN:  V.  Strahan 

ATTN:  J.  Blandford 

ATTN:  C.  Kleiner 

ATTN:  K.  Hull 

ATTN:  J.  Pickel ,  Code  03I-BBQ1 
ATTN:  V.  De  Martino 
ATTN:  R.  Pane  holy 
ATTN:  V.  Michel 


Rockwell  Corp 

ATTN:  TIC,  BADE 
ATTN:  T,  Tates 

Sanders  Associates,  Inc 
ATTN:  L.  Brodeur 
ATTN:  M.  Aitel 

Science  Applications,  Inc 
ATTN:  V.  Verbinski 
ATTN:  J.  Spratt 
ATTN:  V.  Orphan 
ATTN:  R.  Fitzwilson 
ATTN:  J.  Beyster 
ATTN:  J.  Naber 
ATTN:  0.  Strobel 
ATTN:  0.  Millward 
ATTN:  0.  Long 
ATTN;  L.  Scott 

Science  Appl ications,  Inc 
ATTN:  C.  Cheek 
ATTN:  J.  Swirczynski 

Science  Applications,  Inc 
ATTN:  W.  Chadsey 
ATTN:  J.  Wallace 

Science  Applications,  Inc 
ATTN:  D.  Stribling 

Scientific  Research  Assoc,  Inc 
ATTN:  H.  Grubin 

Singer  Co 

ATTN:  R.  Spiegel 
ATTN:  J.  Brinkman 
ATTN:  J.  Laduca 
ATTN:  Tech  Info  Ctr 

Sperry  Corp 

ATTN:  Engrg  Lab 

Sperry  Corp 

ATTN:  J.  Inda 

Sperry  Flight  Systems 
ATTN:  0.  Schow 

Sperry  Rand  Corp 

ATTN:  R.  Viola 
ATTN:  C.  Craig 
ATTN:  F.  Scaravagl ione 
ATTN;  P.  Maraffino 


Rockwell  International  Corp 
ATTN:  0.  Stevens 
ATTN:  TIC,  0/41-092,  AJ01 

Rockwell  International  Corp 
ATTN:  A.  Langenfield 
ATTN:  TIC,  106-216 


SRI  International 

ATTN:  A.  Whitson 

SRI  International 

ATTN:  A.  Padgett 

Sundstrand  Corp 

ATTN:  Research  Dept 
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DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 
Sylvanla  Systems  Group  TRW  Electronics  &  Defense  Sector 


ATTN 

L. 

Paupl  is 

ATTN 

A.  Witteles,  MS  Rl/2144 

ATTN 

C. 

Thornhill 

ATTN 

D,  Clement 

ATTN 

L. 

Blaisdel  1 

ATTN 

J.  Bell 

ATTN 

W. 

Dunnet 

ATTN 

W.  Willis 

ATTN 

H,  Holloway 

Sylvania  Systems  Group 

ATTN 

Vulnerability  4  Hardness  Lab 

ATTN 

C. 

Ramsbottom 

ATTN 

P.  Gardner 

ATTN 

H. 

U1  Iman 

ATTN 

W.  Rowan 

ATTN 

P. 

Fredrickson 

ATTN 

P.  Guilfoyle 

ATTN 

H  4  V  Group 

ATTN 

H.  Hennecke 

ATTN 

P.  Reid,  MS  R6/2541 

Strategic  Systems  Div 

ATTN 

H.  Volmerange,  Rl/1126 

ATTN 

J. 

Waldron 

ATTN 

Tech  Info  Ctr 

ATTN 

R.  Kingsland 

Systron-Donner  Corp 

ATTN 

M.  Ash 

ATTN 

J. 

Indel icato 

ATTN 

F.  Friedt 

2  cy  ATTN 

R.  Plebuch 

Teledyne  Brown  Engrg 

2  cy  ATTN 

0.  Adams 

ATTN 

J. 

Me Swain 

ATTN 

D. 

Guice 

Vought  Corp 

ATTN 

T. 

Henderson 

ATTN 

Library 

ATTN 

Tech  Data  Ctr 

Teledyne  Systems  Company 

ATTN 

R.  Tomme 

ATTN 

R. 

Suhrke 

Westinghouse  Electric  Corp 

Texas  Instruments,  Inc 

ATT" 

MS  330,  D.  Grimes 

ATTN 

R. 

McGrath 

ATTN 

N,  Bluzer 

ATTN 

R. 

Stehl in 

ATTN 

E.  Vitek,  MS  320D 

ATTN 

T. 

Cheek,  MS  3143 

ATTN 

L.  McPherson 

ATTN 

F. 

Poblenz,  MS 

3143 

ATTN 

H.  Kalapaca ,  MS  3330 

ATTN 

E. 

Jeffrey,  MS 

961 

ATTN 

J.  Cricchi 

ATTN 

D. 

Manus 

ATTN 

MS  3330 

ATTN 

R. 

Carroll,  MS 

3143 

Westinghouse  Electric  Corp 

TRW  Electronics 

4  Defense  Sector 

ATTN 

S.  Wood 

ATTN 

J. 

Gorman 

ATTN 

F. 

Fay 

ATTN 

c. 

Blasnek 

ATTN 

R. 

Kitter 
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